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PREFACE 


This  report  summarizes  the  work  performed  at  The  Pennsylvania  State 
University  during  the  period  from  1  January  1982  to  29  January  1986,  under  Grant 
No.  AFOSR-82-0013.  The  Principal  Investigator  was  Carlo  G.  Pantano,  Faculty 
Associates  were  Jack  J.  Mecholsky  and  Richard  C.  Bradt,  Research  Associate  was 
Cheryl  A.  Houser,  and  Graduate-Student  Research  Assistants  included  J.  Chesson, 
A.  Tesar,  D.  Coon,  M.  Murtagh,  A.  Gonzalez  and  A.  Phelps. 

A  major  objective  of  this  program  was  to  evaluate  the  slow  crack  growth 
behavior  of  a  systematic  compositional  series  of  glasses  and  to  correlate  these 
data  with  independent  studies  of  the  surface  chemistry.  This  objective  was,  in 
fact,  accomplished  using  the  glass  system  Na20* xAI ^ ( 3-x )SiC>2 •  This 
endeavor  led  to  new  models  for  understanding  stress  corrosion  susceptibility, 
chemical  effects  during  crack  growth  and  the  static  fatigue  limit.  After  a 
general  discussion  of  chemical  effects  in  fracture  (Chapter  I),  these  new  models 
are  described  in  Chapters  II,  III,  IV  and  V  of  this  report.  The  development  of 
these  new  concepts  would  not  have  been  possible  without  the  careful  measurement 
and  documentation  of  trends  throughout  a  family  of  glasses.  A  comparable  study 
of  f luorozi rconate  glasses  was  initiated  during  the  third  year  of  this  program; 
the  results  obtained  to  date  are  highlighted  in  Chapter  VI.  Finally,  Chapter 
VII  presents  the  first  reported  study  of  the  mechanical  behavior  -  and  its 
dependence  upon  environment  -  of  a  silica  sol/gel  material. 
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SURFACE  CHEMISTRY  IN  RELATION  TO  THE  STRENGTH 
AND  FRACTURE  OF  SILICATE  GLASSES 


by 


Carlo  G.  Pantano 

Department  of  Materials  Science  and  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


in:  Strength  of  Inorganic  Glass,  (C.  R.  Kirkjian,  Editor,  Plenum 
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SURFACE  CHEMISTRY  IN  RELATION  TO  THE  STRENGTH  AND  FRACTURE  OF  SILICATE  GLASSES 


by 


Carlo  G.  Pantano 

Department  of  Materials  Science  &  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Introduction 

The  strength  and  fracture  behavior  of  silica  and  silicate  glasses  can  be 
influenced  by  surface  chemistry  in  many  different  ways.  It  may  be  appropriate 
to  start  this  discussion  by  attempting  to  identify  the  most  important 
relationships. 


(1)  Surface  chemical  phenomena  can  result  in  compositional  and  structural 
changes  at  glass  surfaces  during  high  temperature  fabrication 
processes.  In  the  case  of  silica  glass,  this  is  limited  to  chemical 
adsorption  of  water  and  other  gaseous  species,  as  well  as  surface 
segregation  and/or  volatization  of  water  and  impurities  from  the  bulk  of 
the  glass.  In  the  case  of  multicomponent  glasses,  more  dramatic 
composition  changes  may  occur  due  to  the  mobility  and  reactivity  of  the 
modifier  species.  The  extent  of  chemical  and  structural  modification  of 
the  surface  is  a  complex  function  of  the  bulk  composition  of  the  glass 
and  the  thermal  and  atmospheric  history  of  the  glass  surface. 

Obviously,  any  high  temperature  chemical  or  structural  alteration  of  the 
surface  will  influence  subsequent  adsorption  and  surface  chemical 
behavior  in  a  manner  which  cannot  be  predicted  simply  from  a  knowledge 
of  the  bulk  glass  composition.  It  will  be  shown,  for  example,  that  for 
a  given  glass  composition,  some  high  temperature  fabrication  processes 
enhance  the  hydration  resistance  of  the  surface  while  others  degrade  it. 

(2)  The  chemisorption  of  gaseous  species  at  surfaces  and  crack-tips 
under  residual  or  applied  stresses  can  contribute  to  the  rupture  of 
interatomic  bonds.  This  phenomena  is,  of  course,  fundamental  to  slow 
crack  growth  and  will  not  be  treated  in  any  detail  here.  However,  it  is 
worth  pointing  out  that  slow  crack  growth  does  exhibit  a  composition 
dependence.  Thus,  the  strength  behavior  of  real  glass  surfaces  will  be 
influenced  by  any  in-depth  modification  of  the  surface  during  high 
temperature  fabrication. 
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(3)  Chemical  reactions  between  the  glass  surface  and  adsorbed  species  can 
lead  to  broken  bonds,  hydration  films  and  other  structural 
heterogenities  at  the  surface  under  ambient  conditions.  These  reactions 
will  undoubtedly  be  localized  at  high  energy  surface  sites,  especially 
those  associated  with  strained  bonds.  And  while  all  the  bonds  at  the 
surface  may  be  strained  to  some  extent,  the  strain  may  be  accentuated  at 
regions  associated  with  bulk  topological  features,  phase  separated 
regions,  crystallites  or  impurities.  The  distinction  between  this 
phenomena  and  slow  crack  growth  is  perhaps  philosophical  since  both 
involve  bond  breaking  chemical  reactions  at  high  energy  surface  sites. 
However,  it  has  a  special  significance  in  the  sense  that  it  may  be  a 
realistic  source  of  flaws,  or  flaw-precursors,  and  may  in  some  instances 
occur  at  high  energy  sites  not  associated  with  strained  bonds  per  se. 

The  alkali  species  in  multicomponent  glasses,  for  example,  will 
preferentially  react  with  water  and  carbon  dioxide.  The  reaction  and 
distribution  of  surface  reaction  products  and  corrosion  pits  will  be 
non-uniform  to  an  extent  determined  by  the  composition  and  thermal 
history  of  the  glass. 

In  general,  one  expects  these  kinds  of  interactions  to  be  auto-catalytic 
in  the  sense  that  any  subsequent  chemical  reactions  will  also  be 
localized  at  these  sites.  For  example,  the  reaction  between  a  strained 
siloxane  linkage  (=Si-0-Si=)  and  water  not  only  leads  to  a  broken  bond, 
but  the  silanol  groups  (=Si-OH)  so-produced  will  enhance  the  physical 
adsorption  of  water  at  that  same  site.  This  relatively  mobile  water  is 
then  available  to  react  with  other  siloxane  linkages  in  that  region 
whether  they  are  strained  or  not. 

(4)  The  mechanically  induced  'flaws'  at  glass  surfaces  are  high  energy  sites 
where  the  adsorption  and/or  condensation  of  mobile  water  can  result  in 
corrosion  reactions.  That  is,  the  chemistry  of  the  condensed  water  can 
be  modified  dramatically  by  the  dissociation  of  surface  silanol  groups 
and  ion-exchange  reactions  with  modifier  species  due  to  the  high  ratio 
of  surface  area  to  solution  volume.  Thus,  the  geometry  of  the  flaw  may 
be  modified  by  gel  formation  and  swelling,  non-uniform  dissolution,  and 
etching. 

(5)  The  physical  adsorption  of  mobile  water  or  other  gaseous  species  at 
pristine  surfaces,  at  corapositionally-modif ied  surfaces,  on  surface 
reaction  products,  within  pre-existing  flaws  or  at  the  interface  between 
the  glass  surface  and  organic  coatings  provides  a  reservoir  of  reactants 
for  chemical  interactions  upon  the  application  of  stress.  It  is  we  1 1 
known  from  a  wide  variety  of  fracture  mechanics  studies  that  the 
transport  of  water  or  other  reactants  can  limit  slow  crack  growth.  It 
is  possible  that  the  most  significant  influence  of  hydrated  surface 
layers  -  especially  within  microcracks  -  is  their  capacity  to  'store' 
physically  adsorbed  water  until  a  stress  is  applied.  This  'reservoir 
of  reactants'  becomes  especially  significant  when  it  exists  at  the 
interface  between  a  glass  fiber  and  organic  coating. 


Obviously,  one  cannot  be  too  dogmatic  about  these  classifications  since 
they  are  all  dependent  upon  the  same  fundamental  phenomena.  However,  these 
classifications  do  facilitate  the  design  of  critical  experiments,  and  perhaps 
most  importantly,  the  interpretation  of  strength  behavior  for  commercial  glass 


products. 

It  becomes  immediately  evident  that  melt  surfaces  and  fracture  surfaces  are 
most  appropriate  for  experimental  investigation.  The  melt  surfaces  of  glass  are 
among  the  strongest  surfaces  available  and  the  most  prevalent  in  practice.  Most 
importantly,  though,  it  is  the  strength  behavior  of  these  surfaces  which  is  most 
dependent  upon  chemical  interactions.  Understanding  the  way  in  which  these 
surfaces  evolve  at  high  temperature,  during  cooling,  and  during  exposure  to  and 
service  in  the  ambient  atmosphere,  is  extremely  important  but  not  well 
established.  Fracture  surfaces,  on  the  other  hand,  are  most  useful  as  'model' 
surfaces  for  studying  the  evolution  of  melt  surfaces  and  the  behavior  of  cracks 
and  crack-tips  in  the  absence  stress. 

In  this  paper,  we  will  not  consider  the  stress  enhanced  surface  chemical 
reactions  associated  with  slow  crack  growth  phenomena  since  these  are  addressed 
independently  in  these  proceedings.  Rather,  the  hydroxy lation  of  silica  and 
silicate  glasses,  the  characteristics  of  melt  surfaces,  and  the  chemistry  and 
structure  of  clean  fracture  surfaces  will  be  emphasized.  Finally  the 
relationships  between  surface  chemistry  and  fracture  behavior  will  be 
considered.  Although  this  discussion  will  refer  to  the  use  of  a  wide  variety  of 
analytical  techniques  for  studying  glass  surfaces,  the  details  associated  with 
each  of  these  methods  will  not  be  specifically  addressed.  The  reader  is 
referred  to  a  number  of  review  articles  where  the  theory,  instrumentation, 
applicability  and  shortcomings  of  these  methods  are  presented  with  regard  to 
glass  surface  studies^--’. 
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HYDROXYL  GROUPS  AMD  HYDRATION  AT  GLASS  SURFACES 

Silica  Glass 

The  chemisorption  of  water  -commonly  termed  hydroxylation  or  hydroiysis- 
dominates  the  surface  chemical  behavior  of  clean  or  pristine  silicate  glasses. 
There  are  a  number  of  sites  where  the  chemisorption  of  water  may  occur;  viz 


=  Si+  +  H20— >  zSi-OH  +  H+  Cl) 

=  Si-0~  +  H20— aSi-OH  +  0H_  (2) 

= Si-O-Si  =+  H20— 2(=  Si-OH)  (3) 

=  Si-OX  +  H20—~ Si-OH  +  XOH  where  X  «  Li,  Na,  K,  Ca,  Ba,  Pb,  etc.  (4) 


These  reactions  are  fundamental  to  the  surface  chemistry  and  fracture  behavior 
of  silicate  glasses  for  two  important  reasons: 

(1)  regardless  of  how  or  where  the  chemisorption  takes  place,  a  silanol 
group  -SiOH-  is  produced,  where  physical  adsorption  and  clustering  of 
water  and  other  polar  molecules  will  be  localized,  and 

(2)  hydroxylation  of  siloxane  linkages  (HSi-0-Si=)  breaks  those  bonds  which 
comprise  the  structural  backbone  of  silicate  glasses. 

There  has  been  considerable  investigation  of  the  behavior  of  surface  silanol 

groups,  but  vanishingly  little  investigation  of  their  creation  at  high  or  low 

temperatures  on  clean  surfaces. 

Much  of  the  work  concerning  the  behavior  of  silanol  groups  has  been 
confined  to  high  surface  area  silicas.  These  include  silica  gel,  porous  'Vycor' 
glass  and  silica  'fume'  (e.g.  Cab-O-Sil).  In  addition  to  their  technological 
application  as  dessicants  and  packing  for  chromatographic  columns,  the  primary 
rationale  for  the  extensive  study  of  high  surface  area  silicas  is  their 
adaptability  to  infrared  absorption  analysis  .  This  technique  not  only 
provides  a  qualitative  and  quantitative  measure  of  the  hydroxyl  groups  and 
molecular  water  at  the  surface,  but  also  provides  an  indication  of  their 
configuration. 
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Extensive  reviews  are  already  available  concerning  infrared  analysis  and 
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surface  chemistry  of  high  surface  area  silicas  .  Here,  only  the  most 
important  observations  are  presented.  The  major  objective  of  this  section  is  to 
bring  out  those  observations  that  are  most  relevant  to  the  strength  behavior  of 
silica  and  silicate  glasses. 

It  is  well  established  that  silica  surfaces  at  ambient  conditions  are 
covered  to  varying  degrees  with  physically  adsorbed  molecular  water  (often 
termed  water  of  hydration).  However,  it  is  the  underlying  chemisorbed  water 
(often  termed  hydroxyl  or  silanol  groups)  which  is  of  most  concern.  There  are 
two  silanol  configurations  which  can  be  readily  distinguished : vicinal  and 
isolated;  viz, 


vicinal 


H  „--H 
O'''  0 


isolated 


Si  Si 

/|\  /1\ 


/°\ 

Si  Si 

/|\  /1\ 


/°\ 

Si  Si 

A\  /1\ 


1 

f 

fi 

K 
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The  vicinal  hydroxyl  groups  are  hydrogen  bonded  Co  one  another  due  Co  their 
close  proximity  across  pores,  cracks  and  breaks  in  the  silica  ring  structures. 
And  in  some  instances,  it  is  possible  to  distinguish  between  weakly-hydrogen 
bonded  and  strongly-hydrogen  bonded  vicinal  hydroxyl  groups.  The  point  is  that 
the  concentration  and  configuration  of  the  silanol  groups  determine  the  capacity 
for  physical  adsorption  of  molecular  water;  i.e.  physical  adsorption  of 
molecular  water  occurs  predominantly  at  the  vicinal  silanol  groups: 


H 

O 


/ 


/|\ 


/°\  I  I 

Si  Si  Si  Si 

/|\  /|\  /|\  /|\ 


The  molecular  water  of  hydration,  the  vicinal  silanol  groups  and  the  isolated 
silanol  groups  are  readily  identified  with  infrared  absorption  spectroscopy  as 
shown  in  Figure  1. 
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The  molecular  water  is  physically  adsorbed  to  the  silica  surface  through 
hydrogen  bonds  with  the  vicinal  silanol  groups.  Of  course,  multilayer 
adsorption  and  clustering  of  a  molecular  water  can  occur  through  intermolecular 
hydrogen  bonding  of  the  water  molecules  themselves.  Nevertheless,  upon  vacuum 
evacuation  or  moderate  heating,  the  molecular  water  is  readily  eliminated  from 
the  surface  (see  figure  1).  On  the  other  hand,  removal  of  the  underlying 
chemisorbed  water  is  more  difficult.  Temperatures  of  the  order  500C  and  higher 
are  required  to  affect  the  chemisorbed  water  (see  figure  1).  Under  these 
conditions,  the  vicinal  hydroxyl  groups  can  condense  to  liberate  molecular  water 
and  thereby  produce  a  siloxane  linkage  (sSi-O-Sii);  viz 

I 

H  ~  H 
O-  "  '  O 

|  i  heat  0  +H,0 

I  I—*"/  \ 

Si  Si  Si  Si 

A\  A\  A\  /|\ 

Obviously,  the  strongly  hydrogen  bonded  groups  condense  more  readily  that  the 
weakly  hydrogen  bonded  groups.  Most  significant,  though,  is  the  fact  that  the 
siloxane  linkages  formed  by  the  condensation  can  be  'strained1  upon  cooling  to 
ambient  temperatures.  The  amount  of  strain  depends  upon  the  original 
configuration  of  the  vicinal  hydroxyls  and  the  maximum  temperature  of 
dehydration.  Naturally,  there  is  no  mechanism  for  the  removal  of  the  isolated 
silanol  groups  during  a  high  temperature  dehydration  (see  Figure  1).  That  is, 
at  temperatures  of  the  order  850C-900C,  the  concentration  of  surface  silanol 
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groups  is  less  Chan  1  per  100  ft2  and  chese  remaining  silanol  groups  are 
of  the  freely  vibrating  kind. 

Since  the  silanol  groups  -  particularly  the  vicinal  variety  -  are 
responsible  for  the  physical  adsorption  of  molecular  water,  it  is  not  surprising 
that  silicas  dehydrated  at  elevated  temperatures  can  be  extremely  hydrophobic. 

Of  course,  they  will  re-hydrate  upon  exposure  to  water  at  a  rate  dependent  upon 
the  water  vapor  partial  pressure,  the  amount  of  ‘strain’  in  the  siloxane 
linkages  and  the  configuration  of  the  siloxane  linkages  (see  Figure  1).  It  is 
shown  below  that  the  application  of  an  external  stress  accelerates  Che 
re-hydration. 

The  density  of  the  silanol  groups  found  on  a  silica  glass  surface  obviously 
depends  upon  the  thermal  and  atmospheric  history  of  the  sample  but  is  typically 
4-8  per  100  ft2.  It  is  worth  pointing  out  that  if  one  calculates  the 
minimum  number  of  terminal  sites  (sSi+  or  =SiO-)  required  for  the  very 
existence  of  a  surface  on  silica  glass,  one  expects  about  2-3  silanols  per 
100ft2  to  be  formed  via  reactions  (1)  and  (2).  The  additional  silanols 
often  found  at  silica  surfaces  may  be  due  to  roughness  or  porosity,  but  are 
probably  due  to  the  hydrolysis  of  siloxane  bridges  (reaction  3).  Unfortunately, 
a  measure  of  the  initial  and/or  equilibrium  silanol  density  at  a  melt  or 
fracture  surface  of  silica  glass  has  never  been  reported  under  any  conditions. 

Since  most  of  the  experimental  data  used  to  formulate  the  conclusions 
summarized  above  were  obtained  using  high  surface  area  silicas,  one  must 
question  the  significance  of  these  observations  with  regard  to  the  surfaces  cf 
bulk  or  massive  glasses.  The  internal  microporosity  of  the  silica  gel  and 
porous  glass  are  of  some  concern.  And  although  the  silica  fume  -notably 
Cab-O-Sil-  is  a  compact  particulate  silica  which  contains  no  internal  porosity, 
its  preparation  for  infared  absorption  analysis  requires  that  it  be  mixed  with 


KBr  and/or  be  pressed  inco  a  powder  compact  with  extreme  pressure.  Thus,  the 
studies  of  Bershtein^ and  co-workers  are  extremely  significant.  Here, 
attenuated  total  reflectance  (ATR)  infrared  spectroscopy  was  used  to  monitor 
hydroxyl  groups  at  the  surface  of  massive  silica  glass  specimens.  In  contrast 
to  the  use  of  high  surface  area  microporous  glasses,  silica  gels  or  powder 
compacts  (Figure  2a),  the  multiple  ATR  element  provides  a  more  realistic  glass 
surface  for  study  (Figure  2b).  Here,  silica  glass  with  a  very  low  bulk  water 
must  be  used.  However,  any  adsorption  of  water  on  the  surface  is  readily 
detected  due  to  absorption  of  the  infrared  radiation  at  each  internal 
reflection.  Moreover,  it  is  possible  to  hydrate  these  massive  surfaces  under  an 
applied  load,  as  well  as,  to  subject  them  to  strength  measurements.  This  work 
has  already  been  reviewed  in  the  literature  and  need  not  be  covered  in  detail 
here^;  however,  a  summary  of  the  most  significant  findings  is  appropriate. 

These  ATR  studies  showed  that  the  surfaces  of  fused  silica  -after  hydration 
at  temperatures  of  the  order  70°C  to  600°C  -  contained  both  physically 
adsorbed  molecular  water  and  silanol  groups.  And,  the  concentration  of  silanol 
groups  was  much  higher  when  the  hydration  was  carried  out  under  stress. 

Although  the  spectral  resolution  of  the  ATR  data  is  inferior  to  the  direct 
transmission  data  obtained  with  powder  compacts,  the  published  spectra  do 
suggest  the  presence  of  both  vicinal  and  freely  vibrating  hydroxyl  groups. 

Perhaps  of  greater  significance  was  the  conclusion  of  Bershtein,  et 
al.^  that  the  hydrolysis  of  siloxane  linkages  occurs  preferentially  at 
weakened  ^si-O-Si^  sites  in  the  glass  and  that  the  concentration  of  these 
'weakened'  sites  set  a  limit  upon  the  number  of  sites  which  could  be 
hydroxylated.  They  found,  for  example,  that  increasing  the  time  of  hydration 
increased  the  thickness  of  the  hydrated  layer  but  not  the  concentration  of 
hydroxvls  within  it,  typically  3-5  percent  by  weight.  They  showed  further  that 
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hydrolysis  of  the  weakened  bonds  created  'passages'  for  the  intrusion  of 
molecular  water.  The  'passages'  were  observed  even  when  the  hydration  was 
carried  out  in  the  absence  of  an  applied  stress.  Using  mercury  porosimetry, 
they  found  that  these  'passages'  had  widths  less  than  1  urn,  with  the  majority 
possessing  widths  less  than  600  8.  The  'weakened'  bonds  responsible  for 
the  preferential  hydrolysis,  intrusion  of  molecular  water  and  corresponding 
creation  of  defects,  were  suggested  to  be  associated  with  sub-microscopic 
'ordered'  regions  in  the  glass  structure.  These  investigators  suggested  that 
the  presence  of  impurities  -particularly  bulk  water-could  reduce  the  stresses  at 
the  interfaces  or  boundaries  of  these  ordered  regions,  whereas  in  high  purity 
silica  the  non-uniform  or  localized  weakening  of  siloxane  bridges  could  be 
accentuated. 

Finally,  the  Bershtein  group10*11  related  the  strength  of  the  hydrated 
silicas  to  the  measured  distribution  of  hydroxyls  in  the  surface.  They  showed  a 
distinct  reduction  in  strength  due  to  the  non-uniform  hydrolysis  of  siloxane 
bonds,  and  an  increase  in  strength  upon  acid  etching  of  the  hydrated  layers. 
Nonetheless,  they  also  showed  that  by  carrying  out  the  hydrolysis  at  elevated 
temperatures  and  in  the  presence  of  sufficient  water,  the  strength  could  in  fact 
be  improved  with  hydration.  This  they  attributed  to  the  elimination  of 
'weakened'  siloxane  linkages,  the  creation  of  a  'microplastic'  surface  layer, 
and  consequently  a  more  uniform  distribution  of  applied  stress. 

The  work  on  high  surface  area  silicas  and  compact  silica  plates  emphasizes 
Che  importance  of  strained  surface  bonds,  the  configuration  of  silanols  and 
siloxane  linkages,  and  of  course,  the  presence  of  an  externally  applied  stress. 
Recently,  Dunken  and  Hoffman  4  have  attempted  to  shed  light  on  these 
concepts  using  theoretical  'cluster'  calculations.  They  essentially  calculate 
the  probability  of  a  water  molecule  disrupting  siloxane  linkages  as  a  function 


of  Che  geometry  and  composition  of  the  cluster.  Although  this  work,  too,  has 
been  summarized  in  a  recent  review^,  some  of  the  more  important 
observations  are  included  here.  Tor  example,  calculations  on  a 
(Si^O^H^)  cluster  showed  that  the  simultaneous  interaction  of  the 
hydroxyl  portion  of  a  water  molecule  with  a  silicon  atom,  and  the  protonic  end 
of  the  water  molecule  with  the  siloxane  bridge,  can  be  more  favorable  than  an 
interaction  between  the  water  molecule  and  a  silanol.  Moreover,  the 
calculations  showed  that  siloxane  bridges  between  tetrahedra  already  containing 
a  silanol  are  most  susceptible  to  hydrolysis.  With  regard  to  the  geometry  of 
the  cluster,  the  calculations  suggested  that  hydrolysis  of  a  siloxane  linkage 
with  an  angle  <115°  is  more  favorable  than  hydrolysis  of  a  stretched 
siloxane  linkage! 

Some  direct  evidence  for  the  presence  of  strained  siloxane  bonds,  and  their 
preferential  hydrolysis,  has  been  provided  by  Raman  spectroscopy.  The  Raman 
spectra  for  vitreous  silica^,  and  for  dehydrated  silica  gels^’1^, 
exhibit  bands  at  605  cm  *  and  495  cm  *  which  are  assigned  to  the 
vibration  of  three-membered  and  four-membered  planar  rings,  respectively.  The 
three-membered  ring  is  of  particular  significance  because  it  represents  a  high 
energy  configuration  consisting  of  strained  siloxane  bonds.  Krol  and  van  Lierop 
examined  dehydrated  silica  gel  and  found  that  the  existence  of  these  three-  and 
four-membered  rings  correlated  with  the  presence  of  freely  vibrating  silanol 
groups;  i.e.  the  degree  of  freedom  provided  by  the  surface  stabilizes  these  high 
energy  planar-ring  configurations.  These  authors  also  examined  the  re-hydration 
of  silica  gel  with  Raman  spectroscopy^,  Upon  exposing  the  dehydrated  gel 
to  the  ambient  atmosphere,  they  observed  a  complete  disappearance  of  the 
605  cm  ^  band  with  little  or  no  corresponding  attenuation  of  the  band  at  495  cm 
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Thus,  it  can  be  concluded  chat  the  highly  strained  siloxane  bonds  in  the 
three-merabered  rings  were  preferentially  hydrolyzed  by  water. 

Obviously,  the  configuration  of  silica  tetrahedra  at  the  tip  of  a  highly 
stressed  crack  will  also  have  a  strong  influence  on  chemical  reactivity  and 
therefore  upon  slow  crack  growth  phenomena.  Although  there  is  no  way  to 

i 

describe  or  measure  this  characteristic  at  the  present  time,  Wiederhorn*7  ] 

I 

showed  that  the  interaction  between  stressed  crack  tips  and  water  was  very  ! 

i 

different  in  normal  vs  anomalous  glasses.  Since  a  major  difference  between  | 

normal  and  anomalous  glasses  is  the  way  in  which  =Si-0-Sis  linkages  are 
distorted  under  stress,  the  configuration  dependence  of  hydrolysis  reactions  is 
further  supported.  | 

In  general,  then,  we  are  at  the  point  where  additional  understanding  of  the 
hydroxy lation  of  silica  glass  will  require  more  extensive  characterization  of 
the  structure  of  glass  -  in  the  bulk  and  at  the  free  surface.  The  schematic  in  | 

Figure  3a  emphasizes  in  a  more  realistic  way  the  role  of  the  bulk  glass 
structure  and  ring  statistics  in  determining  the  number  of  strained  siloxane 

bridges  at  the  surface,  and  consequently,  the  total  number  and  configuration  of  j 

silanol  groups.  Unfortunately,  the  ability  to  further  develop  these  ideas  is 

limited  due  to  the  lack  of  experimental  techniques  which  are  sensitive  to  the 

strain  and  configuration  associated  with  siloxane  bridges  and  hydroxyl  groups, 

as  well  as  to  the  network  topology.  Nonetheless,  it  is  probably  worth 

emphasizing  the  point  that  the  hydroxyl  groups  are  responsible  for  the 

adsorption  of  water,  and  in  a  sense,  represent  sites  of  preferential  attack 

where  water  molecules  gain  access  to  the  structure.  In  fact,  the  molecular 

water  will  'cluster'  at  these  regions  due  to  hydrogen  bonding  interactions 

between  the  water  molecules  themselves.  In  principle,  of  course,  the  absence  of 

silanol  (or  other  polar)  surface  sites  would  greatly  inhibit  the  adsorption  of 


water. 


The  formation  of  silanol  groups  on  silica  glass  is  usually  limited  to  the 
outermost  monolayer.  In  contrast,  the  hydroxy lation  and  hydration  of  glasses 

containing  modifier  species  -particularly  alkali  ions-  occurs  in 

1  18  19 

depth  »  *  .  The  hydrogen  depth  profiles  in  Figure  4  show  quite  clearly 
that  the  ability  to  hydrate  glass  in  water  is  enhanced  for  alkali  containing 
glasses.  In  general,  the  extent  of  hydrogen  incorporation  correlates  directly 
with  the  alkali  content.  Of  course,  these  profiles  provide  little  indication  of 
the  molecular  form  of  the  hydrogen,  the  configuration  of  the  silanols,  or  the 
uniformity  of  hydration. 

In-depth  hydration  of  alkali  containing  glasses  also  occurs  due  to  the 
•  •  1  20  21 

water  vapor  in  the  ambient  atmosphere  ’  ’  .  Figure  5  shows  schematically 

the  concentration  profiles  observed  for  alkali  and  hydrogen  at  alkali-silicate 
glass  surfaces.  These  reactions  were  proposed  on  the  basis  of  AES  (Auger 
electron  spectroscopy),  XPS  (x-ray  photoelectron  spectroscopy)  and  SIMS 
(secondary  ion  mass  spectroscopy )  studies  of  fracture  surfaces,  sputter-cleaned 
surfaces,  and  a  variety  of  glass  commercial  surfaces,  all  after  exposure  to  the 
ambient  atmosphere. 

It  is  quite  clear  that  the  phya  cal  adsorption  of  molecular  water  is  more 
extensive  on  these  multicomponent  silicate  glasses  than  it  is  on  pure  silica 
glass.  The  reasons  for  this  are  demonstrated  schematically  in  Figures  3b  and 
3c.  In  addition  to  the  silanol  groups  (which  form  due  to  dangling  Si+  and 
SiO  sites),  the  alkali  ions  are  also  capable  of  adsorbing  molecular  water. 

The  adsorbed  molecular  water  can  react  with  the  non-bridging  oxygens  in  the 
sub-surface  to  liberate  the  alkali  ions  and  create  free  hydroxyl  ions.  The 
alkali  ions  are  generally  observed  to  diffuse  out  of  the  hydrated  surface 
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whereas  the  hydroxyl  ions  will  attack  siloxane  bridges  to  further  depolymerize 
the  silicate  network.  The  enhanced  attack  of  siloxane  bonds  under  these 
alkaline  conditions  produces  many  more  silanol  groups.  Thus,  extensive  in-depth 
hydration,  and  perhaps  dissolution,  of  the  surface  is  expected.  The  alkali  ions 
liberated  from  the  structure  will  accumulate  at  the  glass-atmosphere  interface 
where  they  react  with  CO^  to  form  carbonate  species  (see  Figure  5).  The 
carbonate  species  and  surface  alkalinity  are  not  observed  for  pure  silica  glass. 

The  existence  of  modified  surface  layers  on  multicomponent  glasses  as  shown 
in  Figure  5  -  after  even  short  exposure  to  the  ambient  atmosphere  -  seems  to  be 

a  general  phenomena.  It  is  in  this  regard  that  the  observations  of  Stavrinidis 
22 

and  Holloway  are  quite  interesting.  They  observed  a  narrow  region 
(~12nm)  of  low  reflectivity  at  the  interface  created  by  healing  cracks  in  a 
soda-lime-silica  glass.  Since  this  interface  was  formed  only  when  the  crack 
surfaces  were  first  exposed  to  a  humid  atmosphere,  it  is  undoubtedly  due  to  a 
hydrated  surface  film.  Interestingly,  these  hydrated  crack  surfaces  were  able 
to  be  healed. 

Also  significant  in  the  hydration  of  alkali  containing  glasses  is  the 
microporosity  which  can  be  generated  within  the  modified  surface  layer.  In  some 
cases,  the  intrusion  of  molecular  water,  attack  of  siloxane  linkages  by  free 
hydroxyl  and  molecular  water,  and  condensation  of  some  silanol  groups,  can  all 

combine  to  transform  and  swell  the  surface  into  a  microporous  'gel'. 

23  24 

Sewell  ,  and  also  Martin,  Akmc  and  Oh  ,  showed  the  presence  of 

microporosity  on  the  hydrated  surfaces  of  E-glass  and  soda- 1 ime-s i 1 ica  glass 

fibers.  Nonetheless,  it  is  important  to  emphasize  that  the  hydroxy lat ion  of 

alkali  containing  glasses  need  not  result  in  a  'gel'  layer  formation.  Indeed, 

the  alkali  ions  in  the  surface  and  sub-surface  are  replaced  by  hydrogen-bearing 

species  to  form  a  'hydrogen'  glass.  However,  it  is  the  initial  composition,  and 


particularly  the  network  structure,  of  the  surface  which  determines  the  extent 
to  which  it  hydrates  and  is  transformed  into  a  swollen  'gel'  layer  via  the 
'solvation'  of  internal  silanols  by  molecular  water. 

The  dependence  of  hydration  and  'gel'  formation  upon  the  subtle 
characteristics  of  the  glass  network  structure  is  especially  evident  in  the 
binary  alkali-silicates.  One  can  compare  surface  film  formation  on  binary 
alkali  silicate  glasses  whose  molar  alkali  concentration  -  and  therefore 
non-bridging  oxygen  density  -  are  identical.  The  hydrated  surface  films  are 
more  stable  in  the  order  lithia-silica  >  soda-silica  >  potassia  silical 

• 

This  is  due  to  the  lower  partial  molar  volume  for  the  silica  network  in  the 
glasses  with  smaller  alkali  ions.  The  reduced  volume  of  pore  space  in  the 
leached  lithia  silica  glasses  inhibits  the  intrusion  of  water  and  breakdown  of 
the  'silica-rich'  surface  film  which  forms  upon  preferential  release  of  the 
alkali  component.  In  the  potassia  silica  glasses,  on  the  other  hand,  a  large 
pore  volume  makes  the  intrusion  of  water  and  the  hydrolysis  of  siloxane  bridges 
easier;  e.g.  the  ionic  radius  of  K  is  approximately  equal  to  that  of 
H^0+.  This  enhances  the  production  of  a  soluble,  colloidal  silica 
'gel'  film  at  the  surface. 

In  general,  then,  one  expects  the  hydration  of  multicomponent  silicate 

glasses  to  occur  in  stages  as  shown  schematically  in  Figure  6.  The  relative 

rates  of  'hydrogen-glass'  formation  vs  'gel'  formation  will  depend  upon  the 

glass  composition  and  structure,  the  solution  pH  and  the  temperature.  The 

simultaneous  creation  of  mechanical  stresses,  and  the  extent  to  which  they  are 

transmitted  to  the  sub-su..jce,  will  be  dependent  upon  the  details.  However, 

the  spontaneous  cracking  of  glass  due  to  hydration  and  surface  film  formation 
n  5 

has  been  reported  . 
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MELT  SURFACES 

It  has  been  pointed  out  that  the  hydration  of  glass  surfaces  is  very  much 
dependent  upon  composition  and  structure.  During  the  formation  and  cooling  of  a 
melt  surface,  reorientation  of  the  silicate  structure,  segregation  of 
impurities,  volatilization  of  modifiers,  water,  and  other  impurities  and 
chemisorption  phenomena  determine  the  final  characteristics  of  the  glass 
surface.  It  is  generally  observed  that  melt  surfaces  have  better  hydration 
resistance  and  lower  chemical  reactivity  than  ground,  polished  or  fractured 
surfaces.  The  'purification*  of  the  surface  due  to  volatization  of  water, 
impurities  and  alkali,  and  the  more  favorable  molecular  orientation  of  the 
silica  structure,  are  possible  reasons  for  this.  The  presence  of  residual 
stress  at  the  melt  surface  due  to  non-uniform  cooling  or  composition  gradients 
can  also  influence  the  chemical  reactivity.  Otherwise,  very  little  is  known  in 
detail  about  these  phenomena. 

Silica  Glass 

The  direct  surface  analysis  of  silica  'melt'  surfaces  has  not  been 
prevalent  for  a  couple  of  reasons.  Primarily,  it  is  the  chemical  structure  of 
the  silica  surface  and  its  hydroxyl  content  which  is  of  most  interest.  With  the 
exception  of  infared  absorption  spectroscopy,  none  of  the  available  surface 
analysis  techniques  are  especially  sensitive  to  these  characteristics*.  In 
addition,  the  geometry  of  the  sample  becomes  important  for  many  surface 
analyses.  Silica  fibers  which  exhibit  very  high  strengths,  for  example,  would 
be  most  interesting  to  study.  However,  their  dimensions  3nd  geometry  are  not 
appropriate  for  most  surface  analytical  methods.  The  ATR  technique  would  be 
most  appropriate  providing  that  flat,  plane  parallel  sheets  with  pristine  melt 
surfaces  could  be  formed. 
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Another  factor  which  has  contributed  to  the  lack  of  fundamental  studies  of 
silica  melt  surfaces  is  the  difficulty  in  preparing  and  preserving  surfaces  for 
study.  Although  surface  reactions  of  silica  are  slow  in  the  ambient  atmosphere, 
they  do  occur.  Obviously,  they  must  be  prevented  in  order  to  positively 
identify  those  phenomena  which  occur  at  high  temperature. 

One  indication  of  the  chemistry  at  silica  'melt'  surfaces  is  provided  by 
glass  capillary  gas  chromatography  (GC)^.  This  method  is  sensitive  to  the 
adsorption  of  nanogram  quantities  of  gaseous  species  on  the  inside  'melt' 
surfaces  of  glass  capillaries.  Figure  7  compares  the  gas  chromatograms  for 
polar  test  mixtures  analyzed  in  fused  quartz  and  high  purity  silica^.  The 
tails  and  reduced  intensity  of  the  peaks  obtained  in  the  fused  quartz  capillary 
indicate  that  some  of  the  gaseous  species  are  being  adsorbed  whereas  the  high 
purity  silica  provides  an  almost  perfect  chromatogram.  The  point  is  that  even 
the  trace  levels  of  impurities  present  in  the  fused  quartz  have  influenced  the 
surface  chemistry.  Quite  likely,  some  of  the  impurities  and  bulk  hydroxyl 
groups  have  'segregated'  at  the  inside  surface  of  the  capillary  during  the  high 
temperature  manufacture.  In  a  sense,  the  use  of  (GC)  to  study  the 
chemistry  of  pristine  melt  surfaces  is  ideally  suited  for  comparison  with  the 
strength  behavior  of  fibers.  Here,  the  inside  surface  of  a  25  meter  capillary 
is  'sampled'  by  the  polar  test  mixture.  Surely,  the  same  high  energy  sites 
responsible  for  these  adsorption  phenomena  will  play  a  role  in  the  fatigue  of 
glass  fibers. 

It  should  also  be  emphasized  that  both  temporary  and  permanent  stresses  can 
exist  at  melt  surfaces  and  these  will  surely  influence  the  surface  chemistry. 

All  forming  methods  subject  the  molten  glass  to  tensile  stresses  which  are 

maximized  at  the  surface.  Even  after  cooling,  the  surface  will  be  under  stress 

1 3 

due  to  the  reduced  interatomic  distances  expected  there  . 
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Multicomponent  Silicate  Glasses 

Due  to  the  high  mobility  and  reactivity  of  modified  species,  compositional 
and  structural  differences  are  more  prevalent  at  the  melt  surfaces  of 
multicomponent  glasses.  A  considerable  number  of  surface  studies  have  been 
reported  for  the  melt  surfaces  of  fibers,  bottles,  flat  glass  etc*.  In 
general,  it  has  been  found  that  the  composition  profiles  and  hydration 
susceptibility  of  multicomponent  melt  surfaces  can  be  quite  variable. 

Figure  8  presents  SIMS  depth  profiles  for  the  surface  of  an  E-glass  ribbon. 
The  ribbons  were  prepared  with  a  special  bushing  simply  to  provide  flat  surfaces 
amenable  for  study.  These  ribbons  provide  a  convenient  'model'  for 
understanding  the  high  temperature  chemistry  associated  with  fiber  glass 
surfaces.  Although  the  quench  rate  may  be  somewhat  less  for  the  ribbons,  one 
notes  that  this  surface  -prepared  in  air-  is  modified  in  depth.  The  outermost  5 
nm  are  enriched  in  boron,  and  to  some  extent  sodium,  while  the  sub-surface 
region  is  enriched  in  calcium  and  silicon.  These  profiles  -expressed  as  cation 
to  oxygen  ratios-  were  compared  with  the  profiles  measured  for  a  'bulk'  fracture 
surface  of  the  same  E-glass.  All  of  the  ribbon  profiles  matched  the  bulk 
profiles  at  depths  greater  than  40  nm  with  the  exception  of  sodium.  The  Na/O 
ratio  was  2  to  3  times  lower  in  the  outermost  200  to  600  nm  of  the  ribbon  than 
in  the  bulk  E-glass  standard.  Thus,  the  ribbon  surface  exhibits  a  modified 
layer  at  the  outermost  20  nm  -where  sodium  and  boron  are  enriched-  and  a  sodium 
depleted  region  extending  to  600  nm.  Although  the  profiles  in  the  outermost  20 
nm  may  be  due  in  part  to  sputtering  effects  associated  with  the  analysis  , 
the  sodium  depletion  extending  to  600  nm  is  independent  of  sputtering  effects. 

In  contrast,  the  inside  surface  of  an  alumina-borosil icate  glass  tube 
exhibits  a  dealkalized  layer  which  extends  to  only  40  nm.  The  profiles  in 
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Figure  9  show  that  the  degree  of  dealkalization  is  quite  dramatic;  the  sodium 
concentration  was  reduced  by  two  orders  of  magnitude  here,  in  contrast  to  a  2  to 
3  times  reduction  for  the  ribbon  surface.  In  contrast  to  many  other  'commercial 
surfaces' ,  the  inside  surface  of  tubing  exhibits  the  most  extensive  surface 
depletion  of  alkali  (with  the  exception  of  intentionally  dealkalized  products). 
This  is  undoubtedly  related  to  the  fabrication  process  which  involves  a  high 
flow  of  'air'  through  the  hot  glass  tubing  and  thereby  enhances  the 
volatilization  of  alkali  species. 

Of  course,  one  can  intentionally  modify  the  composition  of  glass  surfaces 

through  surface  treatments.  It  has  been  demonstrated,  for  example,  that  the 

* 

inside  surface  of  soda-lime  capillary  can  be  improved  for  (GC)*  by  passing 
sulfur  dioxide  or  dif luoroethane  through  the  preform  and  capillary  during  the 
fabrication^.  The  use  0f  sulfur  dioxide^®  or  dif luoroethane^® 
is  well-known  for  treating  the  internal  surfaces  of  glass  containers.  All  of 
these  treatments  enhance  the  hydration  resistance  of  the  glass  surface.  In 
fact,  it  was  recently  shown  that  the  combined  use  of  sulfur  dioxide  and 
dif luoroethane  could  further  enhance  the  hydration  resistance  of  the  glass 
surface  .  It  is  quite  clear  that  these  modified  surfaces  are  dealkalized 
after  the  high  temperature  treatment.  However,  it  is  also  observed  that  the 
hydrogen  content  in  the  surface  is  reduced.  Therefore,  the  alkali  species  are 
not  simply  replaced  by  hydrogen  during  these  treatments.  It  is  primarily  the 
high  temperature  polymerization  and  relaxation  of  the  silicate  network  near  the 
surface  which  i3  responsible  for  the  improved  hydration  resistance  after  these 
treatments.  Thus,  it  is  important  to  recognize  that  not  only  the  composition, 


but  particularly  the  network  structure  of  the  glass,  may  be  modified  by  the 
treatment  of  surfaces  at  high  temperatures. 
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Obviously  Che  critical  question  is  whether  or  not  the  high  temperature 

modification  of  thin  surface  layers  can  influence  strength  behavior.  Budd,  et 
31 

al.  measured  the  stress  corrosion  susceptibility  for  the  internal 
surfaces  of  glass  bottles  with  and  without  internal  sulfur  treatment.  His 
results  are  duplicated  in  Table  1.  They  indicate  that  modification  of  even  a 

50nm  surface  layer  -at  least  for  a  melt  surface-  can  have  a  dramatic  influence 

32 

on  the  strength  and  dynamic  fatigue  behavior.  Yashchishin  and  Sharagov 
examined  the  impact  toughness  of  plate  glass  treated  at  elevated  temperatures 
with  sulfur  dioxide  and  also  dif luorodichloromethane.  They,  too,  observed  an 
increase  in  the  fracture  resistance  due  to  a  modification  of  thin  surface 
layers.  The  reasons  for  the  increased  strength  and  fatigue  resistance  are  by  no 
means  clear,  although  the  enhanced  hydration  resistance  provided  by  these 
treatment  points  to  the  importance  of  surface  chemistry.  Here,  it  can  only  be 
stated  that  the  surface  composition  and  structure  has  been  modified  by  a  high 
temperature  treatment,  and  at  the  same  time,  a  measurable  improvement  in  the 
strength  and  fatigue  behavior  are  observed. 

CLEAN  SURFACES 

It  is  quite  noticeable  that  little  effort  has  been  directed  at 
understanding  the  fundamental  nature  of  glass  surfaces  via  the  spectroscopic 
study  of  'ideal'  or  'clean'  surfaces.  Here,  the  objective  might  be  to 
understand  the  nature  of  the  dangling  bonds  at  a  fresh  fracture  surface  produced 
in  ultra  high  vacuum  and  the  relative  effects  of  surface  reconstruction, 
adsorption  of  bulk  impurities,  and  adsorption  of  gaseous  species  during  its 
passivation.  In  a  sense,  a  more  fundamental  understanding  of  ’ideal'  surfaces 
is  required  to  fully  define  the  reaction  mechanisms  responsible  for  the 
properties  exhibited  by  'real'  surfaces  during  processing  in,  or  exposure  to, 
particular  atmospheres. 
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Hochstrasser  and  Antonini  crushed  silica  glass  in  an  ultra  high 
vacuum  environment  (“10“*®  torr).  They  used  electron  spin  resonance  (ESR) 
to  characterize  the  paramagnetic  dangling  bonds  and  adsorption-desorption 
measurements  to  examine  interactions  with  CO^.  Silica  fracture  surfaces 
are  expected  to  contain  approximately  2-3  dangling  bonds  per  100 
The  ESR  analysis  detected  only  4 2  to  10Z  of  the  expected  number  of  =Si* 
dangling  bonds,  and  failed  to  detect  the  presence  of  any  =  SiO*  dangling 
bonds.  The  investigators  concluded  that  most  of  the  dangling  bonds  were 
eliminated  by  a  self-reconstruction  process  during  and/or  after  the  fracture 
event.  This  occurred  even  at  temperatures  as  low  as  -150°C.  It  is  most 
interesting  that  the  number  of  dangling  bonds  detected  by  ESR  increased  with  the 
purity  of  the  silica.  This  suggests  that  impurities  in  the  glass  are 
responsible  -at  least  in  part-  for  the  elimination  of  dangling  bonds.  At 
temperatures  of  the  order  300°C,  all  of  the  dangling  bonds  were  eliminated 
by  surface  reconstruction  and  the  out-diffusion  of  impurities. 

a  < 

More  recently,  Lacharrae  et  al.  used  Auger  electron  spectroscopy  to 
examine  the  fracture  surfaces  of  silica  glass  and  binary  sodium  and  potassium 
silicate  glasses.  They  found  a  local  enrichment  of  alkali  at  the  fracture 
surfaces  created  in  a  vacuum  of  10~*°  torr.  They  concluded  that  the 
fracture  of  the  glasses  had  followed  zones  of  local  alkali  enrichment  already 
present  in  the  bulk  glass.  This  interpretation  is  inconsistent  with  the  random 
network  theory  for  glass  structure,  but  interestingly,  it  would  be  consistent 
with  the  cluster  models  of  silicate  glass  structure  where  alkali  ions 
preferentially  "decorate"  the  interfaces  between  clusters.  Thus,  fracture 
surfaces  could  be  enriched  in  alkali  if  the  crack  follows  a  path  between  the 
clusters.  Nonetheless,  other  investigators  have  offered  a  different 
interpretation  for  the  local  enrichment  of  alkali  observed  at  clean  fracture 


surfaces. 
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Pantano,  eC  al.  have  used  ion  scattering  spectrometry  (ISS)  to 
examine  fracture  surfaces  of  alkali-silicate  glasses  created  in  ultra  high 
vacuum.  The  ISS  method  is  unique  because  it  is  sensitive  only  to  the  outermost 
monolayer  of  the  solid.  In  fact,  by  scattering  ions  off  the  surfaces  at  a 
glancing  angle  only  those  atoms  which  'protrude*  from  the  surface  structure  are 
detected.  Unfortunately,  conventional  ISS  is  insensitive  to  hydrogen. 
Nonetheless,  these  data  show  quite  clearly  that  alkali  ions  do,  in  fact, 
dominate  the  outermost  monolayer  of  these  clean  fracture  surfaces.  Under  the 
conditions  where  the  ISS  spectra  for  pure  silica  glass  exhibited  the 
stoichiometric  ratio  of  Si  and  0,  the  spectra  for  the  alkali  silicate  glasses 
exhibit  little  or  no  detectable  Si  and  0.  Obviously,  these  'ideal'  fracture 
surfaces  are  not  a  simple  projection  of  the  bulk  structure.  These  investigators 
have  proposed  chat  the  alkali  ions  assume  new  positions  in  the  outermost 
monolayer  wherein  they  shield  the  silicate  tetrahedra.  It  is  likely  that  this 
atomic  rearrangement  requires  only  relaxation  or  a  'jump'  of  alkali  species  into 
Che  surface  monolayer,  and  not  a  long  range  diffusion  process.  And  while  the 
presence  of  localized  charges  at  the  dangling  surface  bonds  would  enhance  this 
outward  relaxation  of  alkali  ions,  the  presence  of  charged  surface  sites  is  not 
a  necessity.  Any  highly  deformable  cation  will  have  a  tendency  to  relax  into 
the  surface  monolayer  because  it  can  reduce  the  force  field  perpendicular  to  the 

surface  and  provide  a  more  complete  coordination  structure. 

35 

Kelso,  et  al.  compared  the  experimental  ISS  spectra  for  clean 
fracture  surfaces  of  silica  and  potassium  silicate  glass  with  the  results  of  a 
molecular  dynamics  (MD)  simulation  of  a  surface  monolayer.  A  fracture  surface 
was  created  in  the  MD  simulation  by  removing  the  periodic  boundary  conditions  in 
one  dimension  and  then  allowing  the  glass  to  relax  at  300K.  The  MD  calculations 
suggested  that  oxygen  would  dominate  the  outermost  monolayer  of  silica  glass 


while  potassium  and  oxygen  would  dominate  the  outermost  monolayer  of  potassium 
trisilicate.  The  MD  results  were  in  reasonable  agreement  with  the  experimental 
ISS  data.  Interestingly,  the  MD  simulations  showed  a  distinct  depletion  of 
potassium  in  the  sub-surface.  The  profile  was  reminiscent  of  a  'double-layer'. 

The  most  significant  conclusion  to  be  drawn  from  these  studies  of  clean 
surfaces  is  that  in  the  absence  of  external  environmental  species,  the  glass 
surface  will  lower  its  free  energy  state  by  reconstruction  processes.  The 
studies  described  in  a  later  section  of  this  paper  suggest  that  this 
rearrangement  may  occur  during  a  fracture  event.  In  the  case  of  melt  surfaces, 
the  rearrangement  will  be  enhanced  due  to  the  elevated  temperatures,  for 
example,  during  fiber  drawing.  Obviously,  any  attempt  to  'passivate'  surfaces 
by  attaching  chemical  species  to  dangling  bonds  at  the  surface  must  be  carried 
out  simultaneously  with  the  creation  of  the  surface,  or  by  an  internal 
'Gibbsian'  adsorption  process. 

RELATIONSHIPS  BETWEEN  BULK  COMPOSITION  AND  STRUCTURE, 

SURFACE  CHEMISTRY  AND  FRACTURE  OF  GLASS 

It  should  be  quite  clear  at  this  point  that  the  surface  chemistry  is 
dependent  upon  both  the  environmental  conditions  and  bulk  composition  and 
structure  of  the  glass.  Ideally,  then,  one  would  like  to  examine  the  fracture 
of  glass  specimens  whose  surface  chemistry  -and  its  dependence  upon  bulk 
composition  and  structure-  could  be  controlled  and  characterized.  This  could 
provide  very  fundamental  information  about  the  fracture  behavior  of  glasses. 
Indeed,  there  have  been  many  measurements  of  strength  distributions  and  slow 
crack  growth  in  controlled  environments,  and  these  measurements  have  been  made 
on  glasses  of  various  'bulk'  compositions.  But  none  of  these  studies  have  been 
coupled  with  direct  analysis  of,  for  example,  fracture  surfaces  prepared  at 
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sub-critical  velocities  in  various  chemical  environments.  Moreover,  few  of 
these  studies  have  attempted  to  examine  the  strength  and  fracture  behavior  and 
its  dependence  on  surface  chemistry  in  glasses  whose  composition  and  structure 
are  systematically  modified.  Obviously,  these  are  not  easy  tasks,  but  may  be 

necessary  in  order  to  provide  insight  to  the  micromechanisms  of  crack  growth. 

36 

In  this  regard,  Kaplow  attempted  to  prepare  clean  silica  films  under 
ultra  high  vacuum  conditions,  analyze  the  surfaces  directly  with  Auger  electron 
spectroscopy,  and  then  measure  their  strengths  with  a  fracture-bend  device.  He 
proposed  further  to  modify  these  atomically  clean  silica  surfaces  by  exposure  to 
various  chemical  environments,  and  to  correlate  the  resulting  surface 
composition  directly  with  the  measured  strengths.  It  was  his  intention  to 
'passivate'  the  surface  by  reacting  the  dangling=Si+  andHSiO  surface 
sites  in  these  chemical  environments.  Unfortunately,  the  silica  films  he 
prepared  -although  atomically  clean-  contained  thermal  and  structural  flaws 
whose  magnitude  insensitized  the  strength  behavior  to  chemical  effects.  Of 
course,  the  studies  of  Antonni  and  Hochstrasser  ,  and  also  Kelso  , 
have  already  shown  that  these  dangling  bonds  are  rapidly  eliminated  even  under 
ultra-high  vacuum  by  a  self-reconstruction.  Thus,  it  is  unlikely  that  they 

would  have  been  available  for  passivation  under  these  conditions. 

•  24  37 

Martin,  et  al.  and  also  Donnett,  et  al  measured  the  strengths 

and  adsorption  isotherms  for  E-glass  and  soda-lime-silicate  glass  fibers  which 
were  drawn  under  varying  -but  controlled-  conditions.  In  principle,  at  least, 
this  affords  one  the  highest  probability  of  'passivating'  the  surface,  as 
well  as  understanding  the  way  high  temperature  surface  chemistry  influences 
strength  behavior.  Unfortunately,  neither  study  utilized  direct  surface 
analysis  to  characterize  the  composition  and  structure  of  the  fiber  surfaces. 

The  adsorption  isotherms  provide  only  an  indirect  measure  of  these  surface 
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properties,  and  moreover,  the  adsorbates  used  can  very  often  modify  the  surface 
under  study.  Nonetheless,  it  is  worth  emphasizing  that  both  groups  observed 
that  the  development  of  surface  microporosity  was  a  key  factor  in  the 

correlation  between  drawing  atmosphere  and  strength  behavior. 
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Chesson  and  Kerly  have  taken  a  different  approach  in  attempting  to 
correlate  bulk  composition  and  structure,  surface  chemistry  and  fracture 
behavior.  Here,  a  systematic  compositional  series  of  sodium  aluminosilicate 
glasses  have  been  studied  so  that  the  influence  of  various  chemical  and 
structural  features  could  be  evaluated  with  regard  to  fracture  and 
stress  corrosion  behavior.  In  the  system  Na^* xAl^*  (3-x)SiC>2,  the  bulk  density 
of  non-bridging  oxygen  sites  can  be  varied  independent  of  the  sodium  concen¬ 
tration.  Based  upon  the  current  understanding  of  hydration  -  already  summarized 
in  this  paper  -  these  structural  variations  should  cause  systematic  changes  in 
the  surface  pH,  type  of  surface  film  formation,  water  transport  mechanism  and 
dissolution  kinetics. 

Figures  10,  11  and  12  reveal  that  the  chemical  properties  of 
these  glasses  are,  in  fact,  very  sensitive  to  the  Al^O^/Na^O  ratio  (which 
controls  the  non-bridging  oxygen  density).  It  can  be  seen  that  these  glasses 
exhibit  a  wide  range  of  'equilibrium'  pH  values  (Figure  10).  Thus,  the 
crack-tip  chemistry  in  each  of  these  glasses  will  be  very  different  in  aqueous 
environments.  Although  one  might  argue  that  the  pH  at  the  tip  of  a  moving  crack 
may  not  reach  its  'equilibrium'  value,  the  equilibrium  pH  for  a  given  glass 
characterizes  its  'tendency'  to  react  with  water.  Indeed,  the  pH  of  the 
solution  contained  in  tht  crack  tip  may  be  of  little  importance;  rather,  the 
surface  pH  -  that  is,  the  number  of  ionized  silanol  groups  -  may  be  of  more 
significance.  Of  course,  the  implications  of  equilibrium  pH  with  regard  to  the 
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aging  of  surfaces  and  pre-exiscing  flaws  -  even  in  Che  absence  scress  -  are  more 
straightforward. 

The  soda  and  silica  dissolution  rates  are  also  extremely  sensitive  to  the 
Na^O/AljO^  ratio  as  shown  in  Figure  11.  This  is  an  intrinsic  property  and  not 
a  pH  effect.  That  is,  the  dissolution  rates  were  measured  in  water  -  where  the 
pH  could  drift  -  and  at  constant  pH  using  a  buffer.  The  dissolution  rates  were 
slightly  higher  under  constant  pH  conditions,  but  the  dependence 

upon  the  Al20j/Na20  ratio  was  not  affected.  This  strong  dependence  of  dissolution 
rate  upon  the  Al20^/Na20  ratio  -  independent  of  'solution  pH'  -  lends  validity 
to  the  suggested  importance  of  'surface  pH'. 

The  hydrogen  profiles  in  Figure  12  provide  some  insight  to  surface  pH 
effects,  the  transport  of  water  and  surface  film  formation.  They  show  that  the 
glasses  with  x  =  0  and  x  =  .2  form  distinct  hydrated  layers  at  the  surface, 
whereas  the  glasses  with  x  ^  .4  exhibit  reduced  hydrogen  penetration.  More¬ 
over,  the  amount  of  hydrogen  incorporated  in  the  leached  layer  decreases  with 
the  Al2C>2/Na20  ratio.  These  same  glasses  were  examined  with  diffuse  reflectance 
Fourier  transform  infared  spectroscopy.  These  data  revealed  that  the  sodium 
ions  at  non-bridging  oxygen  sites  were  released  at  a  much  higher  rate  than 
those  at  (AlO^)  sites.  Thus,  although  all  of  these  glasses  have  the  same  num¬ 
ber  of  sodium  ions  -  and  therefore  the  same  number  of  sites  where  protonic 
species  may  be  immobilized  -  it  is  primarily  the  non-bridging  oxygen  sites  which 
influence  the  hydration  phenomena  and  surface  pH.  That  is,  it  is  the 
non-bridging  oxygen  site,  and  not  the  sodium  ion  per  se ,  which  is  responsible 
for  the  reaction  with  water. 

The  large  effect  of  the  first  alumina  substitution  (x  =  0.2)  on 
the  chemical  properties  suggests  that  the  Al^  does  not  substitute  randomly 
in  the  Na^o  3Si02  glass  structure.  Rather,  it  occupies  those  sites  where  the 


transport  and  attack  of  water  is  most  significant,  for  example  within  a  topo¬ 
logical  network  of  non-bridging  oxygens.  This  conclusion  is  consistent  with 
the  results  of  Smets  and  Lommen^  who  measured  the  diffusion  coefficient  for  water  (D^ 

during  hydration  of  a  similar  compositional  series  of  aluminosilicate  glasses. 

They  report  that  a  change  in  the  Al203/Na20  ratio  from  0  to  .25  in  Na20.4Si02 
reduced  D„  n  by  a  factor  of  400. 

A  wide  variety  of  mechanical  properties  are  also  being  evaluated  in  this 
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same  series  of  glasses  .  Here,  only  the  fracture  toughness  (KI(,)  and  stress 
stress  intensity  exponent  (n)  are  reported  in  Figures  13  and  14.  The  fracture 
toughness  -  measured  in  liquid  nitrogen  via  the  controlled  microflaw 
technique  -  increased  with  the  first  substitution  of  AljOj  in  Na20*3Si02  and  then 
decreased  with  additional  substitution.  The  initial  increase  in  K  is  consistent 

X.V/ 

with  the  idea  that  the  first  additions  of  A1203  to  Na20'3Si02  do  not  substitute 
randomly  in  the  structure,  but  rather  occupy  the  weakest  sites  -  topological  or 
otherwise  -  in  the  Na  0*3Si0  structure.  The  decrease  in  K  with  further 
polymerization  of  the  network  structure  is  more  surprising  -  particularly  since 
the  elastic  modulus  increases  steadily  with  Al^t^/NajO.  However,  it  is  notable 
that  the  activation  energy  for  sodium  ion  motion  (@T<Tg)  exhibits  the  same 
dependence  upon  the  Al203/Na20  ratio.  The  activation  energies  were  obtained 
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independently  for  this  same  series  of  glasses  via  internal  friction  measurements 
This  correlation  suggests  that  the  motion  of  alkali  ions  may  be  involved  in  the 
propagation  of  cracks  -  at  least  in  the  absence  of  external  chemical  effects. 

This  is  consistent  with  the  observations  made  on  the  alkali-silicate  fracture 
surfaces  created  in  vacuum  where  an  alkali-ion  rearrangement  was  observed.  It 
leads  naturally  to  the  conclusion  that  the  atomic  rearrangement  observed  at 
those  clean  fracture  surfaces  occurred  during  the  crack  propagation  in  ultra 
high  vacuum.  The  fact  that  the  glasses  with  Al,0,/\ao0  =  1  have  fully  polymerized 
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networks  such  as  pure  silica  glass  -  but  exhibit  lower  fracture  toughnesses  than 
pure  silica  glass  -  is  also  consistent  with  the  idea  that  sodium  ions  can  be 
involved  in  the  fracture  process.  Interestingly,  the  notion  that  alkali  ions 
play  a  role  in  fracture  was  first  suggested  by  Cox  in  1969^.  Later, 

Rauschenbach  and  Hinz^  considered  the  stress-induced  segregation  of  sodium 
to  crack-tips  in  glass. 

The  stress  intensity  exponents  in  Figure  IV-14  were  measured  in  deionized 
water  using  the  stressing  rate  technique.  These  data  suggest  that  it  is  the 
combined  effect  of  crack  tip  pH  and  network  structure  -  and  not  the  mere 
presence  of  alkali  -  which  controls  the  susceptibility  to  stress  corrosion. 

That  is,  the  Na^O* A^Oj*  2SiC>2  glass  exhibits  a  susceptibility  to  stress  corrosion  which 
is  comparable  to  SiC^.  It  should  be  recalled  that  the  glasses  have  isostructural 
networks  and  that  the  equilibrium  surface  pH  of  the  aluminosilicate  is  extreme¬ 
ly  low  in  spite  of  the  soda.  Of  course,  additional  studies,  particularly  slow 
crack  growth  experiments  carried  out  in  solutions  buffered  to  the  'equilibrium 
pH' ,  are  needed  to  fully  rationalize  this  behavior. 


SUMMARY 

It  is  quite  clear  that  chemical  reactions  at  glass  surfaces  -  under  ambient 
conditions  as  well  as  during  high  temperature  fabrication  processes  -  will 
modify  the  composition,  molecular  structure  and  morphology  of  the  surface  to 
varying  degrees.  It  is  also  evident  that  the  bulk  glass  composition  and 
structure  will  influence  the  reaction  mechanisms,  and  in  particular,  the  'depth' 
of  interaction.  However,  the  extent  to  which  these  surface  chemical  phenomena 
directly  affect  strength  and  fracture  behavior  has  not  been  established. 

Indeed,  correlations  may  be  drawn  between  the  surface  chemical  reactivity  and 
hydration  susceptibility  of  various  glasses  and  glass  surfaces  and  their 


strength  and  fracture  behavior.  Nonetheless,  the  fundamental  processes  which 
underly  these  apparent  correlations  are  not  accessible  simply  from  a  surface 
analysis.  Even  the  reconstruction  processes  observed  to  occur  at  fracture 
surfaces  of  alkali-containing  glasses  may  be  a  consequence  -  and  not  the  means  - 
by  which  fracture  occurs. 

It  is  likely  that  surface  chemical  reactions  which  occur  under  stress  are 
the  critical  phenomena  influencing  the  strength  and  fracture  behavior  of  glass. 
Of  course,  these  stresses  need  not  be  limited  to  an  applied  stress  but  un¬ 
doubtedly  include  the  intrinsic  'monolayer'  stresses  associated  with  the  very 
existence  of  a  surface,  surface  stresses  due  to  thermal  history  and  dealka¬ 
lization  of  multicomponent  glass  surfaces,  and  the  stressed  bonds  associated 
with  the  topology  and  impurities  in  the  surface  region.  The  important  point  to 
be  made  is  that  the  classical  and  modern-day  spectroscopic  approaches  to  surface 
science  are  insensitive  to  the  presence  of  the  stresses.  That  is,  surface 
analyses  cannot  be  localized  at  microscopic  regions  of  high  stress.  And, 
uniform  surface  stresses  whose  magnitude  is  of  significance  in  fracture  cannot 
be  achieved  in  situ  on  specimens  which  are  convenient  for  surface  studies  (with 
the  possible  exception  of  fibers  and  ribbons).  Even  the  direct  analysis  of 
surfaces  created  in  a  slow  crack  growth  experiment  is  of  questionable  value 
because  the  surface  species  produced  at  fracture  may  be  intermediate  to  the 
surface  species  found  by  an  analysis  after  fracture. 

These  comments  are  not  meant  to  discourage  surface  chemical  studies  for 
understanding  strength  and  fracture  behavior.  To  be  sure,  knowledge  of  the 
composition  of  glass  surfaces  and  their  sensitivity  to  water,  pH  effects  and  the 
bulk  glass  properties  is  necessary.  And  in  a  practical  sense,  these  approaches 
are  required  for  the  characterization  of  surfaces,  surface  treatments  and 
protective  coatings.  However,  the  direct  re lat ionships  which  can  be  drawn 
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between  classical  surface  chemistry  and  surface  analyses,  and  the  fundamental 
micromechanisms  of  fracture,  may  be  limited.  Nevertheless,  it  is  hoped  that 
this  review  provides  a  useful  summary  of  glass  surface  chemistry,  as  it  might 
relate  to  the  strength  and  fracture  of  silicate  glasses. 
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Figure  1 


Figure  2 


Infrared  spectra  for  porous  glass  showing  absorption  bands  for 
molecular  water,  vicinal  hydroxyl  groups  and  freely-vibrating 
hydroxyl  groups;  (A)  heated  to  100°C  (B)  heated  to  300°C  (C)  heated  to 
500°C  (D)  re-hydrated  after  heating  to  500°C  (from  7). 


b  .  Attenuated  Total  Reflectance  (ATR) 


.  Methods  of  infrared  absorption  analysis  which  are  'surface  sensitive' 
(A)  transmission  through  a  microporous  compact  where  many 
surfaces  are  sampled 

(3)  attenuated  total  reflectrance  (ATR)  where  the  surface  is 
sampled  by  multiple  internal  reflections. 
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Figure  3 
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Schematic  representations  of  silicate  glass  surface  structures 
containing  silanol  groups  and  adsorbed  molecular  water  (al  silica 
surface  under  ambient  conditions  (b)  alkali-silicate  surface  under 
ambient  conditions  (c)  alkali-silicate  surface  hydrated  in  water. 
(The  ideal  random  networks  in  these  schematics  were  generated  with 
triangle  rafts . ) 
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Figure  4.  Hydrogen  and  alkali  depth-profiles  obtained  with  sputter-induced 
photon  spectroscopy  (SIPS)  after  soaking  commercial  glasses  in 
water  for  48h  at  90°C  (a)  Corning  code  7190  (trace  alkali) 

(b)  Coming  code  7740  (4%  soda)  (c)  sheet  glass  (12.8%  soda) 

(d)  Corning  code  0160  (11.3%  pottasia,  2.6%  soda)  (e)  Corning 

code  0080  (16.2%  soda)  (f)  K-profile  for  Corning  code  0160.  (from  19) 


Figure  5.  Model  for  the  reaction  between  a  clean  sodium-silicate  glass 
surface  and  water  vapor  in  the  ambient  atmosphere  (A)  sodium 
depth-profile  for  the  clean  surface  (B)  idealized  sodium  depth- 
profile  after  reaction  with  adsorbed  water  and  carbonaceous 
species,  and  the  corresponding  hydrogen  incorporation. 
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Figure  6. 


Surface  layer  formation  on  sodium-silicate  glass.  Hydrogen  glass 
corresponds  to  the  initial  replacement  of  the  alkali  species  with 
hydrogen  to  produce  silanol  groups  at,  and  within,  the  glass  surface. 
The  gel  is  due  to  the  incorporation  of  molecular  water  (at  the 
silanol  groups  in  the  hydrogen  glass)  and  the  subsequent  trans¬ 
formation  of  the  surface  into  a  swollen  layer  of  colloidal  silica. 
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Figure  7.  Glass  capillary  gas  chromatographic  (GC)“  analysis  of  a  nonpolar 
test  mixture  in  fused  quartz  and  nigh-purity  silica;  tails  on  the 
peaks  and  their  reduced  intensity  in  the  fused  quartz  indicates 
'adsorption’  of  these  species  (from  Zb). 
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Figure  8.  Depth  profiles  obtained  with  secondary  ion  mass  spectroscopy 
(SIMS)  at  the  surface  of  an  E-glass  ribbon.  (The  Na/O 
intensity  ratio  measured  on  a  bulk  standard  of  E-glass  is  -'15 
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Figure  9.  Depth  profiles  obtained  with  secondary  ion  mass  spectroscopy 
(SIMS)  at  the  inside  surface  of  commercial  glass  tubing. 
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ABSTRACT 

Crack  velocity,  V,  was  measured  as  a  function  of  stress  intensity,  K  , 
for  a  series  of  sodium-aluminosilicate  glasses  corresponding  to  compositions 
represented  by  Na20 'XAl^O^ ‘O-xJSiO^  with  x  values  of  0.0,  0.2,  0.4,  0.6, 

0.8  and  1.0.  Thus,  the  molar  concentrat ion  of  Na20  was  held  constant 
while  A^O^  was  substituted  for  an  equimolar  amount  of  SiO  .  Although  the 
corrosion  resistance  of  these  glasses  varies  by  at  least  three  orders  of 
magnitude,  the  slopes  (n)  of  the  log  V-log  curves  only  changed  from  ~17  to 
"25.  These  observations  indicate  that  ion  exchange,  leaching  and  dissolution 
are  not  the  primary  factors  which  control  slow  crack  growth. 

Although  the  underlying  mechanism  of  slow  crack  growth  is  a  chemical 
reaction  between  environmental  species  and  crack  tip  bonds,  it  is  suggested  that 
the  rate  of  the  chemical  reaction  is  controlled  through  the  elastic  strain. 

Under  the  assumption  that  the  slope  (n)  of  the  log  V-log  K^  curve  varies  as  the 
stiffness  in  the  crack  tip  region,  an  expression  was  derived  for  the  stiffness 
in  terms  of  Young's  elastic  modulus  (E).  An  excellent  correlation  was  observed 
between  n  and  E  for  sodium  aluminosilicate,  soda-lime-silica  and  fused  silica 
glasses  in  agreement  with  our  hypotheses.  Altogether,  this  implies  that  the 
expression  for  the  relation  between  crack  velocity  and  stress  intensity  could  be 
written  as: 


V  =  V  k?se 
o  I 

where  V  and  n are  constants  and  S  and  E  are  material  parameters  determined  by 
the  elastic  constants.  Thus,  the  composition  dependence  of  the  stress  corrosion 
susceptibility  is  due  to  the  composition  dependence  of  the  elastic  constants  and 
not  the  chemical  corrosion  properties. 
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INTRODUCTION 

The  dependence  of  mechanical  strength  upon  time  and  environment  has  been 
known  for  some  time  (1).  The  mechanism  of  strength  degradation  in  glasses  has 
been  reported  to  be  crack  extension  at  sub-critical  velocities  (2).  The  rate  of 
crack  extension,  and  hence  the  magnitude  of  strength  degradation,  is  governed 
primarily  by  the  applied  stress  and  stress-corrosion  susceptibility.  The 
objective  of  this  study  was  to  establish  the  dependence  of  this  slow-crack 
growth  behavior  upon  the  composition  and  network  structure  of  the  glass. 

BACKGROUND 

It  is  generally  accepted  that  crack  extension  occurs  via  a  chemical 
reaction  between  the  material  and  the  environment  surrounding  it.  Charles  and 
Hillig  (3)  proposed  that  the  rate  at  which  such  a  reaction  proceeds  is  altered 
by  the  application  of  an  external  stress.  Wiederhorn  (2)  suggested  that  the 
external  stress  offsets  a  portion  of  the  required  activation  energy  in  such  a 
manner  to  increase  the  reaction  rate.  Since  the  applied  stress  is  intensified 
at  the  crack  tip,  the  reaction  rate  at  the  tip  would  be  expected  to  be  large 
compared  to  the  surrounding  planar  surfaces.  Under  these  conditions,  Wiederhorn 
proposed  the  following  model: 

V  =  V  exp  [(-E  +  bK_ ) /RT ]  (1) 

O  O  X 

where  V  =  rate  of  crack  extension 

V  =  constant 
o 

E  =  activation  energy  in  the  stress  free  state 
o 

b  =  constant 

Kj  =  stress  intensity 

while  R  and  T  are  the  gas  constant  and  absolute  temperature. 


This  model  has  been  reported  to  describe  experimental  stress  intensity  - 

crack  velocity  data  very  well  (4,5)  and,  consequently,  has  gained  widespread 

acceptance.  Nevertheless,  most  researchers  characterize  the  slow  crack  growth 

behavior  by  the  n-parameter  which  relates  the  stress  intensity  to  the  crack 

* 

velocity  through  the  following  relationship: 

V  =  V  (K  )n  (2) 

O  I 

The  n  parameters  reported  for  silicate  based  glasses  range  from  10  to  40  (6). 

A  number  of  attempts  have  been  made  to  determine  the  feature  or  property  of 
a  glass  which  is  fundamental  to  the  observed  n-parameter.  The  tendency  of  a 
glass  to  react  with  an  aqueous  environment  or  undergo  slow  crack  growth  has  been 
examined  through  the  aqueous  corrosion  behavior  (7,9).  One  of  these  studies 
showed  that  the  leaching  rate  of  sodium  ions  varied  inversely  with  the 
n-parameter  (7).  The  linear  thermal  expansion  coefficient,  which  is  related  to 
the  molecular  structure  of  the  material,  was  also  shown  to  vary  inversely  with 
the  n-parameter  for  many  oxide  glasses.  These  observation  suggest  that  the  two 
properties  which  play  a  major  role  in  determining  the  slow  crack  growth  behavior 
are  chemical  corrosion  and  molecular  structure.  In  the  present  study,  the 
effects  of  these  two  characteristics  were  separated  in  an  effort  to  determine 
which  is  dominant  in  the  slow  crack  growth  process.  Sodium- 

•ff 

It  can  be  argued  that  Eq  (1)  has  more  fundamental  meaning  than 
Eq  (2)  because  it  relates  crack  velocity  to  the  activation  energy  and 
activation  volume;  nonetheless,  it  is  more  common  and  convenient  to 
use  n,  the  slope  of  the  log  velocity-log  stress  intensity  curve,  as  a 
measure  of  the  resistance  to  slow  crack  growth. 
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aluminosilicate  glasses  were  chosen  for  investigation  because  the  structure  and 
corrosion  behavior  can  be  systematically  varied  by  changing  the  soda  to  alumina 
composition  ratio. 


EXPERIMENTAL  PROCEDURE 

Sodium-aluminosilicate  glasses  corresponding  to  compositions  represented  by 
the  general  formula: 


Na2o*x  A1203*(3-x)  Si02  (3) 

were  prepared  with  x  values  of  0.0,  0.2,  0.4,  0.6,  0.8  and  1.0.  In  this  series, 
the  molar  concentration  of  Na2o  was  held  constant  while  Al203  was  substituted 
for  an  equimolar  amount  of  Si02<  These  glasses  are  identified  by  their 
Al203/Na20  molar  ratios;  thus,  Na20*3Si02  is  labeled  0.0  glass  and  Na20  • 
Al203*2Si02  is  labeled  1.0  glass. 

All  of  the  glasses  were  melted  in  platinum  or  irridium  crucibles  using 
silica  (5  Pm  Minusil)  and  reagent  grade  aluminum  hydroxide  and  sodium  carbonate. 
The  high  alumina  glasses  were  crushed  and  remelted  two  or  three  times  to  insure 
homogeneity.  After  annealing,  the  glasses  were  cut  into  plates  approximately 
7.5  cm  x  1.27  cm  x  0.25  cm;  a  groove  was  machined  down  the  center  of  the 
specimen.  The  depth  of  the  groove  was  approximately  one-half  of  the  plate 
thickness  and  the  width  was  0.00075  to  0.001  cm.  After  cutting  and  grooving, 
the  specimens  were  re-annealed  to  remove  any  residual  stresses,  and  then  were 
pre-stressed  to  initiate  a  sharp  crack. 

The  stress  intensity-crack  velocity  relationship  for  each  composition  was 
obtained  by  the  constant  moment  double  cantilever  beam  technique  (CM-DCB)  as 
described  by  Freiman,  Mullville,  and  Mast  (10).  The  aluminum  arms,  through 
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which  the  moment  was  applied,  were  attached  to  the  glass  specimens  with  epoxy. 
The  applied  stress  intensity,  K^,  was  calculated  from  the  relationship  derived 
by  Freiman,  Mullville,  and  Mast: 

Kx  =  [PL]/[2  VTt  3  (4) 

where  P  =  applied  load 

L  =  length  of  moment  arm 
I  =  moment  of  inertia 

t  =*  specimen  thickness  in  the  groove  (1.25  mm). 

The  crack  extension  was  monitored  with  a  40X  travelling  microscope  equipped  with 
a  filar  eyepiece. 

The  crack  growth  behavior  of  each  composition  was  examined  in  both 

_3 

deionized  water  (pH  5)  and  in  10  M  NaOH  solution  (pH  11).  At  least  two,  and 
up  to  five,  specimens  of  each  glass  were  tested  at  each  condition.  Each 
specimen  generally  provided  10-40  individual  data  points.  The  solution  pH  was 
checked  periodically  throughout  each  test. 


RESULTS 


Crack  Velocity  Data  in  Water 


The  results  of  constant  moment  double  cantilever  beam  experiments  in 
deionized  water  (pH  5)  on  0.0,  0.2,  0.4,  0.6,  0.8,  and  1.0  glasses  are  given  in 


Figures  la,  lb,  lc ,  Id,  le,  and  If,  respectively.  The  results  have  been  plotted 
according  to  the  classical  reaction  model  (Equation  1)  of  Wiederhorn  (2).  The 


corresponding  regression  lines  are  depicted  in  each  figure,  as  well  as  the 
values  of  the  slope  (b/RT)  and  the  95 %  confidence  intervals  calculated  from  the 
linearized  version  of  Equation  1: 


ITY  (m/s) 


Ln  CRACK  VELOCITY  (m/s) 


Nq20  0.4Al203  -2.6Si02 
water  (pH  =  5) 


SLOPE  =  255*  3.0 
N  =  22.5*  2.3 


K  (MPa  •  mV2) 


Figure  1c.  Stress  intensity-crack  velocity  data  for  0.2  glass 
in  water. 
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Log  V  =  [(Log  V  )-E*/RT]  +  [b/RT]K  (5) 

o  I  1 

The  n-parameters  of  the  linearized  power  law  relationship: 

Log  V  =  Log  V  +  f n I  Log  K.  (6) 

o  l 

are  also  given  in  the  plots.  In  each  case,  the  data  were  fit  to  a  two  parameter 
regression  model  (11).  For  convenience,  the  slopes  and  n-values  have  been 
listed  in  Table  I. 

All  compositions  exhibited  a  linear  trend  between  the  log  of  the  crack 

1/2 

velocity  and  the  applied  stress  intensity  in  the  range  K  =  .25-. 55  MPa*m  ; 

-9  -4 

the  corresponding  crack  velocities  were  in  the  range  10  to  10  m/s.  The 

1/2 

0.0  glass  exhibited  a  plateau  in  the  range  =  0.2-0.25  MPa>m  at  a 
-9 

velocity  of  10  m/s.  This  type  of  behavior  has  been  reported  to  occur  in  this 
composition  by  others  (9).  All  of  the  glasses  exhibited  thresholds  at  low 
(commonly  termed  the  stress  corrosion  limit  or  static  fatigue  limit).  An 
interpretation  of  these  thresholds  on  the  basis  of  crack  healing  has  been 
reported  (12). 

-3 

Crack  Velocity  Data  m  10  M  NaOH 

Figures  2a,  2b,  2c,  2d,  2e  and  2f  show  the  stress  intensity-crack  velocity 

-3 

data  sets  obtained  in  a  10  M  NaOH  aqueous  solution  (pH  11)  for  0.0,  0.2,  0.4, 
0.6,  0.8,  and  1.0  glasses,  respectively.  The  slopes  and  n-parameters  for  each 
composition  are  also  given. 

None  of  the  glasses  showed  any  significant  change  in  slope  or  n-parameter 

at  this  higher  pH.  However,  the  velocity  plateau  in  the  0.0  glass  occured  at  a 

“8  1  / 
higher  velocity  (10  m/ s )  and  extended  up  to  a  stress  intensity  of  0.38  MPa*ml/ 


CRACK  VELOCITY  (m/s) 


\ 


0.15 

0  25  0.35 

0.45 

0  55  0.6 

K. 

(MPa-m1/2) 

Figure  2a. 

Stress  intensity-crack  velocity  data 
in  NaOH  solution. 

for  0.0  glass 
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Figure  2c.  Stress  intensity-crack  velocity  data  for  0.4  glass 
in  NaOH  solution. 
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Figure  2d.  Stress  intensity-crack  velocity  data  for  0.6  glass 
in  NaOH  solution. 


DISCUSSION 


Figure  3  plots  n-parameters  of  the  cumulative  data  sets  versus  the 
composition  of  the  specimens  expressed  as  the  Al/Na  ratios.  The  n-parameters 
range  from  16  to  23  which  is  typical  of  silica-based  glasses  (6).  The 
n-parameter  increased  with  A^O^  concentration  up  to  an  Al/Na  molar  ratio  of 
0.6.  The  n-parameter  was  invariant  for  compositions  with  Al/Na  molar  ratios  in 
excess  of  0.6. 

The  two  characteristics  which  are  expected  to  play  a  major  role  in 
determining  the  slow  crack  growth  behavior  are  the  glass  structure  and  the 
chemical  corrosion.  First,  the  corrosion  behavior  of  these  glasses  will  be 
addressed,  and  then  the  structural  characteristics  will  be  considered.  Finally, 
a  relationship  will  be  established  between  the  elastic  properties  and  the  slow 
crack  growth  parameter.  It  will  be  shown  that  the  glass  structure  is  the 
determining  factor  in  the  slow  crack  growth  behavior  whereas  corrosion  effects 
are  important  only  at  very  low  K^. 

Corrosion  Effects 

It  is  now  well  accepted  that  the  attack  of  alkali-silicate  glasses  in 
aqueous  solution  is  due  to  two  pH-dependent  mechanisms  (13,14).  In  solutions  of 
pH  less  than  9.8,  ion  exchange  and  hydration  of  the  glass  surface  are  the 
dominant  effects  whereas  in  solutions  of  pH  >7,  dissolution  of  the  glass  occurs. 
The  corrosion  behavior  of  alkali-silicate  glass  can  be  explained  over  a  wide 
range  of  conditions  by  a  combination  of  these  two  mechanisms. 

In  acidic  and  neutral  pH  solutions,  the  primary  corrosion  mechanism  is  one 
of  ion  exchange  between  sodium  ions  in  the  glass  and  hydrogen  or  hydronium  ions 
in  the  solution.  Such  an  exchange  can  be  represented  by  the  following  reaction: 
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Figure  3 
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n-parameter  as  a  function  of  composition  in  the  glass 
system  Na.O. xAl  0  • (3-x)SiO  . 
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I  SiONa  .  +H  «(H-0)  .  „ .  -*■  =  SiOH  . 

glass  2  solution  glass 


+  Na  *(H_0)  .  ..  „ 

2  solution 


(7) 


Although  the  exchange  process  does  not  directly  affect  the  silicate 
network,  the  concentration  of  hydrogen  ions  in  solution  is  lowered  and  the  glass 
surface  is  depleted  of  sodium  ions.  Under  dynamic  conditions,  the  contacting 
solution  can  be  replenished  such  that  a  'pure  water1  solution  is  present  at  the 
surface  at  all  times.  In  this  case,  the  exchange  would  occur  continually  and  be 
limited  only  by  the  interdiffusion  of  hydrogen  species  and  sodium  ions  in  the 
glass  surface.  Conversely,  under  static  conditions  the  attacking  solution  is 
not  replenished  and  the  pH  of  the  contacting  solution  increases  to  an 
equilibrium  value  which  is  in  excess  of  the  original  pH.  The  magnitude  of  this 
equilibrium  pH  is  determined  by  the  composition  of  the  glass  and  the  glass 
surface  area  to  solution  volume  ratio. 

In  high  pH  solutions,  and  to  some  extent  in  all  aqueous  solutions,  hydroxyl 
ions  attack  siloxane  bonds  and  dissolution  of  the  silicate  network  occurs;  e.g. 


B 

s 

i 

81 


(-Si  "  0  ~  SiE)glass  +  (0H  5 


solut ion 


(=Si-0~ ) 


glass 


+  (HO  -  Si=) 


solution 


(8) 


The  silicate  solution  species  whose  solubility  controls  this  dissolution 
reaction  is  monos ilicic  acid  (H^SiO^).  Thus,  the  pH  dependence  of  the 
dissolution  reaction  can  be  represented  in  terms  of  the  equilibrium  constant  for 
its  ionization;  i.e. 


H.SiO,  +  H,SiO  "  +  H+ 


(9) 


V> 
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(10) 


(11) 


The  concentration  of  monosilicic  acid  in  solution  decreases  due  to  ionization, 
and  in  order  to  maintain  equilibrium,  more  silicate  network  is  dissolved.  Since 
the  ionization  becomes  substantial  at  pH  >9.8  (see  equation  11),  so  too  does 
dissolution  of  the  silicate  network. 

Qualitatively,  the  corrosion  mechanism  in  the  sodium-aluminosilicate 
glasses  will  follow  the  general  description  presented  above  for  the  simple 
alkali  silicate  glasses.  However,  the  presence  of  the  aluminate  species 
influences  the  pH  dependence  of  the  dissolution;  i.e.  at  low  pH,  the  soluble 
aluminum  ion  (Al^+)  can  form.  This  enhances  the  solubility  and  dissolution  of 
aluminosilicate  glasses  in  acidic  solutions.  Otherwise,  the  behavior  in  water 
and  basic  solutions  are  comparable. 

Figures  4  and  5  present  the  equilibrum  pH  values  and  dissolution  rate 
constants  for  the  sodium-aluminosilicate  glasses  in  water.  It  is  quite  clear 
that  the  reactivity  of  these  glasses  -  in  water  -  decreases  with  increasing 
Al/Na  ratio.  The  equilibrium  hydrogen  ion  activity  varies  by  4-5  orders  of 
magnitude  through  the  compositional  series,  while  the  sodium  leach  rate  and 
silica  dissolution  rate  varies  by  3  orders  of  magnitude.  All  these 
aluminosilicate  glasses  had  an  identical  concentration  of  soda;  namely,  25  mole 
percent.  Nevertheless,  they  exhibited  wide  variations  in  corrosion  resistance 


Figure  A.  The  equilibrium  of  pH  of  Na  O-xAl  0  •  (3-x)SiO  glasses 

in  contact  with  water  at  30  C;  bulk  glass  specimens  where 
the  surface  area  of  glass-to-solution  volume  ratio  was 
10" 1  cm-1. 


Moles  Na20  per  m  /Hr 


and  crack-growth  susceptibility.  Thus,  concepts  and  theories  which  might 
attribute  these  properties  to  the  concentration  of  sodium  ions,  per  se,  cannot 
be  altogether  correct.  In  the  context  of  this  study,  though,  the  question  is: 
can  this  wide  variation  in  chemical  corrosion  properties  be  used  to  rationalize 
the  crack-growth  behavior? 

At  first  glance,  one  notes  the  inverse  correlation  between  the  n-parameters 
(Fig.  3)  and  the  chemical  reactivity  (Figures  4  and  5);  i.e.  the  n-parameter 
increases  -  at  least  qualitatively  -  as  the  extent  of  pH  change  and  glass 
dissolution  decreases.  A  similar  correlation  was  reported  by  Tummala  (7),  and 
Wiederhorn  (15)  attributed  this  to  changes  in  crack-tip  pH.  That  is,  glasses 
with  higher  alkali  leach  rates  would  lead  to  a  significant  difference  in  pH 
between  the  localized  crack-tip  region  and  the  bulk  solution.  This  is  a 
plausable  suggestion  which  is  clearly  consistent  with  the  data  reported  here. 

The  0.0  and  0.2  glasses  exhibit  higher  equilibrium  pH  values,  and  thus,  one 
would  expect  the  solution  near  the  crack-tip  in  these  glasses  to  be  alkaline 
relative  to  the  bulk  solution.  And  correspondingly,  these  two  glasses  exhibit 
the  lowest  n-values.  Nevertheless,  this  apparent  correlation  is  not  consistent 
with  the  behavior  of  high  Al/Na  ratio  glasses.  The  n-parameter  for  these 
glasses  were  found  to  be  the  same  at  pH  5  and  pH  11.  Yet,  these  glasses 
dissolve  readily  at  pH  11.  Similarly,  the  0.0  and  0.2  glasses  -  whose  crack  tip 
pH  values  are  expected  to  be  controlled  by  the  high  leach  rate  rather  than  the 
bulk  solution  -  showed  differences  in  the  low  region  between  pH  5  and  pH  11. 
Altogether,  these  observations  indicate  that  ion-exchange,  leaching  and 
dissolution  are  not  the  primary  factors  which  control  slow  crack  growth.  These 
observations  were  made  in  a  family  of  glasses  whose  compositions  were 
systematicaly  varied,  and  thus,  the  lack  of  complete  consistency  with  the 
chemical  corrosion  model  lends  further  validity  to  this  conclusion. 


Structure  Effects 


It  is  widely  accepted  that  as  AljO^  is  added  to  soda-silica  glass,  the 
non-bridging  oxygen  concentration  decreases.  There  is  some  controversy 
concerning  the  range  of  Al/Na  ratios  over  which  this  structural  modification 
occurs,  but  one  of  these  studies  revealed  a  complete  elimination  of  the 
non-bridging  oxygens  at  an  Al/Na  ratio  of  .6-. 7  (16). 

Figure  6  shows  the  variation  of  Young's  elastic  modulus  measured  as  a 
function  of  composition  for  the  sodium-aluminosilicate  glass  series.  Notice 
that  the  moduli  increases  in  a  linear  manner  with  increased  A^O^ 
concentration  for  glasses  with  an  Al/Na  ratio  less  than  about  0.5.  The  moduli 
remained  constant  for  the  glasses  with  an  Al/Na  ratio  of  0.6  or  greater.  It  is 
reasonable  to  expect  that  Young's  elastic  modulus  will  reflect  structural  trends 
within  a  compositional  family  of  glasses.  Since  the  network  polymerizes  as  the 
concentration  of  non-bridging  oxygen  decreases,  this  cross-linking  should  render 
the  glass  more  resistant  to  strain  for  a  given  amount  of  stress.  Thus,  the 
observed  increase  in  Young's  modulus  with  Al/Na  is  consistent  with  a 
corresponding  decrease  in  non-bridging  oxygen  concentration.  In  this  sense,  the 
elastic  moduli  represent  an  indirect  measure  of  the  molecular  structure  of  the 
glass  network. 

It  is  evident  that  the  n-parameter  (Figure  3)  and  Young's  elastic  modulus 
(Figure  6)  exhibit  a  very  similar  dependence  upon  the  A^O^  concentration. 

This  suggests  that  the  Slow  crack  growth  behavior  and  elastic  response  of  a 
material  are  closely  related.  Figure  7  plots  the  n-parameters  as  a  function  of 
Young's  elastic  modulus  for  this  set  of  glasses,  as  well  as,  for 
soda-lime-silica  and  fused  silica  glass.  Clearly,  the  n-parameter  increases 
uniformly  with  Young's  modulus. 


YOUNG  S  ELASTIC  MODULUS 


SHEAR.  MODULUS,  G  (G 


15 


Figure  7 


A  FUSED  Si02 

■  Na  0  CqO  SiO 
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The  relationship  between  the  n-parameter  and  Young's 
elastic  modulus  in  Na  ,0 • xA 1  ,0  •  ( 3-x) Si0o  glasses,  as 
well  as,  in  fused  silica  and  soda-lime-silica. 
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Elastic  Strain  and  Slow-Crack  Growth  Susceptability 

Although  the  underlying  mechanism  of  slow  crack  growth  is  a  chemical 
reaction  between  environmental  species  and  crack  tip  bonds,  the  correlations 
developed  in  the  previous  sections  of  this  paper  suggest  that  the  rate  of  that 
chemical  reaction  is  controlled  through  the  elastic  strain.  In  this  section, 
the  expected  effect  of  strain  upon  the  rate  of  chemical  reaction  is  discussed, 
and  the  elastic  properties  are  used  to  define  the  extent  to  which  that  rate  can 
be  altered  under  the  action  of  an  applied  stress.  In  this  manner,  a 
relationship  will  be  developed  between  the  observed  slow  crack  growth  behavior 
and  the  stress-strain  behavior. 

The  rate  limiting  chemical  reaction  in  slow  crack-growth  is  expected  to 

occur  between  the  aqueous  environment  and  the  highest  strength  bonds  in  the 

glass  structure.  Of  course,  this  assumes  the  material  to  be  homogeneous  and 

isotropic  with  no  easy  fracture  path.  Sodium  ions  are  bound  to  either  a 

non-bridging  oxygen  or  an  (AIO^)  aluminate  tetrahedral  unit.  These  bonds  - 

especially  those  at  the  aluminate  sites  -  are  of  relatively  low  strength.  In 
3+  . 4+  . 

contrast,  the  A1  and  Si  ions  are  strongly  bonded  to  oxygen  to  form  the 
network  or  structural  backbone  of  the  glass.  A  simple  estimate  of  the  bond 
strength  can  be  obtained  from  the  charge/coordination  ratio.  Silicon  ions  exist 
only  in  fourfold  coordination  whereas  aluminum  ions  may  exist  in  either  four-or 
sixfold  coordination  depending  upon  the  composition.  The  bond  strength  can  be 
estimated  as : 


silicon 

(CN=4) 

4/4 

=  1 

(12) 

aluminum 

(CN=4) 

3/4 

=  0.75 

(13) 

aluminum 

(CN=6) 

3/6 

=  0.50 

(14) 

where  CN  =  the  coordination  number.  Therefore,  it  would  be  predicted  that  the 
rate  limiting  reaction  would  be  between  water  and  the  silicon-oxygen  bonds. 

In  silicate-based  glasses,  then,  the  crack-growth  reaction  of  primary 
interest  is: 


(ESi-O-SiE  )  +  H2o  t0t->  2(=  Si-OH) 


(15) 


where  E  =  activation  energy.  It  is  possible  to  express  reaction  15  by  the 
following  series  of  sequential  reactions: 


.  1  -  .+ 

=  Si— 0— Si  =  _*  =  Si  -  0  +  Si  , 

E2  + 

H20  H  +  OH  , 

+  E3  . 

=  Si-0  +  H  H  Si-OH,  and 


= Si+  +  OH 


Si-OH 


(16) 

(17) 

(18) 

(19) 


The  energy  necessary  to  break  a  silicon-oxygen  bond  is  given  by  E^,  the 
dissociation  energy  of  a  water  molecule  is  given  by  E2>  and  the  energy  gains 
due  to  satiation  of  the  dangling  surface  charges  by  the  dissociated  environment 
are  given  by  and  E^.  Since  reaction  15  is  equivalent  to  the  sum  of 
reactions  16  through  19,  the  activation  energy,  E  t,  can  be  represented  by  an 


analogous  sum: 


The  reaction  rate  can  be  approximated  by: 


expt-Etot/RT]  (21) 

and  then,  substituting  Equation  20  into  Equation  21  yields: 

rate  «  exp  [-(E^  +  E2  +  E^  +  E^)/RT]  (22) 

Assuming  the  velocity  of  crack  extension  to  be  limited  by  the  rate  at  which  the 
reaction  can  proceed,  the  crack  velocity  can  be  approximated  by: 

V  -  Vq  exp  [-(E1  +  E2  +  E3  +  E4)/RT]  (23) 

The  energy  input  to  the  crack  tip  system  lowers  the  required  activation 
energy.  This  energy  may  be  input  in  the  form  of  mechanical  energy,  thermal 
energy,  etc.  In  the  slow  crack  growth  phenomenon,  the  effect  of  mechanical 
energy  on  the  system  is  of  most  concern.  Since  the  presence  of  H+  and  OH 
species  can  be  expected  to  varying  degrees  under  all  conditions , the  application 
of  external  stress  is  of  little  consequence  in  affecting  the  ability  of  the 
environment  to  dissociate  or  satiate  the  surface  charges.  Rather,  it  is  the 
effect  of  an  external  stress  on  the  energy  necessary  to  break  a  silicon-oxygen 
bond,  Ej,  which  is  of  primary  interest.  This  concept  is  unique  relative  to 
recent  suggestions  (17),  but  is  more  consistent  with  all  of  the  experimental 
observations.  That  is,  in  isostructural  network  glasses  (e.g.,  SiC>2  and 
Na  O’ A1  _0..*  2SiO,, ) ,  as  well  as  glasses  with  widely  varying  corrosion 


susceptibilities  (e.g.,  NajO "xAljOy ( S-xJsiOj) ,  the  stress  corrosion 

susceptibility  is  relatively  insensitive  to  the  H+  or  OH-  activity  (except  in 

the  low  Kj  threshold/plateau  regions). 

Initially,  consider  an  isolated  silicon-oxygen  bond  where  the  strain  is 

simply  a  change  in  the  bond  length,  r,  relative  to  the  equilibrium  position, 

r  .  The  strain,  e,  in  this  bond  can  be  defined  as: 
o 

e  -  [r  -  r  ]/r  =  Ar/r  (24) 

o  o  o 


If  the  strain  is  tensile,  Ar  is  greater  than  0,  and  the  bond  length  will 
increase . 

The  silicon-oxygen  bond  is  predominantely  ionic  in  nature,  and  so,  the 
energy  of  the  lone  ion  pair  is  given  by  the  sum  of  the  coulombic  attractions, 
and  the  repulsion  and  ionization  energies  (18): 


[(Z.  Z0e2)/(4ire  r)]  +  [B/r”]  +  E. 
12  o  ion 

where:  Z. 

l 

e 

t 

o 
r 

B ,  n 

E. 
ion 


=  charge  on  ion  i 
=  electronic  charge 
=  permittivity 

=  ion-ion  separation  distance 
=  constants 
=  ionization  energy 


(25) 


This  relationship  is  shown  graphically  in  Figure  8.  The  equilibrium  separation 


4 


The  energy  to  rupture  a  silicon-oxygen  bond  can  be  defined  as: 


E.  *  E.  -  E 
1  b  o 


(E,  and  E  are  defined  in  Figure  8)  whereas  the  energy  required  to  rupture  a 

D  O 


strained  bond  can  be  defined  by: 


V  -Eb  "Ea 


After  substituting  Equation  26  into  Equation  27: 


V  ■  (El  *  V  -  Ea 


E.  ‘  =  E,  -  AE 
1  1  a 


where:  &E  =  change  in  rupture  energy  under  the  action  of  an  applied  stress. 

a 

When  the  bond  is  strained  from  r^  to  r,  the  repulsion  term  -  which  is 
proportional  to  r  n  -  drops  quickly  to  zero  with  respect  to  the  attraction 
term.  Therefore,  this  term  can  be  ignored  without  imparting  a  significant  error 
in  the  analysis.  The  change  in  energy,  a^  , can  now  be  defined  as: 


aE  =  [E,  -  (k/r  )]  -  [El  -  (kr)] 


3 


where:  k  =  [ Z^Zje^]/ [ 4 it ] .  Equation  31  can  be  rewritten  as: 

AE  -  k[ (r-r  )/rr  ],  (32) 

a  o  o 

and  then  substituting  Equation  24  into  Equation  32  yields: 

AEa=k[e/r]  (33) 

Equation  33  shows  quantitatively  that  as  the  strain  (or  stress)  is  increased, 
the  change  in  bond  rupture  energy  increases.  The  magnitude  of  the  increase  is 
determined  directly  by  the  amount  of  strain. 

This  exercise  is  a  simplified  manner  of  showing  that  it  is  necessary  to  do 
work  on  a  bond  in  order  to  strain  it  and  that  the  work  is  stored  in  the  bond  and 
used  to  reduce  the  effective  bond  rupture  energy.  Elasticity  theory  defines  the 
amount  of  work  done  on  the  bond  by: 


o 


where  Y  =  stored  elastic  strain  energy.  Thus,  the  next  step  in  the  derivation 
is  to  define  the  stored  elastic  strain  energy  at  the  crack-tip  in  terms  of  the 
elastic  constants. 

For  the  isolated  ion  pair,  the  stress  necessary  to  strain  a  bond  from  0  to  e 
is  uniquely  defined  by  the  properties  of  the  bond,  whereas  the  bond  of  interest 
in  slow  crack-growth  is  bound  within  an  infinite  glass  network.  The  amount  of 
stress  necessary  to  strain  that  bond  from  0  to  e  is  no  longer  uniquely  defined 
by  the  bond,  but  is  determined  by  the  material  surrounding  the  bond.  In  other 
words,  the  amount  of  stored  elastic  strain  energy  available  is  determined  by  the 
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elastic  properties  of  the  surrounding  material  and  not  the  properties  of  the 
bond  itself. 

It  is  reasonable  to  assume  that  the  stored  elastic  strain  energy  in  the 
vicinity  of  the  crack  tip  is  available  to  reduce  the  effective  bond  rupture 
energy.  The  work  necessary  to  produce  a  strain  can  be  expressed  as  (19): 

Y  =  ( l/2)c . j  e .  t  ^  (35) 

where:  =  stiffness  matrix.  Considering  only  the  unidimensional  strain, 

e  2>  Equation  35  becomes: 

Y  =  cn  (  e2)2  (36) 


The  magnitude  of  the  stored  energy  is  determined  by  the  strain  and  stiffness. 

Qualitatively,  the  n-parameter  should  be  related  to  the  magnitude  of  the 
stored  energy  since  glass  can  support  only  small  strains.  Thus,  assuming  that 
the  strain  in  the  crack  tip  region  does  not  vary  significantly  from  glass  to 
glass,  the  difference  in  stored  energy  would  be  determined  primarily  by 
differences  in  the  stiffness.  Therefore,  the  n-parameter  should  vary  as  c. 


'll' 


The  elastic  moduli  can  be  used  to  express  either  the  stiffness  (s..)  or 

ij 


the  compliance  (c„).  For  an  isotropic  material  (19): 


sn  =  1/E 


(37) 


s 


12 


'>!  E 


(38) 


c,,  =  [stl  +  s ,„]/[( s ,  .  -  s,„)(s,,  +  2s, „)] 


(39) 


where:  v  =  Poisson's  ratio.  Using  these  relationships,  it  is  possible  to 

express  c^  as: 


cn  =  S  •  E  (40) 

where:  S  =  [1  —  v ] /  [ C 1  +v)(l_2v)] 

Since  the  n-parameter  is  expected  to  vary  as  c^,  and  c^  *  S*E,  then  a 
relationship  between  n  and  S *E  exists.  Figure  9  shows  the  relationship  between 
the  n-parameter  and  the  predicted  function  (Eq  40)  of  the  elastic  modulus  (S»E). 
Given  the  fact  that  the  prediction  is  semi-quantitative,  the  agreement  is  very 
good.  Thus,  one  could  re-write  the  expression  for  the  crack  velocity  as: 

V  =  V  (KT)nSE  (41) 

o  l 

where  Vq  and  n  are  constants  and  S  and  E  are  material  parameters.  Notice  that 
silica  and  soda-lime-silica  data  also  fit  the  correlation  in  Figure  9. 

Therefore,  n  may  be  a  universal  constant  for  oxide  glasses  in  which  case  the 
material  parameters  S  and  E  can  be  used  to  predict  the  composition  dependence  of 
slow  crack  growth  via  Eq  (41).  The  analysis  developed  does  provide  insight  as 
to  why  the  stress-strain  behavior  of  a  material  is  intimately  related  to  its 
slow  crack  growth  behavior.  It  leads  naturally  to  the  conclusion  that  the 
composition  dependence  of  the  stress  corrosion  susceptibility  is  due  to 
composition  dependence  of  the  elastic  properties  and  not  the  chemical  corrosion 


properties . 


PARAMETER 


S 


i 


Figure  9.  The  relationship  between  n-parameter  and  the  predicted 
function  of  the  elastic  moduli  for  Na^O’xAl^O  • (3-x)SiO 
glasses,  as  well  as,  for  fused  silica  and  soda-lime-silic. 
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CONCLUSIONS 

The  region  1  slow  crack  growth  behavior  was  characterized  by  the 
n-parameter  which  varied  from  16  to  23  and  varied  linearly  with  alumina 
concentration  for  0.0,  0.2,  and  0.4  glasses.  However,  for  0.6,  0.8,  and  1.0 
glasses  the  n-parameter  was  independent  of  alumina  concentration  and  equivalent 
to  the  n-parameter  for  the  0.4  glass.  The  invariance  of  the  n-parameter  with 
variations  in  the  environmental  pH  suggest  that  corrosion  plays  only  a  minor 
role  in  determining  the  slow  crack  growth  behavior.  It  was  proposed  that  the 
elastic  properties,  and  hence  the  glass  structure,  determines  the  n-parameter. 

The  fundamental  mechanism  of  slow  crack  growth  was  assumed  to  be  a  chemical 
reaction  between  the  glass  and  the  aqueous  environment  in  contact  with  it.  It 
was  proposed  that  the  activation  energy  of  such  a  reaction  would  be  lowered  by 
an  amount  equivalent  to  the  strain  energy  stored  in  the  crack  tip  vicinity.  It 
was  shown  that  the  n-parameter  is  directly  related  to  the  magnitude  of  this 
stored  strain  energy  which,  in  turn,  is  determined  by  the  elastic  properties. 

In  this  manner,  the  elastic  properties  were  related  directly  to  the  slow  crack 
growth  behavior.  A  good  correlation  between  the  n-parameter  and  elastic 
properties  was  observed  for  sodium-aluminosilicate,  soda-lime-silica,  and  fused 
silica  glasses. 


ACKNOWLEDGEMENT 

The  authors  gratefully  acknowledge  the  financial  support  of  the  Air  Force 


Office  of  Scientific  Research  (AFOSR-82-0013) . 


87 


REFERENCES 


1 


1.  R.  E.  Mould  and  R.  D.  Southwick,  "Strength  and  Static  Fatigue  of  Abraded 
Glass  Under  Ambient  Conditions,"  J.  Am.  Ceram.  Soc.,  42  [12],  582  (1959). 

2.  S.  M.  Wiederhorn,  "A  Chemical  Interpretation  of  Static  Fatigue,"  J.  Am. 
Ceram.  Soci,  55  [2],  81  (1972). 

3.  W.  B.  Hillig  and  R.  J.  Charles,  "Surfaces,  Stress-Dependent  Surface 
Reactions,  and  Strength,"  High  Strength  Materials,  pp.  682-701,  edited  by 
W.  F.  Zackey,  Wiley  Publishing  Co.,  New  York,  New  York,  1965. 

4.  S.  M.  Wiederhorn  and  L.  H.  Bolz,  "Stress  Corrosion  and  Static  Fatigue  of 
Glass,"  J.  Am.  Ceram.  Soc.,  53  [10],  543  (1970). 

5.  S.  M.  Wiederhorn,  "Crack  Growth  as  an  Interpretation  of  Static  Fatigue," 
Journal  of  Non-Crystalline  Solids,  19,  169  (1975). 

6.  S.  W.  Freiman,  "Fracture  Mechanics  of  Glass,"  Glass  Science  and  Technology: 
Elasticity  and  Strength  in  Glasses,  Vol.  15,  pp.  21-78,  Edited  by  D.  R. 
Uhlman  and  N.  J.  Kreidl,  Academic  Press,  New  York,  New  York,  1980. 

7.  R.  R.  Tummala,  "Stress  Corrosion  Resistance  Compared  with  Thermal  Expansion 
and  Chemical  Durability  of  Glasses,"  Glass  Technology,  17  [4],  145  (1976). 

8.  H.  M.  Cekirge,  W.  R.  Tyson  and  A.  S.  Krausz,  "Static  Corrosion  and  Static 
Fatigue  of  Glass,"  J.  Am.  Ceram.  Soc.,  59  [5-6],  265  (1976). 

9.  C.  J.  Simmons  and  S.  W.  Freiman,  "Effect  of  Corrosion  Processes  on 
Subcritical  Crack  Growth  in  Glass,"  J.  Am.  Ceram.  Soc.,  64  [11],  683 
(1981). 

10.  S.  W.  Freiman,  D.  R.  Mulville  and  P.  W.  Mast,  "Crack  Propagation  Studies  in 
Brittle  Materials,"  Journal  of  Materials  Science,  8  [11],  1527  (1973). 

11.  W.  Wasserman  and  J.  Neter,  Applied  Linear  Statistical  Models,  Richard  D. 
Irwin,  Inc.,  Honewood ,  Illinois,  1978. 

12.  D.  N.  Coon,  J.  J.  Mecholsky  and  C.  G.  Pantano,  "A  Static  Fatigue  Limit  due 
to  Crack-Healing  in  Na20«3SiC>2  Glass,"  to  appear. 

13.  Corrosion  of  Glass  by  D.  E.  Clark,  C.  G.  Pantano  and  L.  L.  Hench  (Books  for 
Industry,  New  York,  1979). 

14.  T.  M.  El-Shamy,  J.  Lewis  and  R.  W.  Douglas,  "Dependence  on  the  pH  of  the 
Decomposition  of  Glasses  by  Aqueous  Solutions,"  Glass  Technology,  B(3),  81 
(  1972). 

15.  S.  M.  Wiederhorn  and  H.  Johnson,  "Effect  of  Electrolyte  pH  on  Crack 
Propagation  in  Glass,"  J.  Am.  Ceram.  Soc.,  56  [4],  192  (1973). 


ir.'-v- 


I. 

► 

L.' 


A 


fcr 


■ 


.V.VA^Vy-tV'^Tyi.11 '.".«■  ,wTs.»  Jj.f  v.'M.  j*  ■<  ijw  y  j  ^gfgr  g«  gag? 


88 


16.  B.M.J.  Smets  and  T.P.A.  Lommen,  "The  Incorporation  of  Aluminum  Oxide  and 
Boron  Oxide  in  Sodium  Silicate  Glasses,  Studied  by  X-ray  Photoelectron 
Spectroscopy,"  Physics  and  Chemistry  of  Glasses,  22  [6],  158  (1981). 

17.  T.  A.  Michalske  and  S.  W.  Freiman,  "A  Molecular  Interpretation  of  Stress 
Corrosion  in  Silica,"  Nature,  295  [11],  511  (1982). 

18.  W.  D.  Kingery,  H.  K.  Bowen  and  D.  R.  Uhlmann,  Introduction  to  Ceramics, 
Wiley  Co.,  New  York,  New  York,  1976. 

19.  J.  F.  Nye,  Physical  Properties  of  Crystals,  Oxford  University  Press,  1979. 


■  -V. 


1-25 


89 


CHAPTER  III 


SLOW  CRACK  GROWTH  BEHAVIOR  IN  THE  THRESHOLD  REGION 


D.  N.  Coon,  J.  J.  Mecholsky  and  C.  G.  Pantano 
Department  of  Materials  Science  and  Engineering 
Pennsylvania  State  University 
University  Park,  PA  16802 


Submitted  for  publication  to  J,  American  Ceramic  Society. 


*  *  ■  ’  »  '  .  "  fc  *  a  *  ■  *  a'*  »  *  *  •  a^*  kN  ‘v  a*»  , 


90 


Slow  Crack  Growth  Behavior  in  the  Threshold  Region 

by 

D.  N.  Coon,'*'  J.  J.  Mecholsky  and  C.  G.  Pantano 
Department  of  Materials  Science  and  Engineering 
Pennsylvania  State  University 
University  Park,  PA  16802 

ABSTRACT 

Crack  velocity  was  measured  as  a  function  of  stress  intensity  for 

Na^O -3Si02  glass  in  several  solutions  varying  in  pH  from  1  to  11.  At 

1/2  —8 

pH  =  11,  a  velocity  plateau  was  observed  at  K^<.35  MPa’m  and  v  s10 

1/2  —9 

m/sec;  in  a  solution  of  pH  5,  the  plateau  shifted  to  K<.30  MPa*m  '  and  v<10 

m/sec.  When  AlCl^  was  added  to  the  solution,  the  plateau  disappeared  and  a 

1/2 

threshold  was  observed  at  K=.40-.45  MPa*m  .  Similarly, 

sodium-aluminosilicate  glasses  exhibited  a  threshold  rather  than  a  velocity 
plateau . 

This  behavior  is  explained  in  terms  of  interactions  between  the  fracture 
surfaces  wherein  ionized  surface  sites,  specifically  adsorbed  aluminum  cations, 
or  hydrogen  bonds  create  forces  of  repulsion  and  attraction.  These  observations 
represent  the  first  experimental  evidence  that  processes  behind  a  propagating 
crack-tip  can  dramatically  affect  the  crack  extension  rate.  Since  these  effects 
are  limited  to  the  low  velocity/low  stress  intensity  regime,  they  are  relevant 
to  the  practical  strength  and  fatigue  of  glass. 

Now  with  EG&G  Idaho,  Idaho  Falls,  ID. 
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INTRODUCTION 

Wiederhorn  (l)  has  proposed  chat  Region  1  slow  crack  growth  can  be 

explained  by  a  stress  enhanced  chemical  reaction  between  the  environment  and 

glass.  These  theories  predict  that  the  log  of  the  crack  velocity  varies 

linearly  with  the  applied  stress  intensity.  However,  it  has  also  been  reported 

that  slope  changes  are  observed  at  low  stress  intensities  (2).  The  objective  of 

this  study  was  to  examine  the  slow  crack  growth  behavior  of  a  simple  NajOOSiOj 

glass  at  low  stress  intensity,  ,  and  to  determine  the  effects  of  specific 

+  3+ 

adsorbates  -  namely  H  and  Al  -  upon  the  presence  or  absence  of  velocity 
plateaus  and  thresholds.  The  crack  growth  characteristics  of  the  binary 
sodium-silicate  glass  in  the  environment  containing  aluminum  was  compared  to  the 
behavior  of  sodium-aluminosilicate  glasses  in  water. 

BACKGROUND 

Simmons  and  Freiman  (2)  were  the  first  to  report  the  existence  of  low 

velocity  plateaus  in  the  V-K^  diagram  for  glass.  They  observed  these  plateaus 

in  binary  soda-silica  glasses  at  K^  levels  which  depended  upon  pH,  although  they 

1/2 

were  typically  in  the  range  .30-. 40  MPa-m  .  These  investigators  suggested 
that  the  velocity  plateaus  were  due  to  corrosion;  i.e.,  the  crack  grows  at  a 
velocity  which  approximates  the  rate  of  silica  dissolution.  Nevertheless,  their 
data  showed  that  the  plateau  shifted  to  lower  velocity  with  increasing  pH 
(whereas  the  corrosion  rate  increases  with  increasing  pH).  Moreover,  the 
plateau  occurred  at  the  same  velocity  in  both  25%  Na^CWSZ  SiO^  and  33% 

Na2~67%  SiO^  glasses  even  though  the  corrosion  rates  of  these  glasses  are  very 
different  at  any  pH. 

More  recently,  Lawn  and  co-workers  (3,4)  have  attempted  to  explain  slow 
crack  growth  on  the  basis  of  mechanisms  occurring  behind,  and  not  at,  the  crack 
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tip.  They  argue  that  environmental-fracture  surface  interactions  can  cause  the 
neutralization  of  charged  surfaces,  or  conversely,  the  charging  of  neutral 
surfaces.  Electrostatic  interactions  result  when  two  charged  surfaces  are  in 
close  proximity  to  one  another.  Surfaces  of  like  charge  would  repel,  whereas 
surfaces  in  the  presence  of  electrolyte  species  could  form  adsorbed  double 


I 

» 

I 

► 

l* 

i 


layers  whose  overlap  or  cationic/anionic  bridging  could  lead  to  attractive 
forces.  Since  crack  surfaces  are  separated  by  only  a  small  distance  at  low  K 
the  authors  suggested  the  existence  of  large  electrostatic  interactions  between 
the  crack  walls.  They  have  predicted  three  distinct  types  of  stress 
intensity-crack  velocity  behavior  depending  on  the  electrical  nature  of  the 
crack  surfaces  (Figure  l).  There  are  two  features  of  their  model  which  are 
relevant  here.  First,  at  low  levels,  constant  velocity  plateaus  can  be 
explained  by  electrostatic  interactions.  Second,  above  some  level,  shown  as 
KA  in  Figure  1,  the  slow  crack  growth  behavior  becomes  independent  of  the 
electrical  nature  of  the  crack  surfaces.  This  behavior  at  higher  is  in 
agreement  with  the  work  of  Wiederhorn  and  Johnson  (5).  They  measured  the  Region 
I  slow  crack  growth  behavior  of  soda-1 ime-s i 1 ica  glass  in  various  aqueous 
environments  which  produced  surfaces  of  known  -  but  varying  -  zeta  potential. 

The  zeta  potential  was  assumed  to  be  a  measure  of  the  electrical  charge  at  the 
glass  surface.  They  observed  no  effect  of  zeta  potential  on  slow  crack  growth 
behavior.  This  suggests  that  at  velocities  in  excess  of  10  ^  m/s,  the  slow 
crack  growth  behavior  of  glass  is  independent  of  electrostatic  crack  wall 
interactions.  However,  electrostatic  effects  may  be  important  at  lower 
velocities  than  were  considered  by  Wiederhorn  and  Johnson  (5). 

The  possibility  of  estimating  the  relative  surface  charge  due  to  ionized 


hydroxyls  on  a  silicate  glass  surface  has  been  reported  (6).  The  ionization 
process  is  represented  bv: 


STRESS  INTENSITY 


Figure  1 


Schematic  of  crack  velocity  as  a  function  of  stress 
intensity  for  cracks  with  interacting  surfaces.  The 
threshold  behavior  (A)  represents  fracture  surfaces 
with  an  attractive  interaction,  whereas  the  plateau 
(C)  is  due  to  a  repulsive  interaction.  The  extension 
of  Region  I  -  (B)  -  occurs  when  surface  forces  are 
insufficient  to  influence  the  crack  growth. 
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SiOH  £  SiO-  +  H+, 


(1) 


in  which  case  the  fraction  of  hydroxyls  ionized  as  a  function  of  pH  is  given  by: 


[SiO~ 1  . rt-9. 8  +  pH 

TsTohT  =  10 


(2) 


These  calculations  are  given  in  Table  I.  A  high  pH  environment  produces  a 
surface  of  high  net  negative  charge.  It  is  in  this  sense  that  the  surface 
charge  of  a  silicate  glass  is  pH  dependent.  Of  more  interest  here  is  the  fact 
that,  in  principle,  the  surface  charge  in  a  crack  can  be  varied  through  control 
of  the  solution  pH  and/or  the  solution  composition  during  a  crack-growth 
experiment . 

The  existence  of  negatively  charged  silica  surfaces  in  solutions  of  high  pH 
has  been  experimentally  verified  using  streaming  potential  and  other  methods 
(7).  These  studies  also  showed  that  the  presence  of  multivalent  cations  in 
solution,  for  example  trivalent  aluminum  (Al+^),  can  neutralize  or  even  reverse 
the  charges. 

EXPERIMENTAL  PROCEDURE 

Specimens  of  Na20*3Si02  glass,  as  well  as  a  series  of  sodium 
aluminosilicate  glasses  (Na  0- xAl .0,- ( 3-x)Si0_) ,  were  cut  into  plates 
approximately  7.5  cm  x  1.27  cm  x  0.25  cm  with  a  groove  down  the  center  of  the 
plate.  The  groove  depth  was  approximately  one-half  of  the  plate  thickness  and 
the  groove  width  was  0.00075  to  0.001  cm.  Slow  crack  growth  behavior  was 
determined  by  the  constant  moment  double  cantilever  beam  technique  as  described 
by  Freiman,  Mullville,  and  Mast  (8).  The  slow  crack  growth  behavior  of  these 
specimens  was  examined  in  water  (pH  5),  HC1  solution  (pH  1),  NaOH  solution  (pH 
11),  and  a  NaOH  solution  to  which  0.4  g/1  of  AlClj*6H^0  had  been  added. 
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RESULTS 

Figure  2  shows  the  crack  velocity  data  sets  of  the  NajO'SSiC^  glass  tested 

1/2 

in  solutions  of  pH  1,  5,  and  11.  Above  stress  intensities  of  0.35  MPa*m  '  and 
-8 

velocities  of  10  m/s,  no  difference  in  the  slow  crack  growth  behavior  was 
observed.  All  measurements  simply  scatter  about  the  line  predicted  by  the  data 
obtained  in  water.  At  low  velocity,  a  stress  independent  plateau  was  observed 

_o 

in  solutions  of  pH  5  and  11.  The  plateau  was  observed  at  a  velocity  of  10  m/s 

in  the  pH  11  solution  and  10  m/s  in  the  pH  5  solution.  No  plateau  was 

observed  at  pH  1 .  The  pH  dependence  of  the  V-K^  diagrams  obtained  here 

contradict  those  reported  by  Simmons  and  Freiman  (2),  but  interestingly,  these 

data  are  quite  consistent  with  their  corrosion  model. 

Figure  3  shows  two  sets  of  V-K^  data.  One  was  observed  for  NajO’SSiC^ 
tested  in  pH  11  solution  while  the  other  was  observed  in  a  pH  11  solution  to 
which  A1C1^«6H20  had  been  added  prior  to  testing.  In  the  solution  containing 
A1  ions,  no  plateau  was  observed. 

The  series  of  sodium-aluminosilicate  glasses  [Na^O* xAl jO^ ' ( 3-x)Si02 
with  x=0.2,  0.4,  0.6,  0.8,  and  1.0],  as  well  as  commercial  borosilicate  and 
soda-lime  glasses  [which  contain  a  fraction  or  up  to  2  percent  Al^O^]  do  not 
exhibit  low  velocity  plateaus  under  any  conditions  (2,9).  Figure  4  presents  the 
crack  velocity  sets  for  the  Na.,0  * .  2A 1  *  2 . 8S  i  0,,  and  Na.,0 '  A1  2O ^  *  2S  i  C>2 

glasses  at  pH  11.  These  should  be  compared  with  the  V-K  d  i.igr  ims  in  Figure  3. 

DISCUSSION 

The  surface  charge  of  silicate  glasses  is  known  to  be  r>H  dependent.  It  has 
already  been  pointed-out  that  the  surface  is  predominantly  negative  m  sol  it  ions 
of  high  pH,  but  that  the  presence  of  Al  *  ^  ions  in  s.>!  »t  ;  »n  •  a  1  '  .  -  r  :  ne  sut  f  -ice 
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charge.  First,  the  effect  of  solution  pH  upon  the  low  velocity  plateaus  will  be 
discussed,  and  then  the  effects  of  A1  ions  in  solution  are  considered. 

The  slow  crack  growth  behavior  of  Na2o  *3Si02  glass  was  determined  in 
solutions  of  pH  1 ,  5,  and  11  (Figure  2).  Above  stress  intensities  of 

1  j  n 

approximately  0.35  MPa*m  '  ,  no  effect  of  solution  pH  on  crack  growth  was 

1/2 

observed.  However,  at  stress  intensities  below  0.35  MPa*m  '  ,  the  behavior  of 
the  material  changed  dramatically  with  pH.  A  constant  velocity  plateau  was 
observed  when  the  experiments  were  performed  in  pH  5  and  pH  11.  The  plateau 
occurred  over  a  larger  stress  intensity  range  and  was  observed  at  a  higher 
velocity  in  the  pH  11  solution.  The  plateau  occurred  over  a  smaller  range  of 
stress  intensities  and  was  observed  at  lower  velocities  in  water.  When  the 
solution  pH  was  lowered  to  1,  the  plateau  was  not  observed. 

The  Na2o*3Si02  glass  is  known  to  chemically  alter  the  aqueous  solution  in 

contact  with  it  (6).  That  is,  the  alkali  is  released  into  solution  and  thereby 

increases  the  pH.  The  pH  usually  increases  up  to  a  point  termed  'the 

equilibrium  pH'  (pH^),  but  the  pHg(j  value  depends  upon  the  glass  composition 

and  surface-area  of  glass  to  solution-volume  ratio.  Obviously,  the  pH^  will  be 

achieved  quite  rapidly  within  the  confines  of  the  crack-tip,  especially  in  the 

low  velocity  regions.  The  equilibrium  pH  of  Na20‘3Si02  has  been  reported  (9), 

and  using  that  data,  one  can  predict  that  the  crack  tip  solution  will  increase 

in  pH  by  at  least  4  units  for  pH  =1  and  pH^  a  5  (where  pH  is  the  initial  or 

o  o  o 

'bulk'  solution  pH).  Since  pHQ  =  11  is  well  above  the  equilibrium  pH,  however, 
no  significant  change  in  the  crack  tip  pH  was  anticipated  under  these  alkaline 
conditions.  Using  these  predictions  of  the  crack-tip  pH  (5,  9  and  11  for 
solution  pHQ  values  of  1,  5  and  11,  respectively),  one  can  estimate  the 
fraction  of  ionized  surface  sites  at  the  crack-tip.  A  crack  tip  solution  of  pH 
11  would  result  in  a  surface  whose  silanol  species  are  completely  ionized.  In 


this  case,  both  fracture  surfaces  would  be  highly  negative  and,  there  would  be  a 
'strong'  electrostatic  repulsion  (although  the  force  is  clearly  insufficient  to 
propagate  the  crack  at  ^<.15  MPa.ra  ).  A  crack  tip  solution  of  pH  9  (pH  5 
bulk  solution)  would  result  in  a  surface  where  the  silanols  are  only  partially 
ionized.  In  this  case,  the  fracture  surfaces  would  possess  a  small  negative 
charge.  There  would  be  a  'weaker'  electrostatic  repulsion  and  thus,  the  plateau 
shifts  to  lower  velocity  and  the  threshold  moves  to  higher  K^.  a  crack  tip 
solution  of  pH  5  would  result  in  a  surface  which  is  almost  completely 
hydroxy lated .  Both  fracture  surfaces  would  be  neutrally  charged  and  no 
repulsive  forces  will  exist.  Nevertheless,  there  may  be  an  interaction  between 
these  neutrally-charged,  but  fully  hydroxylated  surfaces.  In  fact,  it  is 
possible  that  the  threshold  observed  here,  as  well  as  in  the  other  glasses,  is 
due  to  hydrogen  bonding  and/or  double-layer  overlap  behind  the  crack-tip  which 
creates  an  attractive  force  between  the  fracture  surfaces.  These  thresholds  are 
usually  observed  in  glasses  containing  modifier  species  where  surface  silanol 
densities  and/or  double-layer  formations  are  expected  to  be  most  extensive.  It 
is  likely  that  the  attractive  force  due  to  hydrogen  bonding,  in  particular,  is 
significant  only  when  ionic  forces  are  minimized  or  when  very  low  conditions 
bring  the  fracture  surfaces  into  close  proximity.  In  this  sense,  the  thresholds 
observed  at  low  may  be  due  _  in  general  -  to  hydrogen-bonding  between 

surface  silanols.  The  fact  that  a  threshold  (or  plateau)  has  not  been  reported 
for  silica  can  be  attributed  to  the  low  surface  silanol  density  of  the  fully 
cross-linked  SiO^  glass;  thus,  Region  I  crack  growth  extends  to  very  low  Kj . 

Any  plateaus  and  thresholds  for  silica  glass,  if  they  exist,  would  occur  at  a 
crack  velocity  below  the  detection  limit  of  most  experimental  measurements. 

The  low  velocity  plateau  was  not  observed  when  the  Na20*3Si02  glass  was 
tested  in  a  solution  containing  lCT^  moles/1  A1C13»6H20;  the  aluminum  chloride 


was  added  directly  to  the  pH  11  solution  in  order  to  keep  the  alkali  ionic 
strength  a  constant.  The  pH  of  this  solution  dropped  to  5  due  to  addition  of 
the  aluminum  chloride  species.  Nevertheless,  according  to  Horn  and  Onoda  (7), 
the  presence  of  this  aluminum  species  at  10“^  M/1  and  pH  5  changes  the  zeta 
potential  from  -100  mV  to  125  mV;  alone,  the  pH  change  from  11  to  5  would  reduce 
the  value  of  the  zeta  potential  slightly  {“20  mV),  but  could  not  reverse  the 
sign.  The  elimination  of  the  low  velocity  plateau  with  AlCl^gi^O  addition  - 
and  the  consequent  threshold  behavior  -  is  consistent  with  the  concept  that  the 
plateaus  result  from  negatively-charged  fracture  surfaces;  i.e.,  Al+^  can 
neutralize  and/or  reverse  the  charge  by  specific  adsorption  upon  the  fracture 
surfaces  at  and  behind  the  crack-tip.  It  is  likely  that  the  presence  of  a 
threshold  under  these  conditions  is  due  to  Al+^  'cation  bridging'  across  the 
crack  and/or  overlap  of  the  double-layers  associated  with  the  adsorbed  aluminum 
species;  in  a  sense,  there  is  little  fundamental  difference  between  these  two 
effects.  These  data,  and  their  interpretation,  are  clearly  consistent  with  the 
recent  models  put  forth  by  Lawn  and  co-workers  (3,4). 

It  is  especially  noteworthy  that  the  sodium  aluminosilicate  glasses  -  at 
pH  11  -  do  not  exhibit  the  low  velocity  plateaus  expected  under  these 
conditions.  Thus,  it  seems  likely  that  the  aluminum  in  the  glass  plays  a 
dynamic  role  during  the  crack  growth  process  in  the  threshold  regime.  The 
absence  of  low  velocity  plateaus  in  the  sodium-aluminosilicate  glasses  (Figure 
4),  as  well  as  in  commercial  glasses  containing  alumina  and  probably  in  any 
glasses  containing  multivalent  cations,  may  be  due  to  'steady-state' 
dissolution/adsorpt ion  phenomena  at  and  behind  the  crack-tip.  There  is  always 
some  dissolution  of  the  glass  within  the  crack,  and  thus  the  presence  of  Al+^  or 
other  multivalent  ions,  can  be  expected  even  when  they  are  not  intentionally 


introduced  into  the  crack-growth  environment.  This  represents,  perhaps,  the 
only  practical  approach  to  chemical  passivation  of  slow  crack  growth  in  glass. 

CONCLUSIONS 

The  occurrence  of  low  velocity  plateaus  and/or  thresholds  in  sodium 
silicate  glasses  was  shown  to  be  dependent  upon  environmental  pH  and  the 
presence  of  Al3+  ions  in  solution.  This  behavior  was  shown  to  be  consistent 
with  an  interaction  between  the  fracture  surfaces .  Such  behavior  represents  the 
first  experimental  evidence  that  processes  behind  a  propagating  crack  tip  can 
dramatically  affect  the  crack  extension  rate. 
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ABSTRACT 

A  model  is  developed  which  can  explain  the  static  fatigue  limit  in  terms  of 
a  crack  healing  mechanism.  The  model  assumes  that  a  high  water-content  material 
can  form  at  the  crack-tip  due  to  healing  when  K^XK^  (the  static  fatigue 
limit).  This  'healed'  region  is  resistant  to  crack  propagation  not  because  it 
is  blunted,  but  rather,  because  it  is  deformable.  In  this  study,  the  model  was 
tested  by  aging  Na20-3SiC>2  glass  at  K^Kj  ,  and  then  measuring  the  time 
delay  to  crack  extension  upon  reloading  at  Kj>Kjq.  The  time  delays  could  be 
predicted  as  a  function  of  reload  stress  intensity  by  modeling  the  strain  rate 
in  terms  of  a  viscous  deformation  and  by  assuming  the  existence  of  some  critical 
strain  for  crack  re-propagation. 
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INTRODUCTION 

Originally,  crack  extension  studies  were  envisioned  as  a  route  to  the 
prediction  of  service  lifetimes  of  structural  materials.  The  existence  of  a 
static  fatigue  limit,  which  is  a  stress  level  below  which  cracks  do  not  extend, 
would  be  of  fundamental  importance  in  the  design  of  these  materials.  Typically, 
a  static  fatigue  limit  is  observed  in  crack  extension  studies  of  glass  at  stress 
levels  approximately  one-third  the  fracture  toughness.  Clearly,  the 
determination  of  the  mechanism  by  which  this  fatigue  limit  occurs  has  been  the 
stimulus  for  much  experimental  activity. 


BACKGROUND 

The  two  theories  which  are  most  widely  accepted  for  explaining  the  static 
fatigue  limit  in  glass  involve  crack-tip  blunting  (1,2)  and  crack  healing  (3,4). 
More  recently,  it  has  been  suggested  that  interactions  between  the  fracture 
surfaces  behind  the  crack  tip  may  contribute  to  the  fatigue  limit,  or 
'threshold,'  for  crack  growth  (5,6).  The  blunting  and  healing  models  are  most 
pertinent  to  the  arguments  presented  in  this  study. 

Charles  and  Hillig  (1)  suggested  that  the  driving  force  for  blunting  could 
be  envisioned  as  being  composed  of  two  contributing  factors.  First,  as  the 
crack-tip  radius  of  curvature  increases  (i.e.,  the  tip  becomes  more  blunt),  the 
driving  force  for  blunting  decreases.  Secondly,  as  the  corrosiveness  of  the 
environment  increases  (i.e.,  the  tendency  of  the  environment  to  attack  the 
material  increases),  the  driving  force  for  blunting  increases.  They  modeled 
these  trends  as: 

Ebu„,  ■  (VI/S  <» 


109 


where  E  *  driving  force  for  blunting 

o  lunt 

D  3  crack  tip  radius 

V  =  molar  volume  of  the  material  being  corroded 
m 

Y  3  surface  energy  of  the  material  with  respect  to  the  environment. 

An  externally  applied  stress  would  provide  a  driving  force  for  crack  extension 
which  would  act  in  opposition  to  the  crack  tip  blunting  force.  This  driving 
force  would  be  high  for  large  values  of  applied  stress  and  small  crack  tip 
radii  . 

A  simple  energy  balance  of  the  driving  forces  determines  whether  the  tip 
will  blunt  or  extend.  When  the  two  driving  forces  are  exactly  equal,  the  crack 
tip  will  neither  extend  nor  blunt.  This  is  traditionally  termed  the  static 
fatigue  limit,  kIq.  For  all  values  of  greater  than  KIo,  the  driving 
force  for  extension  is  greater  than  the  driving  force  for  blunting,  and  the 
crack  tip  extends  at  some  finite  rate  maintaining  a  constant  tip  radius.  In 
contrast,  for  values  of  less  than  Kj0,  the  driving  force  for  blunting  is 
greater  than  the  driving  force  for  extension,  and  the  crack  tip  blunts. 

However,  this  blunting  will  not  continue  indefinitely  to  produce  a  crack-tip  of 
infinite  radius.  Both  driving  forces  are  dependent  upon  the  tip  radius.  As  the 
crack-tip  radius  increases,  the  driving  force  for  blunting  decreases  by  p  ^ 


while  the  crack  extension  driving  force  decreases  by  p 


The  two  driving 


forces  will  become  equal  at  some  p'  (which  is  greater  than  the  original  radius) 
and  then,  the  crack  extension  ceases. 

Michalske  (2)  used  the  model  described  above  to  explain  the  delay  in  crack 
extension  which  is  observed  after  aging  a  crack-tip  under  load.  He  showed  that 
double  cantilever  beam  specimens  of  soda-lime-silica  glass  aged  at  stress  levels 
below  the  fatigue  limit  exhibited  a  hysteresis  upon  reloading  to  levels  above 
the  fatigue  limit.  He  hypothesized  that  aging  specimens  at  a  of  0.225 


no 


1/2 

MPa*m  ,  which  is  less  than  the  observed  static  fatigue  limit  of  approximately 
1/2 

0.3  MPa'm  ,  would  result  in  crack  tip  blunting.  He  then  showed  that  when 
this  specimen  was  reloaded  to  levels  where  measurable  crack  extension  had 
been  observed  before  aging,  the  cracks  did  not  extend  without  an  appreciable 
time  delay.  This  time  delay  was  attributed  to  the  time  necessary  to  sharpen  a 
crack  tip  which  had  blunted  during  aging.  He  further  supported  his  model 
through  fractographic  examination  of  the  aged  fracture  surfaces.  To  date,  the 
work  of  Michalske  represents  the  strongest  evidence  of  crack  tip  blunting  in 
glass  at  room  temperature. 

Fuller  and  Thomson  (3)  have  developed  an  atomistic  model  for  crack 
extension.  They  assume  a  two-dimensional  lattice  of  two  parallel  infinite 
chains  of  atoms.  The  separation  of  these  two  chains  of  atomically  bonded  atoms 
is  used  to  describe  crack  growth.  Their  derivation  of  the  cohesive  energy  is  an 
expression  containing  three  terms: 


cohes ive 


2yv 


2^  l 

j=n+l 


where  n  35  the  number  of  bonds  in  the  chain 
c  *  the  interatomic  distance 

Y  *  the  surface  energy  or  fracture  surface  energy 

u  =  the  displacement  distance  of  a  given  atom  from  its  equilibrium  position 
x  =  constant . 

The  first  term  simply  accounts  for  the  energy  of  formation  of  two  new  surfaces 
where  Yg  is  the  fracture  surface  energy  per  unit  length  and  nc  is  the  crack 
length.  The  second  term  accounts  for  the  strain  energy  stored  in  the  crack-tip 
bond,  which  is  the  only  bond  considered  to  undergo  a  non-linear  deformation. 

The  third  term  accounts  for  the  strain  energy  in  all  the  bonds  ahead  of  the 


crack  tip  bond.  These  bonds  support  linear  deformation  only.  After  performing 
an  energy  balance.  Fuller  and  Thomson  showed  that: 

crack  extension  force  =  £  restoring  forces  (3) 

The  crack  extension  force  was  taken  to  be  proportional  to  K^.  Thus,  at  K^'s 
greater  than  K+  (where  K+  is  analogous  to  K^),  the  crack  propagates 
catastrophically  since  crack  extension  would  result  in  a  lowered  crack  driving 
force.  For  K^'s  between  K+  and  K  ,  the  crack  can  exist  at  two  possible 
lengths  corresponding  to  equivalent  driving  forces.  In  this  region,  stable 
crack  growth  is  predicted.  For  K^.'s  less  than  K_  (where  K_  is  analagous  to 
Kjo),  the  only  means  for  lowering  the  driving  force  is  a  reduction  in  crack 
length.  Hence,  Fuller  and  Thomson  proposed  crack  healing  to  be  a  viable 
mechanism  to  explain  the  fatigue  limit.  Cekirge,  Tyson,  and  Krausz  (4)  used  the 
concept  of  crack  healing  to  predict  a  fatigue  limit. 

The  objective  of  this  study  was  to  examine  the  crack  extension  behavior  of 
aged  Na20*3Si02  specimens  and  to  analyze  those  observations  based  upon  a  crack 
healing  mechanism. 

EXPERIMENTAL  PROCEDURE 

The  methods  of  specimen  preparation  and  crack-velocity  measurement  (using 
the  constant  moment  double  cantilever  beam  technique)  have  been  previously 

reported  (7).  Since  the  fatigue  limit  in  this  material  occurs  at  approximately 

1/2  1/2 

0.2  MPa>m  ,  the  specimens  were  aged  at  =  .142  MPa«m  '  for  two  hours  in 

deionized  water.  The  specimen  was  then  reloaded  to  stress  intensities  which 

produced  an  observable  crack  velocity.  The  crack  extension  after  aging  was 

recorded  as  a  function  of  and  time. 
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RESULTS  AND  DISCUSSION 

The  crack  extension  behavior  of  the  aged  NajOOSiOj  glass  specimens  is 

given  in  Figure  1.  After  aging,  time  delays  to  crack  elongation  of  0  to  1700 

seconds  were  observed.  The  corresponding  stress  intensities  were  between  .43 

1/2 

and  .38  MPa*m  ,  respectively. 

The  possibility  that  crack  healing  had  caused  the  observed  time  delays  - 
rather  than  crack  blunting  -  was  considered.  In  the  case  of  crack  healing,  no 
geometrical  alteration  of  the  crack  tip  is  expected.  Consequently,  one  would 
expect  the  cracks  to  extend  at  the  same  velocity  before  and  after  aging.  The 
magnitude  of  that  velocity  would  be  determined  by  the  stress  intensity  applied 
to  the  crack.  Figure  2  shows  the  velocity  of  the  crack  after  reloading  plotted 
against  the  velocity  before  aging.  The  1:1  correlation  line  is  also  plotted  to 
show  that  the  cracks  did  extend  at  the  expected  velocity.  Another  indication 
that  crack  blunting  was  not  responsible  for  the  effect  was  the  comparable  delay 
times  observed  for  a  series  of  sodium-aluminosilicate  glasses  (7),  whose 
corrosion  resistance  (and  hence  resistance  to  blunting)  should  greatly  exceed 
the  Na20*3Si02  glass. 

The  mechanism  by  which  crack  healing  could  produce  the  observed  time  delays 
was  envisioned  as  follows.  During  the  aging  process,  the  crack  heals  some  small 
amount.  Since  the  aqueous  environment  chemically  reacts  with  the  crack 
surfaces,  it  is  only  reasonable  to  anticipate  that  the  aqueous  species  will  be 
incorporated  within  the  healed  region.  Although  the  nature  of  this  healed 
region  is  unknown,  this  model  considers  the  healed  crack-tip  to  be  a  region  of 
high  water  content.  It  has  been  reported  that  glasses  of  high  water  content 
(12  vt%)  can  deform  plastically  (8).  It  has  also  been  shown  that  the  viscosity 
of  a  glass  is  lowered  by  the  incorporation  of  water  (9).  Assuming  that  the 
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Figure  1.  The  crack  elongation  behavior  of  Na90-3Si09 
after  aging  under  load  in  deionized  ‘‘water . 
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healed  material  is  deformable,  and  thus  can  exhibit  some  finite  viscosity,  a 
permanent  deformation  or  'strain'  can  account  for  the  time  delay. 

The  strain  rate  can  be  evaluated  if  we  assume  the  material  deforms  in  a 
viscous  nature 
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where  de/dt  =  strain  rate,  and 
r\  *  viscosity. 

The  process  which  is  expected  to  occur  at  the  crack-tip  if  a  deformable  material 
is  present  is  shown  in  Figure  3.  It  is  assumed  that  the  healed  region  begins  to 
deform  immediately  after  reloading.  This  deformation  continued  until  enough 
strain  is  developed  to  stretch  the  bonds  in  the  virgin  materials.  Once  the 
strain  in  these  bonds  reaches  some  limit,  the  crack  begins  to  re-propagate .  The 
experimetal  variables  which  can  be  controlled  are  aging  stress  intensity,  aging 
time,  and  aging  environment.  The  aging  stress  intensity  determines  the  driving 
force  for  crack  extension  during  the  aging  process.  According  to  the  Fuller  and 
Thomson  lattice  fracture  model  (3),  less  healing  would  be  expected  with  an 
increase  in  the  aging  stress  intensity.  In  these  experiments,  the  aging  stress 
intensity,  aging  time  and  aging  environment  were  held  constant.  Thus,  an 
invariant  amount  of  healing  would  be  expected,  and  thereby,  the  critical  amount 
of  strain  in  the  healed  region  would  be  constant.  However,  the  strain  rate  is 
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Figure  3.  Schematic  representation  of  a  crack  tip  which  has  healed 
to  form  a  deformable,  high-water  content  material  (the 
blackened  areas),  and  its  evolution  upon  reloading;  the 
times,  t,  correspond  to  'time  after  reloading.' 


predicted  to  be  dependent  upon  the  'reloading'  stress  intensity.  In  this 
manner,  the  delay  time  would  vary  with  reloading  stress  intensity. 
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In  order  to  quantify  this  effect.  Equation  1  can  be  used  by  assuming  a  crack 
orientation  such  as  shown  in  Figure  3: 


e2  "  de2/dt  =  t ox/ (-6n) ]  +  [a2/3n]  +  [o3/-6n] 


(5) 


Principal  stresses  can  be  expressed  as  a  simple  multiple  of  such  that: 


£ 


°l  ’  V2 


(6) 


Z 


where  is  a  scaling  factor.  Substituting  this  into  Equation  5  yields: 


a 


de2/dt  -  c^[(l/3n)  -  (Qj/bn)  -  (Q3/6n)) 


(7) 


The  stresses  can  be  expressed  as  a  function  of  the  applied  stress  intensity 
such  that: 


o.  =  A.K. 
l  i  1 


(8) 


Combining  Equations  7  and  8  yields  the  desired  result: 


dc2/dt  =  Kr(A[(l/3n)  -  (Qj/6ri)  -  ( Q3/6n)  ] ) 


(9) 
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B  is  Che  accumulation  of  parameters  relating  to  viscosity,  crack  orientation, 
and  crack  length.  These  parameters  will  be  constant  for  any  particular  material 
and  crack  system. 

This  model  shows  that  some  critical  amount  of  strain  (e  ^  )  is  necessary 

for  crack  repropagation.  This  strain  is  attained  at  some  critical  time  (t  .  ) 

crit 

which  depends  upon  the  reloading  stress  intensity  (assuming  the  other  parameters 
such  as  glass  composition,  aging  time,  aging  stress  intensity  and  environment 
are  constant).  This  results  in  the  following  boundary  conditions: 


at  t 
at  t 


0,  e2  =  0 

tcrit’  e2  G crit 


Using  these  conditions  to  evaluate  Equation  10  yields: 


crit 
d  e„ 


B  *K, 
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dt 
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Since  (K  )  ^  ^  *s  a  constant  for  each  individual  aging  experiment,  Equation 

11  can  be  integrated  to  yield: 


(KI  reload ' 


[B/e 


crit 


1  t 


crit 


(12) 


Figure  4  shows  the  aging  data  for  the  Na20*3Si02  glass  plotted  as 

predicted  in  Equation  12.  The  data  reflects  the  expected  trends  reasonably 

well.  A  linear  relationship  existed  for  the  (K_)  .  ,  values  where  a  finite 

I  reload 

delay  time  could  be  observed. 

Figure  5  presents  the  data  reported  by  Michalske  (2),  but  plotted  according 
to  the  expression  derived  in  Equation  (12).  The  trends  are  analagous  to  the 
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Figure  5.  The  delay  time  to  crack  elongation  reported  by  Michalske  (2) 
for  soda-lime-silica  glass;  here,  it  is  plotted  according 
to  equation  12. 


ones  obtained  in  this  study.  It  should  be  mentioned  that  the  slopes  (B/e 
of  the  linear  portions  in  Michalske's  data  and  the  data  obtained  in  this  study 
are  different.  This  is  to  be  expected  since  experimental  variables  (aging  time 
and  aging  stress  intensity)  and  composition  (viscosity  and  restoring  force)  were 
not  the  same  for  the  two  studies. 

CONCLUSIONS 

This  paper  suggests  that  alternatives  are  available  to  explain  the  sort  of 
crack  growth  behavior  previously  attributed  to  crack  tip  blunting.  The  aging 
experiments  presented  here  are  consistent  with  a  crack  healing  mechanism  where 
the  healed  material  is  modified.  This  healed  material  is  more  resistant  to 
crack-propagat ion  upon  reloading  because  it  can  undergo  a  plastic  deformation. 
The  presence  of  crack  healing  during  aging  supports  models  of  a  static  fatigue 
limit  resulting  from  crack  healing.  It  is  suggested  that  the  static  fatigue 
limits  observed  in  crack  extension  studies  can  be  attributed  to  comparable 
healing  mechanisms  thereby  eliminating  the  necessity  of  considering  crack  tip 
geometry  changes. 
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STRESS  CORROSION  OF  SODIUM- ALUMINOSILICATE  GLASSES:  A  COMPARISON  OF 

SMALL  AND  LARGE  CRACKS 

D.  Coon,+  J.  J.  Mecholsky  and  C.  G.  Pantano 
Department  of  Materials  Science  and  Engineering 
The  Pennsylvania  State  University 

Abstract 

The  double  cantilever  beam  and  dynamic  fatigue  techniques  were  used  to 
measure  slow  crack  growth  effects  in  a  series  of  sodium-aluminosilicate  glasses 
Na20*x  AljO^’O-x)  SiOj  with  x  varying  between  0  and  1.  The  constant  moment 
double  cantilever  beam  (CM-DCB)  technique  with  a  large  crack  measures  the 
threshold  Region  0,  and  Regions  I,  II  and  III  separately  while  the  dynamic 
fatigue,  or  stressing  rate,  technique  integrates  the  effect  of  Regions  0,  I,  II 
and  III  using  small  indentation  cracks.  Thus,  although  there  exists,  in 
principle,  a  direct  correspondence  between  slow-crack  growth  measurements  and 
fatigue  tests,  the  measured  n  values  are  not  always  the  same.  The  present  study 
shows  that  dynamic  fatigue  and  CM-DCB  analyses  yield  different  n-values;  and, 
there  were  larger  errors  in  n  values  in  the  dynamic  fatigue  analysis.  The  pH 
had  little  effect  on  the  measured  n  values  in  either  case,  but  surprisingly,  the 
strength  in  the  dynamic  fatigue  tests  increased  at  low  stressing  rates  in  a  pH 
11  solution.  This  increase  was  shown  to  be  due  to  a  relief  of  residual  stress 
and  not  crack  blunting,  and  to  be  consistent  with  the  low  velocity  plateaus 
observed  in  pH  11  solutions  with  the  CM-DCB  measurements.  Thus,  the  absolute 
differences  between  n-values  are  due  to  multiregion  crack  growth. 
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INTRODUCTION 


Several  investigators  have  attempted  to  experimentally  show  a  direct 
correlation  between  the  stress  corrosion  susceptibility  parameter,  n,  as 
determined  by  techniques  that  involve  large  cracks  such  as  the  constant  moment 
double  cantilever  beam  (CM-DCB)  technique  and  those  that  involve  small  cracks 
such  as  static  fatigue  (delayed  failure)  or  dynamic  fatigue  (stressing  rate)  for 
soda-lime-silica  glass.  In  the  case  of  static  fatigue,  Fuller^  showed  a 
one-to-one  mathematical  correspondence  between  the  velocity-stress  intensity 
relationship  and  the  slope  of  the  time  to  failure-applied  stress  curve. 

However,  the  mathematical  correspondence  between  dynamic  fatigue  and  the  V-K 
relationship  is  much  more  difficult  to  evaluate. 

Wiederhorn  et  al.®  suggested  that  differences  between  large  crack 
techniques  and  small  crack  techniques  could  result  from  differences  in  (1) 
residual  stresses  at  the  origin  of  fracture,  (2)  crack-tip  blunting,  and  (3) 
geometry  differences  between  large  and  small  cracks.  Residual  stresses  would  be 
expected  to  affect  the  measurements  to  a  greater  degree  on  small  cracks  than 
large  cracks.  If  crack  blunting  exists,  i.e.  crack  tip  rounding  due  to 
corrosion  processes,  then  it  would  be  expected  to  affect  large  cracks  more  than 
small  cracks.  Finally,  the  geometry  factor  in  crack  propagation  would  be 
different  for  a  straight  crack  as  in  CM-DCB  than  for  the  small  flaws  or 
semi -el li pt ical  indentation  cracks  in  fatigue  tests. 

The  present  study  uses  CM-DCB  and  dynamic  fatigue  to  examine  the  effect  of 
small  and  large  cracks  on  crack  propagation  in  a  series  of  Na20*x  Al^O'O-x) 

SiC^  glasses  with  x  =  0,  .2,  .4,  .6,  .8  and  1.0  as  a  function  of  environment. 
Thus,  it  is  possible  to  compare  crack  propagation  in  glasses  that  are  quite 
susceptible  to  corrosion  (x=0)  with  glasses  that  are  not  very  susceptible  to 
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corrosion  (x  ^  0*6).  In  order  to  farther  investigate  the  effect  of  varying 
corrosion  behavior,  the  crack  propagation  is  evaluated  in  solutions  of  pH“5  and 
pH=*ll . 

It  is  shown  that  different  n  values  are  obtained  for  the  CM-DCB  and 
stressing  rate  techniques  even  after  accounting  for  residual  stress  effects. 

The  absolute  values  of  n  determined  by  the  dynamic  fatigue  were  generally  lower 
than  those  determined  by  the  CM-DCB  technique.  Crack  blunting  is  eliminated  as 
a  reason  for  the  difference.  However,  the  strength  behavior  for  x=0  glass  at 
slow  loading  rates  was  shown  to  be  completely  consistent  with  the  CM-DCB  V-K^. 
results  at  low  velocity-low  stress  intensity  region.  Thus,  although  the 
mathematical  formulism  to  account  for  multiregion  crack-growth  does  not  yet 
exist,  it  is  possible  -  in  principle  -  to  translate  V-K  results  to  dynamic 
fatigue,  and  vice-versa. 


BACKGROUND 


The  relationship  between  crack  velocity  V  and  (Mode  I)  stress  intensity, 
is  often  represented  by  a  power  law  function:^’^ 

V  =  A(KI)n  (1) 


Although  there  is  no  fundamental  reason  for  using  this  relationship,  it  is 

mathematically  convenient  and  is  considered  a  valid  indication  of  stress 

*  , 

corrosion  behavior.  The  slope,  n,  of  the  log  v-log  plot  is  considered  a 


I 

i 

as 

* 


♦Although  the  slope  of  the  log  crack  velocity-stress  intensity  has  more 
fundamental  meaning,  i.e.,  it  is  related  to  the  activation  energy  and 
activation  volume,  it  is  more  common  for  the  slope  of  the  log 
velocity-log  stress  intensity  to  be  represented  as  a  measure  of  the 
resistance  to  slow  crack  growth. 
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characteristic  of  stress  corrosion  and  a  measure  of  the  resistance  to  slow  crack 
growth: 

log  V  =  log  A  +  n [ log  Kj]  (2) 

where:  V  =  crack  velocity 

A  =  constant 

Kj  =  applied  stress  intensity 

Direct  measurement  of  crack  velocity  at  known  applied  stress  intensities  permits 
the  determination  of  n.  However,  difficulties  in  observing  a  moving  crack  tip 
generally  limit  the  usefulness  of  this  'direct  measurement'  method  to  room 
temperature  studies. 

Evans^  has  proposed  a  method  by  which  n  can  be  obtained  from  strength 
measurements.  Assuming  the  rate  of  crack  extension  is  described  by  Equation  1, 
the  following  relationship  can  be  derived: 

log  Oj  =  log  A  +  [  1/n+l  ]  logo  (3) 

where  =  fracture  stress 

A  *  constant 
o  =  stressing  rate 

By  experimentally  determining  the  variation  of  fracture  stress  with  stressing 
rate,  n  is  easily  obtained.  This  is  known  as  the  dynamic  fatigue  technique. 

Such  an  'indirect  method'  eliminates  the  need  for  direct  crack  velocity 
measurement.  The  dynamic  fatigue  method  assumes  that  the  applied  stress  is  the 
only  stress  acting  on  the  strength  controlling  flaw.  Typically,  the  strength 
controlling  flaw  is  introduced  by  indentation.  Since  it  is  well  known  that 
residual  stresses  remain  after  the  indentation  process,  these  stresses  must  be 
removed  in  order  to  apply  the  Evan's  analysis.^ 


Recently  Fuller,  Lawn  and  Cook  have  reported  a  mathematical  correction 
for  the  residual  stresses  which  can  be  introduced  after  the  test.  They  derived 
the  following  expression: 

log  <j£  =  log  B  +  [l/n'  +  l]  log  a  (4) 

where  B  =  constant 

n'  =  apparent  n 

The  apparent  n'  (i.e.,  the  slope  of  the  dynamic  fatigue  curve)  is  related  to  the 
"true"  n  by: 

n  =  4/3  n'  -  2/3  (5) 

8 

The  analysis  of  Fuller,  Lawn,  and  Cook  has  refined  the  original  concepts  of 
Evans  7  such  that  removal  of  residual  stress  prior  to  testing  is  unnecessary. 


EXPERIMENTAL  PROCEDURE 

Glasses  in  the  family: 

Na20*xAl20^»(3-x)Si02  (*>) 

were  prepared  with  values  of  x  =  0.0,  0.2,  0.4,  0.6,  0.8  and  1.0.  In  this 
series  the  molar  concentration  of  Na20  was  held  constant  while  A120.j  was 
substituted  for  an  equimolar  amount  of  SiO 2>  These  glasses  are  herein 
identified  by  their  A^O^/Na^  molar  ratios  or  x-values.  In  this  manner, 
Na20-3Si02  was  labeled  0.0  glass  and  Na20*  A1  2  Si02  was  labeled  1.0  glass. 

The  specimen  preparation  and  testing  procedures  for  the  constant  moment 

9 

double  cantilever  beam  analysis  have  been  previously  reported  .  The  glasses  for 
dynamic  fatigue  analysis  were  cut  into  rectangular  bars  approximately  2.5  cm  x 
0.25  cm  x  0.25  cm.  The  bars  were  ground  with  600  grit  SiC  paper  such  that 
opposite  sides  were  parallel,  and  adjoining  sides  perpendicular .  One  face  was 


polished  with  15  urn  diamond  paste.  All  specimens  were  annealed  after  polishing 


The  dynamic  fatigue  specimens  were  indented  in  the  center  of  the  polished 

face  using  a  Vicker's  indenter.  The  indentation  load  was  constant  (8.73  N)  for 

all  compositions  investigated.  A  surface  flaw  approximately  100  pm  in  depth 

a 

resulted.  Consistent  with  the  analysis  of  Fuller,  Lawn  and  Cook  ,  no  attempt  to 
remove  the  residual  stress  around  the  indentation  was  attempted. 

Five  specimens  of  each  composition  were  broken  immediately  after 

indentation  in  four  point  bending  at  constant  crosshead  speeds  of  1.61  x  10™° 

—7  —7  —8 

m/s,  7.2  x  10  m/s ,  1.69  x  10  m/s,  and  4.23  x  10  m/s  on  an  Instron  tensile 

testing  machine.  The  polished  face  was  the  tensile  surface  in  each  case.  All 

compositions  were  tested  in  deionized  water  (pH  5).  In  addition,  specimens  of 

0.0,  0.6  and  1.0  glass  were  tested  in  NaOH  solution  (pH  11).  The  size  of  the 

fracture  initiating  cracks  were  obtained  by  direct  measurement  with  an  optical 

microscope. 


RESULTS 

The  results  of  the  constant  moment  double  cantilever  beam  (CM-DCB)  analysis 

q 

have  been  reported  .  The  reader  is  referred  to  that  study  for  the  experimental 
stress  intensity  -  crack  velocity  data  which  will  not  be  reproduced  here. 

Figure  1  summarizes  the  slopes  of  the  log  v-log  plots  determined  from  that 
CM-DCB  data. 

The  results  of  dynamic  fatigue  experiments  on  the  0.0,  0.2,  0.4,  0.6,  0.8, 

and  1.0  glasses  are  shown  in  Figures  2,  3,  4,  5,  6  and  7  respectively.  The 

n-values  calculated  from  these  data  are  also  given  in  Figure  1.  Figure  8  shows 

the  variation  in  the  critical  flaw  size  for  the  0.0  glass  tested  in  water  and 

NaOH  solution  as  a  function  of  stressing  rate.  Table  I  lists  the  fracture 

stresses  and  critical  flaw  sizes  for  the  0.0  glass  tested  in  water  and  NaOH 

-8 

solutions  at  a  crosshead  speed  of  4,23  x  10  m/s. 
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Figure  2a,  Dynamic  fatigue  plot  for  Na  0‘3Si09  in  deionized  water 
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Figure  2b.  Dynamic  fatigue  plot  for  Na.0*3Si0»  in  pH  11  NaOH  solution 
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Figure  3.  Dynamic  fatigue  data  for  Na^O’ 2 . SSiO^  in 
deionized  water. 
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Figure  4 


Dynamic  fatigue  data  for  Na.O.AAUO  .2.6SiO 
in  deionized  water.  J  2 


Figure  5.  Dynamic  fatigue  data  for  Na^O •  '  2  .  ASiO^ 

in  deionized  water. 


Figure  7.  Dynamic  fatigue  data  for  Na^’ A120  •  2Si02 
in  deionized  water.  '  4 


DISCUSSION 


The  slope,  n,  as  determined  by  the  CM-DCB  technique  ranged  from  16  to  23 

which  is  typical  of  silicate  glasses.^®  Figure  1  shows  that  n  increased 

linearly  with  Al^O^  concentration  for  0.0,  0.2,  and  0.4  glasses.  However,  n  for 

0.6,  0.8  and  1.0  glasses  was  invariant  and  equivalent  to  that  of  the  0.4  glass. 

These  trends  were  determined  by  'direct1  observation  of  crack  tip  extension  for 

large  cracks  and,  was  reported  to  be  readily  explained  by  the  elastic 
o 

properties . 

8 

The  Fuller,  Lawn  and  Cook  method  of  dynamic  fatigue  resulted  in  the 
variation  of  n  with  composition  as  shown  in  Figure  1.  The  n-value  increased  in 
a  non-linear  manner  with  alumina  concentration  over  a  range  of  10  to  27. 

Clearly,  the  n-values  obtained  by  the  two  techniques  are  not  identical  over  the 
compositional  range. 

It  is  also  noteworthy  that  the  error  observed  in  the  dynamic  fatigue 
analysis  is  approximately  300%  greater  than  that  observed  for  the  CM-DCB 
analysis.  This  results  from  the  manner  in  which  the  dynamic  fatigue  data  must 
be  statistically  analyzed.  In  dynamic  fatigue  analysis: 


slope 

The  error  in  n'  is  defined  by: 

error  (n')  =  error(slope) 

(slope) 

Since  the  n'  value  is  typically  much  greater  than  one, 
than  one.  Equation  8  results  in  a  large  error  in  n'. 
in  the  "true"  n  value  is  defined  by: 

error(n)  =  1 . 33 [error (n ') ] 


(7) 

(8) 

the  slope  is  much  less 
Subsequently,  the  error 

(9) 


Thus  the  dynamic  fatigue  technique  utilizing  indentation  analysis  is  inherently 
less  sensitive  to  measuring  n  than  the  CM-DCB  technique.  Nevertheless,  the 
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difference  in  n-values  between  the  two  measurement  techniques  are  outside  the 
range  of  experimental  error. 

Multi-Region  Crack  Growth 

In  addition  to  the  larger  error  in  the  indentation-dynamic  fatigue 
analysis,  there  is  an  inherent  difference  between  dynamic  fatigue  data  and  CM- 
DCB  data.  Dynamic  fatigue  data  is  taken  over  a  range  of  stress  intensity  levels 
from  0  up  to  a  critical  value,  K^.  Thus  the  crack  spends  time  in  the  various 

regimes  for  different  periods  of  time.  Depending  on  the  initial  crack  size 
the  jump  from  zero  to  the  initial  stress  intensity  is  abrupt  and  continues  at  a 
different  rate  depending  on  the  stressing  rate  and  the  initial  crack  size.  Thus 
the  final  failure  condition  is  an  integrated  effect  over  almost  the  entire  V-K^ 
curve  including  Regions  I,  II  and  III  and  0  (threshold).  The  process  is  more 
complicated  with  the  residual  stress  present  due  to  the  indentation,  but  in 
principle  the  same  conclusion  holds.  The  slope  one  measures  in  dynamic  fatigue, 
then,  is  really  a  combination  of  Regions  0,  I,  II,  and  III.  However,  the  slope 
reported  for  the  V-K^  data  obtained  with  the  CM-DCB  technique  is  only  Region  I! 

Fuller"*  showed  that  for  static  fatigue  experiments,  the  individual  points 
in  the  time-to-fai lure  plots  could  be  related  to  specific  points  on  the  V-K^ 
curve  (many  of  which  were  actually  in  the  threshold  region).  However,  the  sum 
total  of  the  points  on  the  static  fatigue  curve  did  not  result  in  a  slope 
corresponding  to  the  n  value  of  Region  I  in  the  V-K  data.  A  comparable  effect 
is  expected  in  dynamic  fatigue  data,  but  with  the  important  difference  that  the 
stress  is  varying  along  with  the  stress  intensity.  This  is  a  complication  that 
renders  a  detailed  mathematical  analysis  difficult,  but  the  phenomenology  is 
essentially  identical.  In  general,  then,  one  cannot  expect  to  derive  the  same  n 
value  from  the  two  types  of  experiments  in  glasses  which  exhibit  multiregion 
crack  growth. 


-  -”--  W"-  A."-.  » 
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Residual  Stress  Effects 


In  general,  the  n-values  were  insensitive  to  pH,  but  the  0.0  glass  tested 

,-3 


in  10  M  NaOH  solution  exhibited  anomalous  behavior  at  low  stressing  rates 
(Figure  2b).  The  fracture  stress  increased  with  decreasing  stressing  rate  as 
opposed  to  the  decrease  observed  when  this  composition  was  tested  in  water. 

This  type  of  upturn  in  fracture  stress  at  low  stressing  rates  has  been  observed 
by  Yavuz,  Tressler,  and  Minford^  for  SiC  tested  at  high  temperatures.  They 

hypothesized  the  behavior  to  result  from  crack  tip  blunting. 

.12 


Other  work  with  these  glasses  has  shown  that  the  time  to  repropagate  the 
crack  on  a  CM-DCB  specimen  after  aging  under  water  is  approximately  the  same  for 
the  x  =  0,  x  =  0.6  and  x  =  1.0  glasses.  Since  these  glasses  have  enormously 
different  corrosion  behaviors,  these  results  suggest  that  crack  blunting  during 
aging  cannot  explain  the  delay  time  before  crack  propagation;  i.e.,  one  would 
expect  different  amounts  of  blunting  due  to  the  large  differences  in  corrosion 
behavior  for  these  glasses.  Thus,  an  alternate  explanation  is  required  for  the 
upturn  in  strength  at  low  stressing  rates.  This  alternate  explanation  can  be 
the  relief  of  residual  stress. 

Anstis  and  co-workers  have  shown  that  an  indentation  crack  system 
consists  of  two  crack  systems:  radial  and  lateral.  Extension  of  the  lateral 
cracks  was  shown  to  be  responsible  for  residual  stress  relaxation  whereas 
extension  of  radial  cracks  was  shown  to  be  responsible  for  the  strength  decrease 


in  bend  bars.  Since  the  total  stress  intensity,  K  ,  experienced  by 

cot 

specimen  consists  of  two  components: 


KI  =  Ki  +  KI 

tot  applied  residual 


(10) 
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the  critical  condition  is  met  when  K  reaches  the  fracture  toughness,  K_  . 

I .  „  Ic 

tot 

If  K  is  reduced  due  to  lateral  crack  extension  prior  to  failure, 

res idual 

must  necessarily  increase.  In  this  manner,  the  applied  stress  at 

applied 

failure  is  larger  if  the  residual  stress  is  relaxed  before  failure. 

A  more  distinct  low  velocity  plateau  was  observed  for  the  0.0  glass  tested 
at  pH  11  than  when  tested  in  water  (Figure  2b).  Moreover,  the  threshold 
decreased  to  a  lower  value  in  pH  11  solution  (NaOH).  The  change  in  residual 
stress  intensity  is  determined  by  the  change  in  lateral  crack  size: 


K 


I 


res idual 


(11) 


where:  xD  =  constant  =  0.016 
P  =  indentation  load 
C.  =  initial  crack  size 

l 

Cj  =  final  crack  size 

The  extension  of  the  lateral  cracks  would  be  expected  to  progress  at  a  much 

faster  rate  in  a  pH  11  solution  relative  to  water  because  dissolution,  and  hence 

corrosion,  occurs  much  more  rapidly  in  this  solution.  Also,  lateral  crack 

extension  is  expected  at  lower  levels  of  .  Therefore,  greater 

residual 

reductions  in  residual  stress,  and  larger  lateral  cracks  at  failure,  are 

anticipated  in  the  0.0  glass  tested  in  pH  11  relative  to  water.  This  is  a  rate 

dependent  phenomenon  since  a  certain  amount  of  time  is  required  for  crack 

extension,  and  hence,  residual  stress  relaxation.  Thus,  it  is  reasonable  that 

it  is  only  observed  at  the  low  stressing  rates. 

The  applied  stress  intensity  is  determined  by  both  the  crack  length,  c,  and 

applied  stress,  a  : 

app 


Kj  =  Yaapp"\/c"” 


(12) 


where  Y  =  geometric  and  location  constant. 

It  is  apparent  that  for  a  given  level  of  applied  stress,  a  larger  (or  longer) 
crack  results  in  a  larger  stress  intensity.  It  is  hypothesized  that  the  long 
lateral  cracks  at  pH  11  grow  along  with  the  radial  cracks  during  strength 
testing.  In  contrast,  the  'short'  lateral  cracks  obtained  in  water  are  stable 
while  the  radial  cracks  grow.  Thus,  we  predict  higher  fracture  stresses  and 
lower  critical  flaw  sizes  in  pH  11.  The  higher  fracture  stresses  have  already 
been  demonstrated.  The  critical  flaw  sizes  are  shown  in  Figure  8.  Although  the 
flaw  sizes  are  comparable  for  the  three  higher  stressing  rates,  the  flaw  sizes 
are  indeed  smaller  at  the  low  stressing  rate  for  pH  11  compared  to  water.  This 
supports  the  concept  that  the  observed  upturn  in  fracture  stress  in  pH  11  is  due 
solely  to  residual  indentation  stress  effects  and  is  not  caused  by  blunting. 

Additional  evidence  that  residual  stress  was  responsible  for  the  upturn  was 
obtained  by  removing  the  residual  stress  prior  to  testing  via  a  post-indentat  '.on 
anneal  cycle.  Table  I  shows  that  both  the  fracture  stress  and  flaw  size  were 
comparable  for  pH  5  and  11  when  tested  after  annealing.  Wiederhorn  et  al.^ 
argue  that  it  is  unlikely  that  the  behavior  of  large  cracks  differ  significantly 
from  small  cracks  when  only  Region  I  or  Region  III  crack  growth  is  involved  in 
degrading  the  strength.  Recall  that  Region  1  is  environmentally  dependent  and 
Region  III  is  environmentally  independent  crack  growth.  Region  II,  however,  is 
the  combined  effect  of  Region  I  and  Region  III  crack  growth  behavior,  and  will 
influence  crack  growth  differently  for  small  and  large  cracks.  In  Regon  II,  the 
portion  of  the  crack  away  from  the  surface  moves  slower  than  the  portion  nearer 
the  surface  because  of  water  depletion.  Thus,  the  shape  of  cracks  are  changed 
due  to  the  water  content  near  the  crack  boundary.  Because  of  the  difference  in 
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Che  crack  shapes  betweea  large  (flat)  and  small  (elliptical)  cracks  in  the  two 
tests,  respectively,  physical  constraints  on  the  crack  front  would  be  expected 
to  be  different  for  the  two  types  of  cracks;  i.e.,  crack  growth  in  Region  II  is 
not  expected  to  be  the  same  for  the  two  type  of  cracks.  Although  we  cannot 
demonstrate  unequivocally  the  validity  of  our  hypothesis,  it  seems  reasonable 
that  the  two  cracks  will  behave  differently.  Elliptical  cracks  can  vary  by  as 
much  as  50 %  in  along  the  perimeter  while  at  flat  cracks  is  expected  to  be 
nearly  constant  except  near  the  sides. 


CONCLUSIONS 

The  absolute  values  of  n  determined  by  dynamic  fatigue  were  generally  lower 
than  those  determined  by  the  CM-DCB  technique.  This  is  because  the  CM-DCB  test 
measures  the  velocity  of  a  crack  at  a  given  stress  intensity,  whereas  the 
dynamic  fatigue  test  measures  velocities  along  various  parts  of  the  V-K^  curve 
depending  on  the  initial  crack  size  and  loading  rate.  The  dynamic  fatigue  n 
parameters  were  consistent  with  the  double  cantilever  beam  analyses  in  showing 
that  pH  had  little  effect  on  n,  while  the  observed  'upturn'  in  fracture  stress 
at  low  stressing  rate  observed  for  the  0.0  glass  was  related  to  the  low  velocity 
plateaus  observed  in  the  stress  intensity  -  crack  velocity  curves  through 
corrosion  and  residual  stress.  In  the  threshold  region  and  in  Region  II,  the 
geometry  difference  between  dynamic  fatigue  tests  with  semi-elliptical  small 
cracks  and  the  CM-DCB  test  with  large  straight  cracks  further  accounts  for  some 
of  the  apparent  difference  in  slow  crack  growth  behavior.  Finally,  it  was  shown 
that  dynamic  fatigue  analysis  produced  n  parameters  exhibiting  larger  errors 
when  compared  with  n  parameters  obtained  by  the  constant  moment  double 
cantilever  beam  technique.  This  difference  is  due  to  the  inherent  error  for  the 
indentation-dynamic  fatigue  technique. 


Thus,  the  major  differences  in  n  values  for  the  CM-DCB  and  dynamic  fatigue 
techniques  for  the  sodium-aluminosilicate  glasses  is  due  to  residual  stress  in 
the  small  cracks,  the  integrating  effect  of  multiregion  crack  growth  upon 
dynamic  fatigue  tests,  and  to  a  lesser  extent,  the  consequential  larger  inherent 
error  in  the  indentation-dynamic  fatigue  technique. 
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ABSTRACT 

The  surface  chemical  reactions  which  might  contribute  to  stress 

corrosion  in  f luorozirconate  glasses  have  been  investigated  and  the 

results  are  related  to  the  slow  crack  growth  behavior  of  these  glasses. 

The  surface  studies  indicate  that  aqueous  environments  are  especially 

aggressive.  It  is  found  that  the  penetration  and  bulk  diffusion  of 

molecular  water  precedes  any  hydrolysis  or  dissolution  of  the  glass 

structure.  The  chemisorption  of  water  to  produce  oxide  or  hydroxide 

species  is  not  prevalent  except  in  the  case  of  the  basic  aqueous 

environments.  The  most  notable  feature  of  the  stress  intensity-crack 

velocity  diagrams  is  a  distinct  and  reproducible  stress-intensity 

1/2  . 

threshold  at  K^  *  .15-. 20  MPa-ra  in  the  environments  which  contain 

water.  This  threshold  is  observed  for  both  gaseous  and  liquid 

environments  even  though  the  surface  chemistry  and  corrosion  in  water 

vapor  and  liquid  water  are  very  different.  The  wet  nitrogen  atmospheres 

exhibit  a  velocity  plateau  whose  magnitude  increases  systematically  with 

the  water  activity;  it  approaches  the  plateau  observed  in  pure  liquid 
-3 

water  at  ~10  m/s.  Altogether,  it  appears  that  the  crack  growth 
mechanism  may  be  independent  of  any  chemisorption  reactions  at  the  crack-tip 
Perhaps,  the  formation  of  a  coordinate  bond  between  molecular  water  and 
the  barium  or  zirconium  cations  weakens  the  coulombic  forces  of  attraction 
with  the  fluorine  anions,  and  thereby,  influences  the  energy  necessary  for 
crack  propagation.  This  interaction  also  appears  to  be  the  first  step  in 
the  corrosion  reaction.  It  is  further  suggested  that  the  stress-intensity 
thresholds  and  velocity  plateaus  may  be  related  to  stress-enhanced  water 
diffusion  ahead  of  the  crack  tip  where  it  modifies  both  the  material  and 
stress  distribution. 
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1 .  INTRODUCTION 


It  has  already  been  shown  chat  both  f iuorozirconate  and  barium-thorium 
fluoride  glasses  exhibit  a  susceptibility  to  dynamic  fatigue  and 
slow-crack  growth  (1).  These  data  indicate  that  water  is  an  especially 
active  environment,  but  the  stress  intensity-crack  velocity  (KV)  diagrams 
have  unique  characteristics  relative  to  those  observed  in  silicate 
glasses.  Otherwise,  very  little  has  been  reported  on  the  fracture  and 
stress  corrosion  behavior  of  these  glasses. 

In  a  more  general  sense,  the  stress  corrosion  mechanisms  in  ionic 
materials,  whether  glassy  or  crystalline,  have  not  been  clarified.  In  a 
study  of  slow  crack  growth  in  single-crystalline  alkaline  earth  fluorides 
(CaFj*  SrF2  and  BaF2^’  ^oth  H2°  an<*  were  found  to  enhance  the 
crack-growth  rate  (2).  The  delayed  failure  of  polycrystalline  MgF^  (3) 
and  ZnSe  (4)  were  also  reported,  and  here  too,  water-enhanced  slow  crack 
growth  was  observed.  Although  these  studies  described  the  importance  of 
dislocation  motion  and  microstructure  in  the  fracture  process,  they  did 
not  lead  to  the  establishment  of  a  stress  corrosion  mechanism. 

In  the  most  definitive  studies  of  stress  corrosion  mechanisms  in 
silicate  glasses  (5-7),  the  availability  of  an  extensive  data  base 
concerning  their  surface  chemistry  (adsorption,  hydration/dehydration, 
chemical  reactivity,  etc.)  and  corrosion  behaviors  was  of  great  benefit. 

A  comparable  data  base  for  fluoride  glasses  does  not  exist.  A  number  of 
papers  have  appeared  recently  concerning  the  corrosion  of  fluoride  glasses 
(8-15),  but  the  reaction  mechanisms  are  still  not  well  understood.  The 
adsorption  and  other  surface  chemical  reactions  have  yet  to  be  reported 
for  these  glasses,  although  some  adsorption  studies  of  water  vapor  on 
crystalline  alkaline-earth  fluorides  have  been  reported  (16-21). 

There  is  no  question  that  these  fluoride  glasses  are  susceptible  to 
stress  corrosion,  and  this  will  certainly  limit  many  of  their  potential 
applications  for  infrared  fiber  optics.  In  this  paper  we  report  further 
on  the  surface  chemistry  and  slow  crack  growth  behavior  of  f Iuorozirconate 
glasses.  The  study  of  chemical  effects  at  fracture  surfaces  is  emphasized 
where  methods  including  low-energy  ion  scattering  (ISS),  Auger  electron 
spectroscopy  (AES)  and  secondary  ion  mass  spectroscopy  (SIMS)  are  used  for 
characterization.  The  constant  moment  double  cantilever  beam  technique 
was  employed  to  determine  the  slow  crack  growth  behavior  in  a  wide  variety 
of  liquid  and  gaseous  environments  including  water,  water  vapor,  methanol, 


methanol  vapor,  ammonia,  formamide  and  water-methanol  mixtures.  These 
data  are  used  to  discuss  the  mechanism  of  stress  corrosion  in  fluoride 
glasses . 

2 .  EXPERIMENTAL  PROCEDURE 

All  of  the  double  cantilever  beam  specimens  were  cut  from  a  single 
billet  of  glass.  The  glass  -  whose  composition  is  ,57ZrF^ ‘.SbBaF^  * 
•03LaF^ *04A1F^  (ZBLA)  -  was  fabricated  by  Le  Verre  Fluore.  The  beams  - 
appproximately  6mm  x  30mm  x  3mm  -  were  cut  and  oriented  so  that  the  crack 
followed  a  path  of  constant  thermal  history  through  the  glass.  A  2mm  wide 
groove  -  approximately  2mm  deep  -  was  machined  down  the  center  of  the  beam 
to  guide  the  crack.  The  groove  and  faces  of  the  glass  plates  were 
polished  to  eliminate  residual  stresses  due  to  the  machining  operations. 
The  cracks  were  initiated  by  thermal  shock  with  a  hot  nichrome  wire.  The 
crack  was  then  extended  under  load  into  the  middle  two-thirds  of  the 
specimen  where  the  crack  velocities  were  recorded.  The  constant  moment 
double  cantilever  beam  technique  (22)  was  used  with  the  crack  elongation 
being  monitored  through  a  40X  traveling  microscope.  The  specimens  were 
enclosed  in  a  sealed  chamber  which  could  be  purged  with  a  gaseous 
environment  before  and/or  during  the  test.  The  controlled  humidity 
environments  were  produced  by  mixing  measured  amounts  of  dry  nitrogen  and 
water-saturated  nitrogen  gas.  The  humidity  was  measured  with  a  dew  point 
hygrometer . 

The  clean  surfaces  used  for  adsorption  studies  were  created  by 
fracturing  pre-notched  rods  of  the  glass  in  an  ultra-high  vacuum 
analytical  chamber.  The  chamber  could  be  baked  and  ion-pumped  to  a  base 
pressure  of  less  than  1  x  10  torr  before  the  fracture  to  eliminate 
water  vapor  and  other  residual  gases.  The  species  present  in  the  vacuum 
environment  were  monitored  with  a  quadrupole  residual  gas  analyzer.  The 
chamber  incorporated  a  low-energy  ion  scattering  spectrometer  (ISS)  and  an 
Auger  electron  spectrometer  (AES)  for  analysis  of  the  'clean'  surfaces 
before  and  after  controlled-amounts  of  water  vapor  were  introduced  through 
a  leak  valve.  This  system  was  also  used  for  the  chemical  analysis  of 
fracture  surfaces  created  externally  in  air  and  liquid  environments.  It 
is  important  to  point  out  that  a  careful  evaluation  of  the  electron  beam 
effects  which  occur  during  Auger  analysis  of  these  glasses  was  carried 
out;  this  evaluation  showed  that  the  electron  beam  current  density  must  be 
maintained  at  or  below  6  uA/cm  to  prevent  beam  damage  to  the  specimen. 
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The  glasses  were  hydrated  at  elevated  temperatures  in  a  tube  furnace 
under  flowing  Ar  or  02  gases  which  were  first  saturated  with  water  vapor 
at  30C  or  75C.  These  glasses,  as  well  as  those  hydrated  at  room 
temperature  in  water-methanol  liquid  mixtures,  were  depth  profiled  with 
secondary  ion  mass  spectroscopy  (STMS). 

3.  CRACK-GROWTH  STUDIES 

The  K-V  diagrams  measured  for  ZB LA  in  a  wide  variety  of  environments 
are  summarized  in  Figures  1  and  2.  In  Figure  1,  the  presence  of  a 
well-defined  threshold  -  followed  by  a  velocity  plateau  -  appears  to 
characterize  the  behavior  in  humid  gas  and  in  liquid  water.  The  velocity 
plateau  shifts  to  higher  values  with  increasing  water  activity  and  reaches 
its  maximum  value  in  liquid  water.  The  K-V  diagrams  in  Figure  2  show  that 
ammonia  vapor  behaves  very  much  like  a  humid  gas ,  whereas  the  methanol 
liquid,  methanol  vapor  and  formamide  extend  the  region  of  slow  crack 
growth  to  lower  stress-intensity.  The  absence  of  a  threshold  or  velocity 
plateau  for  the  nonaqueous  environments  indicates  that  they  are  each 
active  in  promoting  crack  growth. 


0.15  0.25  0.35 
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Figure  1.  The  stress  intensity  -  crack  velocitv  diagrams  for  a  fluoro- 
zirconate  glass  in  aqueous  environments;  the  crack  velocity 
in  water  varies  over  the  range  represented  by  the  box. 
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Figure  2.  The  stress  intensity  -  crack  velocity  diagrams  for  a  fluoro- 
zirconate  glass  in  aqueous  and  non-aqueous  environments. 

It  is  interesting  to  note  that  decreasing  the  activity  of  water  vapor 
in  the  inert  nitrogen  gas  phase  systematically  lowers  the  velocity 
plateau,  whereas  decreasing  the  activity  of  liquid  water  by  dilution  with 
methanol  yields  a  straight  line  whose  slope  depends  upon  the  water 
activity  (1).  The  slope  of  this  line  can  be  used  to  define  the  stress 
corrosion  susceptibility  (n)  which  is  “77  in  pure  methanol  and  is  “40  in 
methanol  containing  1-5%  water.  The  absence  of  thresholds  and  plateaus  in 
the  water-methanol  mixtures  suggests  that  these  features  may  be  related  to 
the  formation  of  a  surface  reaction  product  between  adsorbed  or  liquid 
water  and  the  glass.  The  reaction  products  can  be  neutralized  or  carried 
away  by  the  methanol  'solvent'  in  the  liquid  mixtures  but  not  in  the  humid 
gas.  In  pure  water,  these  reaction  products  probably  accumulate  in  the 
crack,  whereas  in  the  pure  methanol  vapor  and  formamide  they  may  not  form. 
Later,  it  is  suggested  that  HF  may  be  a  surface  reaction  product  which  is 
unique  to  the  environments  which  exhibit  thresholds  and  velocity  plateaus. 

It  is  quite  clear  that  water  is  an  aggressive  crack  growth  agent.  This 
is  consistent  with  the  aqueous  corrosion  behavior  of  f luorozirconate 


glasses  where  rapid  and  extensive  surface  hydration,  dissolution  and 
crystalline  reaction  product  formation  are  evident  (8-15).  This  corrosion 
reaction  is  greatly  accelerated  in  acidic  solutions,  but  is  practically 
non-existent  in  basic  solutions.  For  this  reason,  additional  crack  growth 
measurements  were  made  in  a  pH  11  ammonium  hydroxide  solution.  The  K-V 
diagram  in  this  solution  was  comparable  to  the  curve  for  liquid  water. 
Nonetheless,  these  fracture  surfaces  were  free  of  any  visible  degradation 
whereas  in  the  pure  water  environments  a  highly  corroded  surface  is 
typically  observed. 

Additional  crack  growth  measurements  were  made  in  mixtures  of  methanol 
and  either  acid  or  base.  The  pH  of  the  mixtures  was  systematically  varied 
between  1  and  12.  These  data  are  shown  in  Figure  3  and  they,  too, 
indicate  an  insensitivity  of  the  crack  velocity  to  pH.  The  data  points 
are  scattered  about  a  line  whose  slope  gives  an  n  value  of  about  40.  In 
view  of  the  obvious  pH  dependence  of  the  corrosion  and  dissolution 
reactions  (15),  these  observations  suggest  that  the  stress  corrosion 
mechanism  in  fluorozirconate  glasses  may  be  distinct  relative  to  the 
corrosion  mechanism. 


0.15  0.20  0.25  0.30  0.35 
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Figure  3.  The  sub-critical  crack  velocities  measured  in  mixtures  of 
methanol  and  either  water,  acid  or  base;  in  all  cases,  the 
water  content  is  less  than  1%. 
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4.  SURFACE  STUDIES 

The  reaction  o£  clean  fracture  surfaces  exposed  to  water  vapor,  as  well 
as  the  composition  of  fracture  surfaces  created  in  various  environments, 
have  been  investigated  with  ion-scattering  spectrometry  (ISS)  Auger 
electron  spectroscopy  (AES)  and  secondary  ion  mass  spectroscopy  (SIMS). 

The  ISS  technique  is  sensitive  to  the  elemental  composition  of  the 
outermost  atomic  layer  of  a  solid.  It  has  been  used  here  to  determine 
whether  oxidation  and/or  the  chemisorption  of  water  are  prevalent 
reactions  in  the  surface  monolayer  of  clean  fracture  surfaces.  The  upper 
spectra  in  Figure  4  represents  a  clean  surface  created  by  fracturing  a  rod 
of  ZBLA  glass  under  ultra-high  vacuum  conditions.  The  peaks  due  to  F  and 
Al  are  clearly  evident,  whereas  the  peaks  due  to  Zr,  Ba  and  La  exhibit 
some  overlap.  The  fact  that  the  fluorine  signal  is  strong  -  in  spite  of 
its  lower  scattering  cross-section  relation  to  zirconium  and  barium  - 
indicates  that  the  fluorine  ions  shield  the  cations  at  the  surface,  and 
thereby,  dominant  the  sites  in  the  surface  monolayer.  Clean  fracture 
surfaces  of  this  type  have  been  exposed  in-situ  to  air  and  water  vapor, 
and  in  both  cases,  no  oxygen  signal  due  to  oxide,  hydroxide  or  water 
adsorption  was  detected.  In  contrast,  creation  of  the  fracture  surface  in 
air  leads  to  the  adsorption  of  some  oxygen  (or  hydroxyl)  as  shown  in  the 
second  spectra  of  Figure  4.  This  suggests  that  active  surface  adsorption 
sites  may  be  generated  by  the  stress  which  is  present  during  the  fracture 
event.  In  general,  though,  the  extent  of  oxygen  chemisorption  in  the 
outermost  monolayer  of  the  fracture  surfaces  is  minimal. 

A  similar  ISS  analysis  of  water  adsorption  was  reported  by  Strecker 
(21)  for  clean  surfaces  of  single-crystal  CaF^.  He  found  that  water 
adsorption  was  practically  non-existent  until  the  temperature  was  lowered 
to  the  point  where  an  epitaxial  condensation  of  water  was  observed.  His 
data  also  provided  an  indication  of  the  monolayer  sensitivity  of  this 
method;  i.e.,  condensation  of  a  single  monolayer  of  water  was  sufficient 
to  completely  extinguish  the  Ca  and  F  peaks  in  the  ISS  spectra,  and 
correspondingly,  to  generate  a  strong  oxygen  signal.  Thus,  the  oxygen  or 
water  coverage  represented  in  the  middle  spectra  of  Figure  4  is  small 
indeed . 

The  lower  spectra  in  Figure  4  represents  a  surface  created  by  fracture 
in  water.  Due  to  the  aggressive  reaction  between  f luorozirconate  glass 
and  water,  though,  one  must  recognize  that  the  surfaces  have  been  corroded 
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The  ISS  studies  have  shown  that  neither  the  chemisorption  of  water,  nor 
the  formation  of  oxides  and  hydroxides,  are  prevalent  reactions  at  the 
surface  of  ZBLA.  However,  the  heavy  Zr  and  Ba  cations  are  responsible  for 
the  high  background  in  the  low  energy  region  of  the  ISS  spectra,  and  for 
this  reason,  the  sensitivity  to  oxygen  is  less  than  ideal.  And  since  the 
low  energy  ion  scattering  occurs  almost  exclusively  in  the  outermost 
atomic  layer,  it  is  also  possible  that  ion  scattering  by  oxygen,  hydroxyl 
or  water  species  is  physically  shielded  at  the  surface  by  the  heavier 
cations.  Thus,  a  corresponding  set  of  surface  analytical  studies  were 
attempted  with  Auger  electron  spectroscopy  (AES).  The  AES  technique 
averages  the  elemental  composition  in  the  outermost  3. 0-5.0  nm  of  the 
surface;  moreover,  the  AES  method  is  exceedingly  sensitive  to  oxygen. 

The  Auger  spectra  in  Figure  5  show  that  oxygen  can  be  detected  by  AES 
even  on  the  fracture  surfaces  created  in  vacuum,  but  it  is  not  yet  known 
whether  this  oxygen  is  intrinsic  to  the  glass  or  was  adsorbed  from  the 
residual  vacuum  via  electron-beam  stimulated  adsorption  (23).  In  any 
event,  it  corresponds  to  a  relatively  low  concentration  (~5-6  atomic 
percent).  After  in-situ  exposure  to  “1800  Langmuirs  of  water  vapor,  the 
Auger  spectra  showed  no  substantial  increase  in  oxygen,  but  the  zirconium 
signal  was  eliminated  and  some  increase  in  the  fluorine  signal  was 
observed.  The  attenuation  of  the  zirconium  signal  could  be  due  to 
shielding  by  some  adsorbed  species,  but  it  is  unlikely  that  this  species 
is  water  because  it  would  have  produced  a  significant  increase  in  the 
oxygen  signal.  The  fact  that  hydrogen  is  undetectable  by  AES  suggests 
that  a  hydrogen  species  may  be  involved;  The  observed  increase  in 
fluorine  concentration  leads  one  to  consider  HF  surface  reaction  products. 
The  HF  could  form  due  to  a  reaction  between  ZBLA  and  H^O  with  the 
resulting  hydroxyl  specie  diffusing  further  into  the  glass. 

The  other  spectra  in  Figure  5  was  obtained  on  a  fracture  surface 
created  in  air.  One  notes  again  that  little  oxygen  is  present  on  the 
surface,  but  a  carbon  signal-presumably  an  adsorbed  hydrocarbon-is 
observed.  After  the  analysis  of  six  fracture  surfaces  created  in  air,  it 
became  apparent  that  the  oxygen  content  could  vary  between  5%  and  10%. 

The  C  and  0  ratio  was  not  constant,  and  so,  hydrocarbons,  rather  than  CO 
or  CO 2>  may  be  involved.  In  general,  though,  the  spectra  in  Figure  5 
further  support  the  contention  that  these  fracture  surfaces  are  not 
especially  reactive  in  humid  atmospheres.  The  oxygen  species  -  whatever 
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Figure  5.  Auger  spectra  for  fracture  surfaces  of  ZBLA;  the  center 

spectra  represents  the  clean  vacuum-fractured  surface  after 
exposure  to  1800  Langmuirs  of  water. 
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Figure  6.  Auger  spectra  for  surfaces  of  ZBLA  which  were  created  by 

fracture  in  liquid  environments;  the  spectra  in  Figures  5„ 
and  6  were  obtained  using  a  3KeV  electron  beam  at  'buA/cm^ 
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they  are  -  do  not  cover  the  surface,  but  rather,  are  present  at  a 
relatively  low  concentration  (at  least  relative  to  the  anion  content  of 
the  ZBLA  glass). 

The  Auger  spectra  in  Figure  6  represent  the  influence  of  some  liquid 
environments  upon  the  surface  composition.  In  all  cases,  the  glasses  were 
fractured  under  the  liquid,  and  then,  were  immediately  extracted  and  dried 
in  a  stream  of  dry  nitrogen.  The  specimen  fractured  under  water  shows  a 
significant  attenuation  of  the  barium  and  zirconium  signals,  some  carbon, 
and  background  levels  of  oxygen;  the  obvious  decrease  in  the  [Ba,Zr]-to- 
[F]  ratio,  again,  points  to  the  possible  formation  of  HF  reaction  products 
on  the  surface.  In  contrast,  the  glasses  fractured  in  methanol  and 
ammonium  hydroxide  show  substantial  increases  in  oxygen  content.  The 
oxygen  concentration  represented  in  these  two  spectra  is  approximately  20 
atomic  percent.  The  fact  that  nitrogen  is  not  detected  at  the  surface 
created  in  the  ammonium  hydroxide  indicates  that  specific  adsorption  of 
oxygen  or  hydroxyl  has  occurred.  The  presence  of  these  oxygen  levels  at 
the  surface  has  not  attenuated  the  fluorine  signal,  and  so,  the  exchange 
of  F  and  OH  seems  unlikely  even  in  the  highly  basic  hydroxide  solution.  A 
noticeable  increase  in  carbon  content  was  always  observed  on  fracture 
surfaces  created  in  the  ammonium  hyroxide  which  may  be  due  to  enhanced 
adsorption  of  hydrocarbon,  or  possibly  C0/C02  adsorption,  at  oxidized  or 
hydroxylated  surface  sites.  The  surfaces  created  under  methanol 
consistently  showed  increases  in  both  the  oxygen  and  carbon,  and  in 
general,  the  C  to  0  ratio  was  usually  larger  than  that  observed  for  all 
other  surface  conditions.  This  suggests  a  strong  physical  adsorption  of 
the  methanol  species,  especially  since  the  exposure  to  ultra-high  vacuum 
has  not  triggered  their  desorption.  In  the  case  of  alkali-chlorides, 
bromides  and  iodides,  the  physical  adsorption  of  methanol  has  been 
observed  and  attributed  to  hydrogen  bonding  with  the  surface  anions  (24). 

The  SIMS  technique  has  been  used  to  depth  profile  a  wide  variety  of 
ZBLA  glass  surfaces  after  exposure  to  humid  and  aqueous  environments.  In 
contrast  to  AES  and  SIMS,  this  method  is  relatively  insensitive  to  the 
outermost  atomic  layer(s)  of  the  solid,  but  is  ideally  suited  for 
profiling  .1  to  10.0pm  into  the  surface.  The  fact  that  these  glasses  do 
corrode  in  aqueous  environments,  but  show  little  reaction  at  the 
glass/solution  interface,  suggests  that  the  reactants  penetrate  the 
material  and  breakdown  the  structure  from  within.  Since  it  was  already 
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known  that  Che  hydration  of  ZB LA  is  very  slow  at  ambient  temperatures  - 
even  when  the  relative  humidity  is  as  high  as  95%  -  hydration  experiments 
were  carried  out  at  200C  and  275C  in  water-saturated  Ar  or  0£. 

Figure  7  shows  that  both  oxygen  and  hydrogen  are  found  to  penetrate 
nearly  5,000.0  nm  into  the  surface  within  4  hours  at  200C.  It  is 
especially  noteworthy  that  the  distribution  of  glass  constituents  has  not 
been  influenced  by  the  water  penetration.  In  multicomponent  silicate 
glasses,  hydration  in  water  vapor  usually  leads  to  a  redistribution  of 
modifier  species  and  the  growth  of  a  distinct  surface  layer.  The 
temperature,  time  and  carrier  gas  (Ar  vs  0^)  were  varied,  and  in  all 
cases,  only  the  apparent  penetration  of  H,  0  and  in  some  instances  C,  were 
revealed  in  the  profiles.  This  is  clearly  consistent  with  the  ISS  and  AES 
surface  studies  in  air  and  water  vapor  which  showed  only  minimal  surface 
reactions.  A  quantitative  analysis  of  the  SIMS  data  has  not  yet  been 
completed,  but  the  surface  concentrations  of  water  appear  to  be  quite  low 
after  the  exposure  to  humid  gases  (or  dilute  water-methanol  mixtures  -  see 
below);  thus,  the  oxygen  associated  with  the  molecular  water  observed  in 
SIMS  is  probably  at,  or  below,  the  detection  limit  in  AES  and  ISS(~ 100  ppm) . 


Figure  7.  The  in-depth  concentration  profiles  observed  after  an  elevated 

temperature  (200C)  exposure  of  a  f luorozirconate  glass  to  argon- 
which  was  first  saturated  with  water  at  30C;  profiles  were 
obtained  using  secondary  ion  mass  spectroscopy. 


SIMS  was  also  used  Co  examine  ZB  LA  surfaces  after  exposure  Co  liquid 
environments .  Figure  8  summarizes  Che  hydrogen  profiles  obtained  after 
ten  minute  exposures  at  30C  to  acid,  water  and  base.  The  absence  of  any 
in-depth  hydration  (or  any  other  in-depth  effect)  in  ammonium  hydroxide  is 
clearly  evident  and  is  consistent  with  dissolution  studies  that  show  the 
relative  inertness  of  these  glasses  to  high  pH  environments.  In  contrast, 
the  acid  solution  -  within  only  ten  minutes  -  generated  a  visible  film  on 
the  surface,  and  the  hydrogen  profile  in  Figure  8  indicates  that  this 
layer  is  about  4  um  thick.  The  uniform  hydrogen  content,  and  the 
corresponding  profiles  for  other  constituents  of  the  glass,  reveal  that 
this  is  a  layer  or  film  of  surface  reaction  products.  The  scanning 
electron  microscope  (SEM)  was  also  used  to  examine  this  surface;  a  layer 
containing  dense  agglomerates  and  clumps  of  fluoride  crystals  were 
observed.  The  glass  which  was  exposed  to  water  also  exhibited  a  visible 
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Figure  8.  The  in-depth  concentration  profiles  for  hydrogen  in  the 

surface  layer  of  f luorozirconate  glasses  after  exposure  to 
water,  acid  and  base;  the  profiles  were  obtained  with 
secondary  ion  mass  spectroscopy. 
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Figure  9.  The  in-depth  concentration  profiles  for  a  f luorozirconate 

glass  after  exposure  to  water-methanol  mixtures;  the  profiles 
were  obtained  with  secondary  ion  mass  spectroscopy. 


film  on  Che  surface.  The  corresponding  hydrogen-profile  in  Figure  8  shows 
chat  this  layer  is  about  .5-1.0  ym  in  thickness.  Although  the  profile 
suggests  that  the  hydrogen  content  is  greater  after  the  water  corrosion, 
than  after  the  acid  exposure,  this  is  more  likely  due  to  a  density 
difference  between  the  films  produced  in  water  versus  acid.  The 
'acid-generated'  surface  layer  is  loosely  adherent  and  often  flakes-off 
during  drying  of  the  sample,  whereas  the  'water-generated'  surface  layer 
is  usually  more  tenacious. 

Since  the  mechanism  of  aqueous  attack  -  rather  than  the  formation  of 
surface  layers  -  is  of  primary  relevance  to  stress  corrosion,  it  again 
became  necessary  to  lower  the  water  activity  by  dilution  with  methanol  so 
that  the  earliest  stages  of  the  reaction  could  be  investigated.  Figure  9 
shows  two  sets  of  SIMS  depth  profiles;  one  set  was  obtained  after  a  4  hour 
hydration  of  ZBLA  in  50%  HjO  -  50%  CH^OH  and  the  other  after  a  4  hour 
hydration  in  lOZl^O-SOZCH^OH.  Although  other  mixtures  and  conditions 
were  examined,  these  two  samples  bracket  the  critical  range  of  water 
activity.  In  addition  to  profiling  the  hydrated  surfaces,  the  solutions 
were  analyzed  spectrochemically  for  dissolved  constituents  of  the  glass 
after  the  hydration. 

The  profiles  in  Figure  9  look  very  much  like  those  due  to  water  vapor 
hydration  at  elevated  temperatures  (Figure  7).  Namely,  there  has  been 
little  or  no  attack  of  the  glass  structure  -  only  penetration  of  hydrogen 
and  oxygen;  the  corresponding  infrared  spectra  obtained  here,  as  well  as 
those  reported  elsewhere  (13),  indicate  that  these  species  are  associated 
with  molecular  water.  The  solution  analyses  showed  negligible  quantities 
of  dissolved  fluorine  in  the  50%  HjO  -  50%  CH^OH  mixture.  This 
situation  changes  dramatically  once  the  water  concentration  reaches  70% 

(at  least  for  a  4  hour  exposure).  In  this  case,  a  twenty-fold  increase  in 
dissolved  fluorine  is  observed,  and  the  profiles  show  the  initial  stages 
of  a  distinct  'surface  layer'  formation  due  to  this  dissolution  of  the 
glass  structure.  This  layer  contains  substantial  levels  of  hydrogen  and 
fluorine,  but  the  oxygen  profile  is  still  reminiscent  of  an  underlying 
diffusion  process. 

It  seems  clear  that  the  ISS  and  AES  analyses  of  water  or  acid  corroded 
glasses  detected  surface  reaction  products,  but  these  were  fluorides 
rather  than  oxyf luor ides ,  oxides  or  hydroxides.  They  are  probably 
precursors  to  the  ZrF^  and  BaZrF^  crystals  observed  at  the  surface  after 
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more  extensive  corrosion  in  water  and  acid  (12).  The  SIMS  analyses  show 
that  these  fluoride  reaction  products  will  form  layers  at  the  surface 
whose  rate  of  formation  increases  with  decreasing  pH.  In  addition, 
though,  the  SIMS  mass  spectra  (not  presented  here)  show  peaks  at  m/e  *  20, 
40  and  60  whenever  the  surface  has  been  exposed  to  gaseous  or  liquid  water 

environments.  These  peaks  can  only  be  attributed  to  HF+,  (HF)*  and 

+  #  ^ 

( HF ) ^ ,  and  so  it  is  likely  that  these  surface  layers  contain  hydrogen 

fluoride  (hydrogen  cannot  be  detected  with  AES  or  ISS).  The  important 

point,  though,  is  that  these  surface  layers  are  not  relics  of  the  original 

glass  structure,  but  rather,  are  the  reaction  products  of  its  hydrolytic 

decomposition,  dissolution  and  local  precipitation.  Due  to  the 

non-uniform  nature  of  these  surface  layers,  it  is  not  possible  to 

sputter-profile  through  to  the  underlying  'glass'  and  clearly  define  the 

interaction  that  occurs  at  the  interface.  But,  if  the  layer  formation  is 

limited  by  lowering  the  water  activity  (e.g.,  in  the  water  vapor  or  water 

-  methanol  mixtures),  the  penetration  of  molecular  water  into  the  glass 

can  be  seen. 

These  data  show  that  in  the  case  of  both  water  vapor  (at  200-300C)  and 
liquid  water  the  in-depth  penetration  of  molecular  water  is  more  prevalent 
than  hydrolysis  or  chemisorption  at  the  glass  surface.  Tregoat ,  et  al 
(25)  exposed  a  barium-ytterbium-zinc-thorium  fluoride  glass  to  water  vapor 
at  344°C,  and  in  contrast  to  the  behavior  of  this  ZB LA  glass,  they 
observed  HF  gas  evolution  due  to  this  interaction.  Their  infrared 
spectroscopic  analyses  showed  the  absence  of  molecular  water  vibrations  in 
the  hydrated  surface,  and  thus,  they  proposed  an  interdiffusion  reaction 
between  F  and  OH.  Loehr,  et  al .  (13),  also  using  infrared  spectroscopy, 
did  find  molecular  water  -  and  little  or  no  hydroxide  species  -  in  the 
surfaces  of  glasses  hydrated  in  liquid  water.  The  origin  of  these 
apparent  discrepancies  is  not  known,  but  may  be  related  to  temperature  and 
composition  differences. 

5.  DISCUSSION 

Altogether,  the  surface  studies  indicate  that  the  penetration  and  bulk 
diffusion  of  molecular  water  precedes  any  chemisorption,  hydrolysis  or 
breakdown  of  the  glass  structure.  This  appears  to  be  true  in  humid 
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atmospheres  as  well  as  in  liquid  environments.  Once  the  water  reaches 
some  critical  concentration  in  the  glass  surface,  the  glass  structure 
begins  to  disintegrate.  It  is  likely  that  the  molecular  water  -  being 
somewhat  nucleophilic  -  forms  coordinate  bonds  with  the  Ba  and  Zr  cations. 
This  weakens  the  coulombic  forces  of  attraction  with  the  fluorine  anions, 
and  thereby,  'loosens'  the  structure.  In  the  case  of  pure  zirconium 
fluoride,  mono-hydrates  and  tri-hydrates  can  be  generated  by  exposure  of 
ZrF,  (  p  a  4.65  g/cc)  to  steam  at  100°C-450°C.  These  hydrates  are  lower 
in  density  than  the  original  ZrF^  (p*  3.54  g/cc  and  2.81  g/cc, 
respectively),  and  perhaps  of  greater  significance,  they  form  with  little 
or  no  HF  evolution  (26).  Since  there  is  little  or  no  evidence  for 
ion-exchange  or  preferential  leaching  of  any  components  in  the 
fluorozirconate  glasses,  it  seems  certain  that  in-depth  hydration  by 
molecular  water  is  fundamental  to  their  hydrolytic  decomposition  and 
dissolution;  i.e.  it  is  the  first  step  in  the  corrosion  reaction. 

In  H^O-CH^OH  mixtures  and  humid  atmospheres,  the  water  activity  at 
the  surface  is  low,  and  thus,  the  rate  at  which  molecular  water 
concentrates  in  the  surface  is  limited  by  its  diffusion  into  the  bulk  of 
the  glass.  Here,  the  penetration  of  water  occurs  with  little  or  no 
hydrolysis,  dissolution  or  surface  layer  formation  (see  Figures  7  and  9). 
In  liquid  water,  the  water  activity  at  the  glass/solution  boundary  is 
high.  Thus,  breakdown  and  dissolution  of  the  glass  commences  rapidly,  and 
as  a  result,  the  reaction  products  can  accumulate  and  form  surface  layers. 
The  pH  determines  the  solubility  of  the  glass  constituents  -  once  they  are 
released  from  the  glass  structure  -  and  so  pH  is  critical  to  the  overall 
dissolution  rate  and  extent  of  surface  layer  formation.  At  low  pH,  for 
example,  the  dissolution  rate  is  greatly  enhanced  because  the  solubility 
of  fluorine  increases.  The  extensive  surface  layer  in  this  case  is  due  to 
the  insolubiliy  of  ZrF^  at  fluorine  concentration  greater  than  10-lOOppm; 
this  condition  is  easily  achieved  in  acid  near  the  glass  surface,  and  so 
the  local  precipitation  of  ZrF^  is  extensive.  In  contrast,  the 
dissolution  rate  at  pH  10  is  100-1000  times  less  than  at  pH  2  (15),  and 

under  these  conditions  the  penetration  of  water  is  not  observed.  The  fact 
that  oxidation  and/or  hydroxylation  of  ZBLA  was  observed  in  the  surface 
monolayer  only  for  the  pH  11  ammonium  hydroxide  solution  suggests  that 
these  surface  reaction  products  may  'block'  the  penetration  of  water 
and/or  are  exceedingly  insoluble. 


Some  additional  -  albeit  qualitative  -  evidence  for  a  base-catalyzed 
'passivation'  of  the  ZBLA  surface  was  provided  by  examining  the  corrosion 
of  fracture  surfaces  which  were  pre-treated  in  the  base.  The  fracture 
surfaces  were  created  and  aged  in  pH  10  ammonium  hydroxide  solution  and 
were  then  placed  in  water.  A  set  of  control  samples  were  fractured  in 
air,  as  well  as  under  the  water,  for  comparison.  There  was  no  question 
that  the  fracture  surfaces  created  and  aged  at  pH  10  are  more  resistant  to 
visible  surface  film  formation  than  the  glasses  fractured  in  air  or  under 
water;  i.e.,  a  reproducible  'incubation  time'  can  be  observed  before  the 
glasses  treated  at  pHIO  begin  to  haze  in  the  water. 

The  surface  studies  have  provided  considerable  insight  to  the  stress 
corrosion  mechanism  in  ZBLA.  The  fact  that  the  crack-growth  behavior  is 
relatively  insensitive  to  pH  is  consistent  with  the  idea  that  molecular 
water  is  primarily  responsible  for  breakdown  of  the  glass  structure  in 
aqueous  environments.  It  is  also  clear  that  the  crack  velocities  in  acid, 
base,  water  and  water  vapor,  and  water-methanol  mixtures  depend  in  a 
consistent  way  upon  the  water  activity.  But  since  the  extent  of  surface 
film  formation  and  dissolution  varies  considerably  amongst  these 
environments,  it  can  be  concluded  that  surface  film  formation  and/or 
dissolution  are  not  important  features  of  the  stress  corrosion  reaction. 

In  the  case  of  silicate  glasses,  the  mechanism  of  stress  corrosion  has 
been  attributed  to  a  stress-enhanced  dissociative  chemisorption  at  the 
crack  tip  (7).  This  interaction  occurs  at  the  crack-tip  between  strained 
bonds  and  chemical  species  in  the  environment.  The  theory  and 
experimental  data  support  the  idea  that  chemical  species  which  possess 
both  electron  donor  and  proton  donor  sites  are  a  prerequisite  for  this 
reaction.  This  property  of  chemical  species  has  been  characterized  by 
Gutman  (27)  in  terms  of  so-called  donor  numbers  and  acceptor  numbers. 

Thus,  chemical  environments  with  high  values  for  both  of  these  numbers 

would  be  expected  to  enhance  the  crack-growth  rate  to  the  greatest  extent. 

•  1/2 

In  the  case  of  these  fluoride  glasses  the  crack-velocity  at  “.2  MPa-m 
(above  the  threshold  for  aqueous  environments)  shows  an  excellent 
correlation  with  the  sum  of  those  two  numbers  in  all  of  the  environments 
(28).  This  suggests  that  dissociative  chemisorption  may  also  be 
responsible  for  stress  corrosion  in  these  fluoride  glasses.  However,  the 
surface  studies  of  adsorption  on  clean  fracture  surfaces  provided  little 
evidence  for  the  chemisorption  of  water.  This  indicates  that  the  presence 
of  strained-bonds  may  be  required  for  the  chemisorption  reaction  in  this 


system.  Although  the  ISS  data  in  Figure  4  supports  this  idea,  many  of  the 
other  surfaces  prepared  by  fracture  in  gaseous  and  liquid  environments  did 
not  show  chemisorbed  oxide  a  hydroxide  species  at  the  surface.  It  is  to 
be  emphasized,  though,  that  these  were  all  'fast-fractured'  in  the 
environment,  and  so,  these  surfaces  may  not  be  representative  of  those 
prepared  under  sub-critical  crack  growth  conditions.  Clearly,  the  in-situ 
analysis  of  fracture  surfaces  created  at  low  velocities  would  be  a 
worthwhile  endeavor. 

On  the  other  hand,  the  surface  studies  have  indicated  that  the 
penetration  of  molecular  water  is  the  most  prevalent  interaction  between 
aqueous  environments  and  the  f luorozirconate  glasses.  If  a  weakening  of 
the  ionic  bonds  due  to  the  coordination  of  molecular  water  with  the  Zr  and 
Ba  cations  is  responsible  for  the  sub-critical  crack  growth,  then  the 
crack  velocity  should  scale  with  the  dielectric  constant  or  polarizibility 

of  the  environment.  The  relationship  between  crack  velocity  at  K  *  .2 
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MPa-m  ,  and  the  dielectric  constant  of  the  various  environments,  does 

not  support  this  concept  though.  Thus,  while  hydration  by  molecular  water 

may  be  fundamental  to  the  corrosion  and  dissolution  reactions,  it  may  not 

be  responsible  for  the  bond-breaking  event  in  stress  corrosion. 

Of  particular  significance  are  the  thresholds  (or 
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stress-corrosion  limits)  -  at  about  .20  MPa-m  -  which  are  exhibited  by 
the  HjO  environments  and  the  NH^.  These  environments  also  exhibit 
velocity  plateaus  whose  magnitude  scales  with  the  water  activity.  Since 
the  rapid  penetration  and  diffusion  of  molecular  water  is  prevalent  in 
these  environments,  one  is  led  to  consider  the  possible  effects  of 
'stress-enhanced'  water  diffusion  ahead  of  the  crack-tip.  This  phenomena 
could  influence  the  stress  distribution  at  the  crack-tip  and/or  modify  the 
glass  through  which  the  crack-tip  must  propagate.  Thus,  the  stress 
intensity  and/or  activity  of  water  would  determine  the  extent  of  water 
diffusion  ahead  of  the  crack-tip,  and  quite  possibly,  this  water  of 
hydration  might  dissociate  -  under  stress  -  as  the  crack  propagates.  And 
in  a  qualitative  sense,  an  interaction  between  the  glass  and  environment 
'ahead  of  the  crack-tip'  would  explain  the  erratic  behavior  of  the  cracks 
in  these  glasses.  In  contrast  to  the  behavior  of  silicate  glasses  during 
crack-growth  experiments,  the  cracks  in  even  well-annealed  fluoride 
glasses  seldom  follow  a  path  which  is  perfectly  orthogonal  to  the  applied 
moment.  It  is  also  found  that  once  the  crack  has  arrested  -  e.g.  in  and 
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below  the  threshold  region  of  stress  intensity  -  the  crack  will  not 
re-propagate  when  the  load  is  increased  unless  the  specimen  is  'tapped*  to 
get  it  running. 

6.  SUMMARY 

A  detailed  interpretation  of  the  stress  corrosion  behavior  in 
f luorozirconate  glasses  is  only  just  beginning  to  emerge.  The  origin  of 
the  stress-intensity  thresholds  and  velocity  plateaus  is  not  obvious,  but 
clearly,  an  understanding  of  these  features  could  contribute  to  the 
development  of  a  model  for  stress  corrosion  in  these  materials.  In  this 
regard,  it  is  noteworthy  that  the  environments  where  velocity  thresholds 
and  plateaus  are  observed  also  happen  to  be  those  where  the  potential 
formation  of  HF  reaction  products  can  be  proposed  .  Thus,  the  dependence 
of  crack  growth  rates  upon  the  ionic  strength  of  HF-^O  mixtures  would 
be  useful  information.  In  general,  though,  a  more  substantial  data 
base  concerning  the  surface  chemistry  and  slow-crack  growth  character¬ 
istics  of  these  glasses  will  be  required,  before  a  definitive  model  can 
be  proposed  and  validated.  This  understanding  will  be  necessary  when 
evaluating  the  fatigue  data  yet  to  be  reported  for  fluorozirconate 
optical  fibers,  and  will  be  of  great  benefit  in  the  design  and  develop¬ 
ment  of  more  fatigue-resistant  systems.  And  since  these  fluorozirconate 
materials  represent  a  new  and  unique  class  of  ionical ly-bonded  glasses, 
the  ability  to  compare  their  slow  crack  growth  behavior  with  that  of 
the  more  covalently-bonded  silicate  glasses  could  provide  a  more  funda¬ 
mental  insight  to  stress-corrosion  reactions  in  brittle  solids. 
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ABSTRACT 


The  elastic  moduli  of  several  organometallic  silica  gels  were  measured 
using  pulse-superposition  interferometry  (PSP).  A  series  of  gels  were  prepared 
under  acid-catalyzed  conditions  where  the  molar  ratio  of  water  to 
tetraethoxysilane  in  the  sol  was  varied  (4:1,  12:1,  16:1  and  24.1).  It  was 
found  that  the  elastic  moduli  varied  systematically  with  the  hydrolysis 
conditions  used  to  prepare  the  gels,  as  well  as  with  the  relative  humidity  used 
to  equilibrate  the  gels  before  and  during  the  PSP  measurement.  The  experimental 
PSP  data  for  the  gels  were  evaluated  using  a  composite  model  (matrix  plus  pores) 
and  then  the  elastic  moduli  of  the  gel  matrix  were  determined  by  applying  a 
correction  for  porosity  based  on  a  self  consistent  scheme  (SCS)  approximation. 
The  results  show  that  the  elastic  moduli  of  the  gel  matrix  also  exhibit  a 
systematic  dependence  upon  the  preparation  conditions,  but  they  always  fall 
below  the  corresponding  values  for  fused  silica.  Even  when  excess  water  was 
used  to  hydrolyze  the  TEOS,  this  did  not  lead  to  a  fully  polymerized  silica 
matrix  structure.  It  is  suggested  that  the  lowered  matrix  moduli  values  are  due 
to  the  presence  of  residual  hydroxy  and  ethoxy  groups  which  are  internal  to  a 
polymeric  structure. 


INTRODUCTION 


The  sol/gel  method  of  glass  formation  is  limited  to  the  production  of  only 
small  bodies  because  of  the  cracking  and  fragmentation  which  occurs  during 
processing.  This  is  especially  evident  in  the  case  of  organometallic  sol/gel 
materials;  that  is,  silica  gels  prepared  by  the  hydrolysis  of  tetralkoxysi lanes 
in  non-aqueous  solvents.  It  is  well  recognized  that  the  mechanical  stresses 
generated  during  drying  and  firing  are  responsible  for  the  fracture,  but  there 
have  been  only  a  few  attempts  to  model  them  in  terms  of  the  capillary  forces, 
density  gradients,  and  their  dependence  upon  the  rate  of  drying  or  firing.  The 
development  of  such  models  requires  information  concerning  the  elastic  modulus 
of  the  gel  as  a  function  of  the  textural  parameters  (i.e.,  the  matrix  structure 
and  pore-size  distribution)  and  the  relative  humidity  (i.e.,  the  effects  of 
adsorbed/condensed  water).  And  independent  of  their  role  in  determining  the 
mechanical  integrity  of  these  materials  during  drying  and  firing,  the  structure 
and  mechanical  properties  of  the  'gel  matrix'  are  of  fundamental  interest. 

This  paper  reports  a  systematic  study  concerning  the  effects  of  the 
solution  composition  upon  the  elastic  moduli  of  silica  gels;  the  dependence  of 
the  elastic  moduli  upon  relative  humidity  is  also  presented.  The  structure  of 
these  silicon  alkoxide  gels  consists  of  a  silica-like  matrix  with  microporosity 
whose  volume  fraction,  distribution  and  connectivity  varies  with  the  hydrolysis 
conditions.  The  elastic  moduli  of  the  matrix  phase  were  evaluated  by  modeling 
the  composite  elastic  moduli  in  terras  of  the  volume  fraction  of  the  gel  matrix 
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EXPERIMENTAL  PROCEDURE 

Gel  Preparation 

A  series  of  silicon  alkoxide  gels  (with  water:  TEOS  of  4:1,  12:1,  16:1,  and 
24:1)  were  prepared  systematically  (see  Table  1)  using  TEOS  (SHOCjH^)^) , 
ethanol  (C2H5OH),  water  and  a  1M  hydrochloric  acid  solution  (1).  These 
HjO-to-TEOS  ratios  represent  1,  3,  4,  and  6  times  the  theoretical  quantity  of 
water  required  for  complete  hydrolysis,  respectively,  and  are  hereinafter 
referred  to  as  IX,  3X,  4X  and  6X.  The  solutions  were  cast  in  16  mm  x  150  mm 
pyrex  test  tubes  and  left  to  gel  at  room  temperature.  The  gelation  and  drying 
were  carried  out  over  a  period  of  several  months. 


K- 
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Density 

The  bulk  density  of  the  gels  was  determined  by  Archimedes  and  dimensional 
methods.  The  bulk  densities  measured  by  the  Archimedes  method  were  conducted  in 
water  using  a  wax  immersion  technique  approved  by  ASTM  (2).  Five  determinations 
were  made  on  each  sample.  The  results  are  reported  as  an  average  value  with  95% 
confidence  intervals.  The  bulk  densities  measured  by  dimensional  methods  were 
carried  out  on  cylindrical  gel  samples  with  the  end  faces  polished  flat  and 
parallel  with  one  micron  diamond  paste.  The  diameter  and  thickness  were 
measured  using  a  hand  held  micrometer,  accurate  to  plus  or  minus  5  microns. 
Twenty  determinations  were  made  per  sample  and  reported  as  an  average  with  95% 
confidence  intervals. 

The  matrix  density  was  determined  by  further  drying  the  gels  at  0%  relative 

humidity.  The  dry  weight  was  recorded  and  then  the  gels  were  saturated  with 

deionized  water.  The  volume  of  the  matrix  was  determined  in  the  following 

manner.  The  volume  of  the  water  in  the  pores,  Vu 

H, 


,  is  given  by 
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where  V  is  assumed  equivalent  to  the  volume  of  the  pores  in  the  gel 

h2o 

samples,  A  is  the  saturated  mass,  A  is  the  mass  of  the  matrix,  and  p  n  is 
s  m 

the  density  of  water.  Then, 


where  V  is  the  volume  of  the  matrix  and  V.  is  the  volume  of  the  bulk  gel. 
m  D 


Elastic  Moduli 

The  elastic  moduli  were  determined  using  an  acoustic  interferometry 
technique  known  as  ultrasonic  pulse-superposition  (PSP),  as  described  by 
McSkimin  (3).  PSP  utilizes  an  applied  repetitive  sequence  of  radio  frequency 
pulses-in  the  range  KHz  to  GHz-transmitted  through  a  piezoelectric  transducer 
which  is  bonded  to  the  sample.  The  time  interval  during  which  the  elastic  wave 
propagates  through  the  sample,  from  one  face  to  the  other  and  back,  may  be 
adjusted  to  equal  precisely  the  repitition  rate  of  the  next  RF  pulse.  When  this 
is  achieved,  the  applied  pulses  are  superimposed  upon  the  return  echoes 
producing  a  constructive  interference.  The  echoes  then  build  up  sufficient 
amplitude  to  be  observed  on  an  oscilloscope. 

The  samples  were  prepared  as  flat  parallel  cylinders  5-10  mm  in  diameter 
and  1-3  mm  thick.  X-cut  and  AC-cut  ( 3  mm  and  4.6  mm)  quartz  transducers  were 
used  to  generate  and  determine  the  travel-time  of  the  longitudinal 
(compressional)  and  transverse  (shear)  waves,  respectively.  The  instrument  used 
in  this  work  had  a  frequency  range  of  15-25  MHz,  but  all  of  the  gel  measurements 
were  made  at  20  MHz. 

The  Young's  elastic  modulus  was  determined  from  the  bulk  density  (p)  of  the 


sample  and  Che  corapressional  (V  )  and  shear  wave  (Ve)  velocities  by 
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The  shear  modulus  was  determined  from  the  density  of  the  samples  and  the  shear 
wave  velocities  by 


G  *  pV 


The  bulk  modulus  (K)  was  determined  via 


K  *  PV  2  -  4/  V  2 
H  p  3  s 


and  the  Poisson's  ratio  (v)  was  determined  from  the  ratio  of  the  Young's  elastic 
modulus  and  the  shear  modulus  using 
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Controlled  Humidity  Measurements 


Young's  elastic  modulus,  the  shear  modulus,  the  bulk  modulus,  and  Poisson's 
ratio  were  also  determined  as  a  function  of  the  relative  humidity  at  room 
temperature  using  the  PSP  method.  Here,  the  relative  humidity  was  controlled  by 
mixing  dry  and  saturated  nitrogen  gases  in  a  one  cubic  foot  dessicator.  The 
humidity  levels  were  recorded  by  a  Bacharach  hygrometer  which  is  accurate  to 
plus  or  minus  two  percent.  The  cylindrical  gel  samples  were  equilibrated  at 
increments  of  ten  percent  relative  humidity,  and  although  not  reported,  the 
kinetics  of  equilibration  were  simultaneously  monitored  by  the  PSP  travel-time 
measurements  at  five  minute  intervals  up  to  one  hour. 


Model  for  Composite  Moduli 

A  self  consistency  scheme  (SCS)  was  developed  by  Budiansky,  and  also  by 
Hill  (4,5),  to  characterize  the  elastic  properties  of  isotropic  multi-phase 
composites.  This  model  was  used  here  to  generate  theoretical  elastic  moduli 
values  for  the  saturated  and  dry  gel  composites  assuming  a  'fused'  silica 
matrix.  These  values  were  then  compared  with  the  experimentally  determined 
properties.  The  effective,  or  composite,  moduli  (M  )  are  given  by: 
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where  M  is  either  the  bulk  modulus  (K  )  or  the  shear  modulus  (G  )  for  the 

composite,  M.  is  either  K,  or  G,  for  the  pore++,  M  is  either  K  or  G 
i  xi  m  id  in 

for  the  gel  matrix,  and  n  is  the  volume  fraction  porosity.  F  depends  on  whether 
M  is  K  or  G  ,  and  the  pore  geometry,  as  follows; 


Spherical  Pores 


For  K 


For  G 
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Needle  Pores 


For  K 


F  =  -i  (G,  +  G*) 


In  Che  present  case,  the  pores  are  either  dry  or  saturated  with  water 
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where 
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Due  to  the  nature  of  the  foregoing  relations, K  and  G  must  be  evaluated  by  an 
iterative  procedure  (6). 

The  Young's  elastic  modulus,  shear  modulus,  bulk  modulus,  and  Poisson's 

ratio  of  the  gel  matrix  were  also  evaluated  using  the  SCS  model  since  the 

abovementioned  calculations  showed  the  assumption  of  a  'fused  silica'  matrix  to 

be  false.  The  experimental  values  of  the  shear  and  bulk  modulus  of  the  gel 

composites  were  used  in  this  calculation  and  the  pore  geometry  was  assumed  to  be 

needle-shaped.  It  is  reasonable  to  believe  that  needle-shaped  porosity  is 

meaningful  for  these  gel  materials,  and  indeed,  it  is  shown  that  the  theoretical 

predictions  of  the  composite  properties  where  needle-shapes  are  assumed  more 

closely  match  the  experimental  composite  values.  The  shear  modulus  (G  )  and 

m 

bulk  modulus  (K  )  of  the  gel  matrix  were  determined  using 
tn 


M 

m 


* 

M  F 

F(l-n)  -  M  n 


(13) 


which  follows  directly  from  equation  7  using,  in  this  case,  Kj=G^=0.  The 
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experimental  values  for  M  (i.e.  K  or  G  )  were  used,  and  for  K 
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and  for  G 
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where  g  is  given  by  equation  12.  The  Young's  elastic  modulus  (Em)  and  the 
Poisson's  (  )  ratio  of  the  gel  matrix  were  then  determined  using  the  relations 
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RESULTS  AND  DISCUSSION 

Density 

The  bulk  densities  for  the  IX,  3X,  4X,  and  6X  acid  catalyzed  TEOS  gels  were 
determined  by  the  Archimedes'  method  after  equilibration  at  a  constant  humidity. 
It  is  important  to  emphasize  that  there  was  little  correlation  between  the  gel 
compositions,  the  measured  densities  and  the  measured  wave  velocities  until  the 
gels  were  equilibrated  at  a  controlled  humidity  before  or  during  any 
measurement.  The  attainment  of  equilibrium  could  be  verified  by  measuring  the 
wave  velocities  during  this  equilibration. 

The  average  bulk  densities  as  a  function  of  relative  humidity,  the  average 
matrix  density,  and  the  corresponding  porosity,  are  presented  for  each  gel 
composition  in  Table  2.  The  bulk  densities  were  used  to  determine  the  elastic 
moduli  of  the  gel  composite  as  a  function  of  relative  humidity  while  the  matrix 
densities  were  used  in  the  evaluation  of  the  matrix  elastic  moduli.  The  density 
measured  at  0%  relative  humidity  was  assumed  to  be  the  actual  bulk  density  of 


the  dry  gel. 


Table  2.  The  average  bulk  density  as  a  function  of  relative  humidity,  the  matrix  density 
and  the  percent  porosity  for  IX,  3X,  4X,  and  6X  gel  specimens. 
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Figures  1  and  2  illustrate  the  bulk  and  matrix  density  (at  0%  relative 
humidity)  as  a  function  of  the  molar  ratio  of  water  to  TEOS  used  in  the  gel 
processing.  The  matrix  density  is  linearly  proportional  to  the  molar  ratio  of 
water  to  TEOS  whereas  the  bulk  density  shows  an  inverse  relationship.  This  type 
of  behavior  has  been  reported  by  others  (7-9).  This  trend  has  been  explained  in 
terms  of  the  relative  rates  of  hydrolysis  and  condensation  for  different  ratios 
of  to  TEOS  and  the  corresponding  structural  units  so-produced. 

In  an  acid  catalyzed  silicon  alkoxide  gel,  hydrolysis  proceeds  via 
electrophilic  attack  (8).  At  low  molar  ratios  of  water  to  TEOS  (with  a  constant 
acid  concentration),  highly-extended  polymers  can  be  expected  (1)  The  backbone 
of  this  polymer  is  based  upon  Si-O-Si  bonding,  but  the  'chain'  is  framed  with 
residual  ethoxy  and  hydroxy  groups.  Upon  drying  the  gel,  the  highly-extended, 
but  weakly  cross-linked,  gel  shrinks  upon  itself  with  little  resistance.  Thus, 
these  gels  dry  to  a  relatively  high  bulk  density,  but  the  matrix  density  is  low 
due  to  the  high  residual  organic  component. 

Conversely,  at  high  molar  ratios  of  water  to  TEOS  (again  with  a  constant 
acid  concentration),  the  polymer  units  in  the  gel  contain  a  lower  percentage  of 
residual  ethoxy  groups,  and  correspondingly,  are  denser  and  probably  more  highly 
cross-linked.  During  the  desiccation  and  drying  of  these  gels,  the  more  highly 
cross-linked  polymer  units  (which  can  be  particle-like  at  very  high  l^O/TEOS) 
resist  the  shrinkage  due  to  packing  contraints  and/or  due  to  additional 
polymerization  as  the  system  concentrates.  In  either  case,  the  result  is  a  gel 
of  lower  bulk  density  but  higher  matrix  density. 

Gel  Moduli  and  Their  Humidity  Dependence 


Table  3  presents  the  average  elastic  wave  velocities  (compressional  and 
shear)  as  a  function  of  the  relative  humidity  and  the  sol/gel  composition.  The 
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The  variation  in  the  average  matrix  density  for 
IX,  3X,  AX,  and  6X  acid  catalyzed  TEOS  gels;  the 
(o)  data  point  is  an  independent  density  measure 
ment  for  a  AX  gel  from  reference  (12). 
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ties  (km/sec)  as  a  function  of  relative  humidity 
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elastic  wave  velocities  listed  in  Table  3,  and  the  densities  given  in  Table  2, 
were  used  to  determine  the  Young's  elastic  modulus  (E),  bulk  modulus  (K),  shear 
modulus  (G) ,  and  Poisson’s  ratio.  The  values  for  E,  K,  and  G,  as  functions  of 
relative  humidity  (P/P^)  for  the  IX,  3X,  4X,  and  6X  gels,  are  shown  in  Figures 
3  and  4.  The  elastic  moduli  and  Poisson's  ratios  for  the  'dry'  gels  are 
summarized  in  Table  4. 

In  general,  the  gels  exhibit  the  expected  trend  of  increasing  elastic 

moduli  with  increasing  density.  The  slopes  of  the  E  and  X  vs  P/P  plots  also 

o 

decrease  with  increasing  density,  presumably  due  to  the  decrease  in  pore  volume 
with  increasing  density.  The  non-linearity  in  the  E  and  K  vs  P/PQ  curves  at 
low  relative  humidity  is  probably  due  to  an  adsorbed  water  effect;  and  the  fact 
that  this  inflection  is  enhanced  for  the  higher  porosity  gels  is  consistent  with 
a  dependence  upon  surface  area.  The  shear  modulus  shows  a  weaker  dependence  on 
relative  humidity  because  of  the  negligible  rigidity  of  the  fluid  phase. 


Matrix  Moduli  as  a  Function  of  the  H..O/TEOS  Ratio 


Initially,  the  elastic  moduli  of  the  gel  matrices  were  assumed  to  have  the 
value  for  fused  silica,  and  then  tl i  SCS  model  was  utilized  to  generate 
theoretical  wave  velocities  and  moduli  for  'fused  silica’  gel/composites.  The 
compressional  wave  velocities  for  these  theoretical  'fused  silica' 
gel/composites  are  compared  to  the  experimentally  measured  velocities  of  the 
bulk  gels  in  Figure  5.  The  compressional  wave  velocity  (and  the  elastic  moduli) 
increase  with  density,  in  general,  but  the  dependence  upon  density  is 
quantitatively  distinct  for  the  'fused  silica'  gel/composites  versus  the 
' organometal 1 ic  si lica ' -gel/compos ices.  No  manipulation  of  the  model  or  pore 
geomety  could  be  found  which  would  reduce  the  compressional  wave  velocities  of 


Che  'fused  silica'-gel  to  match  (or  approximate)  chose  for  the  ' organometal lie 
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Table  4.  The  average  elastic  moduli  (GPa)  and  Poisson's  ratio  for 
IX,  3X,  4X  and  6X  gel  specimens  at  0%  relative  humidity. 


Gel  Young's  Elastic  Shear  Bulk  Poisson's 

Modulus*  Modulus*  Modulus*  Ratio* 


Compressional 
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si  lies' -gel.  One  can  only  conclude  that  Che  matrix  phase  in  silica  gels 
prepared  with  TEOS  does  not  exhibit  the  elastic  properties  of  fused  silica. 

Figure  6  presents  the  average  matrix  elastic  moduli  (calculated  according 
to  equations  13-17)  as  a  function  of  the  molar  ratio  of  H^O  to  TEOS  for  the  IX, 
3X,  4X,  and  6X  gels.  It  is  clear  that  increasing  the  HjO/TEOS  ratio  yields  a 
matrix  phase  of  increasing  stiffness.  This  trend  is  consistent  with  the 
aforementioned  polymer  models  for  acid  catalyzed  TEOS  gels  (1,  7-10)  which 
suggest  the  formation  of  more  highly  condensed  silicates  with  increasing 
HjO/TEOS  ratios.  But  these  data  show  clearly  that  although  the  elastic  moduli 
of  the  matrix  phase  approach  the  values  for  fused  silica,  even  the  6X  gel  falls 
short  by  a  considerable  amount.  In  simple  terms,  this  means  that  the  solid 
material  which  exists  after  equilibration  of  the  gels  at  0Z  relative  humidity  is 
not  'silica'.  However,  it  is  not  possible  to  define  a  more  precise  structural 
interpretation  with  these  data  alone.  The  texture  of  the  matrix,  the 
concentration  of  residual  ethoxy  groups,  and  the  extent  of  hydration  all  combine 
to  determine  the  matrix  density  and  moduli. 


SUMMARY  AND  CONCLUSION 

A  series  of  silicate  sol/gel  materials  were  fabricated  with  varying  ratios 
of  water  to  TEOS.  Under  controlled  atmosphere  conditions,  there  is  a  clear 
trend  in  the  density  and  elastic  moduli  of  the  gel  composites  as  a  function  of 
percent  relative  humidity.  The  gels  were  treated  as  two-phase  'gel  composites' 
wherein  the  solid  phase  is  silica,  or  silica-like,  while  the  other  phase  is  a 
f luid-namely ,  dry  air,  moist  air  or  water.  Further  characterization  of  the 
solid  matrix-phase  was  then  carried-out  using  a  self-consistent  scheme 
approximation  (SCS)  for  interpretation  of  the  elastic  properties. 

The  matrix  density  for  each  gel  composition  was  determined  using  the  dry 


Matrix  Elastic  Moduli  (GPa) 
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gel  composites;  i.e.,  gels  equilibrated  at  OZ  relative  humidity.  The  matrix 
density  ranged  from  1.88  gem  ^  to  2.19  gem  ^  for  4:1  to  24:1  gels  (these  may 

_3 

be  compared  to  fused  silica  whose  density  is  2.2007  gem  ).  A  theoretical 
correction  for  porosity,  based  on  the  self  consistent  scheme  (SCS) 
approx imat ion,  was  applied  to  the  gel  composites  to  evaluate  the  elastic  moduli 
of  the  gel  matrix.  The  matrix  values  of  the  Young's  elastic  modulus  varied  from 
13.8  GPa  to  48.8  GPa,  the  shear  modulus  from  5.5  GPa  to  21.5  GPa,  the  bulk 
modulus  from  9.3  GPa  to  23.3  GPa,  and  the  Poisson's  ratio  from  0.24  to  0.16  over 
the  range  of  compositions  studied;  these  may  be  compared  to  the  values  for  fused 
silica  which  are  72.8  GPa,  31.3  GPa,  36.2GPa  and  0,16,  respectively  (11).  It  is 
important  to  emphasize,  though,  that  the  calculation  of  the  matrix  elastic 
moduli  for  the  organometal 1 ic  silica  gel/composites  assumed  a  uniform 
distribution  of  needle-shaped  pores.  However,  it  is  also  noteworthy  that  no 
manipulation  of  the  model  was  sufficient  to  explain  the  organometal lie  silica 
gel/composi te  properties  when  a  fused  silica  matrix  was  assumed. 

The  analysis  of  the  matrix  elastic  moduli  further  supports  the  idea  that  the 
hydrolysis/condensation  of  tetraethoxys i lane  does  not  produce  a  ' fused-si iica ' 
matrix.  The  values  of  the  elastic  moduli  for  the  gel  matrix,  over  the  range  of 
the  compositions  studied,  consistently  fell  below  those  of  fused  silica.  This 
is  presumably  due  to  a  limited  degree  of  condensation  during  the  solut ion-to-gel 
transition.  Thus,  there  may  be  terminal  hydroxy  and  ethoxy  groups  which  not 
only  'frame'  the  matrix  phase,  but  which  are  also  internal  to  it.  Nonetheless, 
the  results  show  a  trend  towards  increasing  stiffness  as  the  water  necessary  for 
hydrolysis  is  increased. 
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